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We investigate the problem of detecting gravitational waves from binaries of nonspinning black 
holes with masses m — 5-20M©, moving on quasicircular orbits, which are arguably the most 
promising sources for first-generation ground-based detectors. We analyze and compare all the cur- 
rently available post-Newtonian approximations for the relativistic two-body dynamics; for these 
binaries, different approximations predict different waveforms. We then construct examples of de- 
tection template families that embed all the approximate models, and that could be used to detect 
the true gravitational- wave signal (but not to characterize accurately its physical parameters). We 
estimate that the fitting factor for our detection families is > 0.95 (corresponding to an event-rate 
\£) , loss < f 5%) and we estimate that the discretization of the template family, for ~ f0 4 templates, 

■ increases the loss to < 20%. 
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1—5 . I. INTRODUCTION 

\o ■ 

A network of broadband ground-based laser interferometers, aimed at detecting gravitational waves (GWs) in the 
I frequency band 10-10 3 Hz, is currently beginning operation and, hopefully, will start the first science runs within this 
J> ■ year (2002). This network consists of the British-German GEO, the American Laser Interferometer Gravitational- wave 
(N| " Observatory (LIGO), the Japanese TAMA and the Italian-French VIRGO (which will begin operating in 2004) Q. 
. The first detection of gravitational waves with LIGO and VIRGO interferometers is likely to come from binary 
black-hole systems where each black hole has a mass |69( of a few M , and the total mass is roughly in the range 
10-40Mq 0, and where the orbit is quasicircular (it is generally assumed that gravitational radiation reaction will 
circularize the orbit by the time the binary is close to the final coalescence 0). It is easy to see why. Assuming 
for simplicity that the GW signal comes from a quadrupole-governed, Newtonian inspiral that ends at a frequency 
** — ■ outside the range of good interferometer sensitivity, the signal-to-noise ratio S/N is cx M 5 / 6 /d (See, e.g., Ref. 3), 
where A4 = Mrj 3 ^ 5 is the chirp mass (with M = mi + ni2 the total mass and rj — mim2/M 2 ), and d is the distance 
I , between the binary and the Earth. Therefore, for a given signal-to-noise detection threshold (see Sec. [nj and for 
equal-mass binaries (rj — 1/4), the larger is the total mass, the larger is the distance d that we are able to probe. [In 
Sec. we shall see how this result is modified when we relax the assumption that the signal ends outside the range 
of good interferometer sensitivity.] 

For example, a black-hole-black hole binary (BBH) of total mass M — 20M© at 100 Mpc gives (roughly) the same 
S/N as a neutron-star-neutron-star binary (BNS) of total mass M = 2.8M Q at 20 Mpc. The expected measured-event 
rate scales as the third power of the probed distance, although of course it depends also on the system's coalescence 
rate per unit volume in the universe. To give some figures, computed using LIGO-I's sensitivity specifications, if we 
assume that BBHs originate from main-sequence binaries |5j, the estimated detection rate per year is < 4 x 10~ 3 -0.6 
at 100 Mpc H0, while if globular clusters are considered as incubators of BBHs the estimated detection rate 
per year is ~ 0.04-0.6 at 100 Mpc 0,0; by contrast, the BNS detection rate per year is in the range 3 x 10 _4 -0.3 
at 20 Mpc [10. The very large cited ranges for the measured-event rates reflect the uncertainty implicit in using 
population-synthesis techniques and extrapolations from the few known galactic BNSs to evaluate the coalescence 
rates of binary systems. [In a recent article [9( , Miller and Hamilton suggest that four-body effects in globular clusters 
might enhance considerably the BBH coalescence rate, brightening the prospects for detection with first-generation 
interferometers; the BBHs involved might have relatively high BH masses (~ 100M Q ) and eccentric orbits, and they 
will not be considered in this paper.] 

The GW signals from standard comparable-mass BBHs with M = 10-40M Q contain only few (50-800) cycles in 
the LIGO-VIRGO frequency band, so we might expect that the task of modeling the signals for the purpose of data 
analysis could be accomplished easily. However, the frequencies of best interferometer sensitivity correspond to GWs 
emitted during the final stages of the inspiral, where the post-Newtonian (PN) expansion [lOj, which for compact 
bodies is essentially an expansion in the characteristic orbital velocity v/c, begins to fail. It follows that these sources 
require a very careful analysis. As the two bodies draw closer, and enter the nonlinear, strong-curvature phase, the 
motion becomes relativistic, and it becomes harder and harder to extract reliable information from the PN series. 
For example, using the Keplerian formula v = (ttM/qw) 1 ^ 3 [where /gw is the GW frequency] and taking /gw = 153 
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Hz [the LIGO-I peak-sensitivity frequency] we get v(M) = 0.14(M/M©) 1 / 3 ; hence, for BNSs v(2.8M ) = 0.2, but for 
BBHs v(20M Q ) = 0.38 and w(40M Q ) = 0.48. 

The final phase of the inspiral (at least when BH spins are negligible) includes the transition from the adiabatic 
inspiral to the plunge, beyond which the motion of the bodies is driven (almost) only by the conservative part of the 
dynamics. Beyond the plunge, the two BHs merge, forming a single rotating BH in a very excited state; this BH 
then eases into its final stationary Kerr state, as the oscillations of its quasinormal modes die out. In this phase the 
gravitational signal will be a superposition of exponentially damped sinusoids (ringdown waveform). For nonspinning 
BBHs, the plunge starts roughly at the innermost stable circular orbit (ISCO) of the BBH. At the ISCO, the GW 
frequency [evaluated in the Schwarzschild test-mass limit as ~ 0.022/M] is /^°(20M Q ) ~ 220 Hz and 

/g S w°( 30M o) ~ 167Hz. These frequencies are well inside the LIGO and VIRGO bands. 

The data analysis of inspiral, merger (or plunge), and ringdown of compact binaries was first investigated by 
Flanagan and Hughes ^l], an d more recently by Damour, Iyer and Sathyaprakash [l^. Flanagan and Hughes 0] 
model the inspiral using the standard quadrupole prediction (see, e.g., Ref. Q), and assume an ending frequency of 
0.02/M (the point where, they argue, PN and numerical-relativity predictions start to deviate by ~ 5% They 
then use a crude argument to estimate upper limits for the total energy radiated in the merger phase (~ 0.1M) and 
in the ringdown phase (~ 0.03M) of maximally-spinning-BBH coalescences. Damour, Iyer and Sathyaprakash [T^ 
study the nonadiabatic PN-resummed model for non spinning BBHs of Refs. E3- EH- where the plunge can be 
seen as a natural continuation of the inspiral [15| rather than a separate phase; the total radiated energy is 0.007M in 
the merger and 0.007M in the ringdown 1171. (All these values for the energy should be also compared with the value, 
0. 025-0. 03M, estimated recently in Ref. 18] for the plunge and ringdown for non spinning BBHs.) When we deal 
with nonadiabatic models, we too shall choose not to separate the various phases. Moreover, because the ringdown 
phase does not give a significant contribution to the signal-to- noise ratio for M < 200M Q [ll], [lj , we shall not include 
it in our investigations. 

BHs could have large spins: various studies 0, |2(j have shown that when this is the case, the time evolution 
of the GW phase and amplitude during the inspiral will be significantly affected by spin-induced modulations and 
irregularities. These effects can become dramatic, if the two BH spins are large and are not aligned or antialigned 
with the orbital angular momentum. There is a considerable chance that the analysis of interferometer data, carried 
out without taking into account spin effects, could miss the signals from spinning BBHs altogether. We shall tackle 
the crucial issue of spin in a separate paper |2l|. 

The purpose of the present paper is to discuss the problem of the failure of the PN expansion during the last 
stages of inspiral for nonspinning BHs, and the possible ways to deal with this failure. This problem is known in 
the literature as the intermediate binary black hole (IBBH) problem [22| . Despite the considerable progress made by 
the numerical-relativity community in recent years jl 'l l23l l24l |25| , a reliable estimate of the waveforms emitted by 
BBHs is still some time ahead (some results for the plunge and ringdown waveforms were obtained very recently 18], 
but they are not very useful for our purposes, because they do not include the last stages of the inspiral before the 
plunge, and their initial data are endowed with large amounts of spurious GWs). To tackle the delicate issue of the 
late orbital evolution of BBHs, various nonperturbative analytical approaches to that evolution (also known as PN 
resummation methods) have been proposed [lj, [15J, [16|, |26( . 

The main features of PN resummation methods can be summarized as follows: (i) they provide an analytic (gauge- 
invariant) resummation of the orbital energy function and gravitational flux function (which, as we shall see in Sec. IIIII 
are the two crucial ingredients to compute the gravitational waveforms in the adiabatic limit); (ii) they can describe 
the motion of the bodies (and provide the gravitational waveform) beyond the adiabatic approximation; and (iii) in 
principle they can be extended to higher PN orders. More importantly, they can provide initial dynamical data for 
the two BHs at the beginning of the plunge (such as their positions and momenta) , which can be used (in principle) 
in numerical relativity to help build the initial gravitational data (the metric and its time derivative) and then to 
evolve the full Einstein equations through the merger phase. However, these resummation methods are based on 
some assumptions that, although plausible, have not been proved: for example, when the orbital energy and the 
gravitational flux functions are derived in the comparable-mass case, it is assumed that they are smooth deformations 
of the analogous quantities in the test-mass limit. Moreover, in the absence of both exact solutions and experimental 
data, we can test the robustness and reliability of the resummation methods only by internal convergence tests. 

In this paper we follow a more conservative point of view. We shall maintain skepticism about waveforms emitted 
by BBH with M — 10-40M Q and evaluated from PN calculations, as well as all other waveforms ever computed for 
the late BBH inspiral and plunge, and we shall develop families of search templates that incorporate this skepticism. 
More specifically, we shall be concerned only with detecting BBH GWs, and not with extracting physical parameters, 
such as masses and spins, from the measured GWs. The rationale for this choice is twofold. First, detection is the 
more urgent problem at a time when GW interferometers are about to start their science runs; second, a viable 
detection strategy must be constrained by the computing power available to process a very long stream of data, while 
the study of detected signals to evaluate physical parameters can concentrate many resources on a small stretch of 
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detector output. In addition, as we shall see in Sec. IVII and briefly discuss in Sec. IVIDI the different PN methods 
will give different parameter estimations for the same waveform, making a full parameter extraction fundamentally 
difficult. 

This is the strategy that we propose: we guess (and hope) that the conjunction of the waveforms from all the 
post-Newtonian models computed to date spans a region in signal space that includes (or almost includes) the true 
signal. We then choose a detection (or effective) template family that approximates very well all the PN expanded and 
resummed models (henceforth denoted as target models). If our guess is correct, the effectualness 26] of the effective 
model in approximating the targets (i.e., its capability of reproducing their signal shapes) should be indicative of its 
effectualness in approximating the true signals. Because our goal is the detection of BBH GWs, we shall not require 
the detection template family to be faithful 26] (i.e., to have a small bias in the estimation of the masses). 

As a backup strategy, we require the detection template family to embed the targets in a signal space of higher 
dimension (i.e., with more parameters), trying to guess the functional directions in which the true signals might lie 
with respect to the targets (of course, this guess is rather delicate!). So, the detection template families constructed 
in this paper cannot be guaranteed to capture the true signal, but they should be considered as indications. 

This paper is organized as follows. In Sec.|nJ we briefly review the theory of matched- filtering GW detections, which 
underlies the searches for GWs from inspiraling binaries. Then in Sees. Hill IIVI and Ivl we present the target models 
and give a detailed analysis of the differences between them, both from the point of view of the orbital dynamics and 
of the gravitational waveforms. More specifically, in Sec. IIIII we introduce the two-body adiabatic models, both PN 
expanded and resummed; in Sec. IIVI we introduce nonadiabatic approximations to the two-body dynamics; and in 
Sec. we discuss the signal-to- noise ratios obtained for the various two-body models. Our proposals for the detection 
template families are discussed in the Fourier domain in Sec. IVII and in the time domain in Sec. I VIII where we also 
build the mismatch metric [2?l |28| for the template banks and use it to evaluate the number of templates needed for 
detection. Section IVIIII summarizes our conclusions. 

Throughout this paper we adopt the LIGO noise curve given in Fig. Q] and Eq. l(2"5|) . and used also in Ref. [l2| . 
Because the noise curve anticipated for VIRGO [see Fig. ^ is quite different (both at low frequencies, and in the 
location of its peak-sensitivity frequency) our results cannot be applied naively to VIRGO. We plan to repeat our 
study for VIRGO in the near future. 

II. THE THEORY OF MATCHED-FILTERING SIGNAL DETECTION 

The technique of matched-filtering detection for GW signals is based on the systematic comparison of the measured 
detector output s with a bank of theoretical signal templates {tti} that represent a good approximation to the class 
of physical signals that we seek to measure. This theory was developed by many authors over the years, who have 
published excellent expositions H3, E^, [H, H^, ISO, H3, E3, H3, H3, |1M H3, E3, E3, H^l hi the following, we summarize 
the main results and equations that are relevant to our purposes, and we establish our notation. 

A. The statistical theory of signal detection 

The detector output s consists of noise n and possibly of a true gravitational signal hi (part of a family {hi} of 
signals generated by different sources for different source parameters, detector orientations, and so on). Although 
we may be able to characterize the properties of the noise in several ways, each separate realization of the noise is 
unpredictable, and it might in principle fool us by hiding a physical signal (hence the risk of a false dismissal) or 
by simulating one (false alarm). Thus, the problem of signal detection is essentially probabilistic. In principle, we 
could try to evaluate the conditional probability P(h\s) that the measured signal s actually contains one of the hi. 
In practice, this is inconvenient, because the evaluation of P(h\s) requires the knowledge of the a priori probability 
that a signal belonging to the family {hi} is present in s. 

What we can do, instead, is to work with a statistic (a functional of s and of the hi) that (for different realizations 
of the noise) will be distributed around low values if the physical si gna l hi is absent, and around high value if the 
signal is present. Thus, we shall establish a decision rule as follows |32j: we will claim a detection if the value of a 
statistic (for a given instance of s and for a specific hi) is higher than a predefined threshold. We can then study 
the probability distribution of the statistic to estimate the probability of false alarm and of false dismissal. The steps 
involved in this statistical study are easily laid down for a generic model of noise, but it is only in the much simplified 
case of normal noise that it is possible to obtain manageable formulas; and while noise will definitely not be normal 
in a real detector, the Gaussian formulas can still provide useful guidelines for the detection problems. Eventually, 
the statistical analysis of detector search runs will be carried out with numerical Montecarlo techniques that make 
use of the measured characteristics of the noise. So throughout this paper we shall always assume Gaussian noise. 
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The statistic that is generally used is based on the symmetric inner product (g, h) between two real signals g and 
h, which represents essentially the cross-correlation between g and h, weighted to emphasize the correlation at the 
frequencies where the detector sensitivity is better. We follow Cutler and Flanagan's conventions |H3 and define 



, , -°° ~9*(f)h(f) [ + °° ~9*(f)h(f) 



— oo 



where 5 n (/), the one-sided noise power spectral density, is given by 



n*(/i)n(/ 2 ) = - / 2 )5„(/i) for AX), (2) 

and S n (fi) = for fx < 0. We then define the signal-to-noise ratio p (for the measured signal s after filtering by hi), 
as 



where the equality follows because (hi,n)(n,hi) = (hi, hi) (see, e.g., |32j). In the case of Gaussian noise, it can be 
proved that this filtering technique is optimal, in the sense that it maximizes the probability of correct detection for 
a given probability of false detection. 

In the case when s = n, and when noise is Gaussian, it is easy to prove that p is a normal variable with a mean of 
zero and a variance of one. If instead s — hi + n, then p is a normal variable with mean y/ (hi, hi) and unit variance. 
The threshold p* for detection is set as a tradeoff between the resulting false-alarm probability, 

/ i r+oo i 

F=\— e-P-' 2 dp = -erfc (p,/y/2) (4) 

V 2vr J pt 2 

(where erfc is the complementary error function |36|). and the probability of correct detection 



V = -erfc [(p, - y/(hi,hi))/V2] (5) 



(the probability of false dismissal is just 1 — 2?). 



B. Template families and extrinsic parameters 

We can now go back to the initial strategy of comparing the measured signal against a bank of Mi templates {u{\ 
that represent a plurality of sources of different typ es and p hysical parameters. For each stretch s of detector output, 
we shall compute the signal-to-noise ratio (s, Ui)/ \J (uj, Ui) for all the Ui, and then apply our rule to decide whether 
the physical signal corresponding to any one of the Ui is actually present within s 4|. Of course, the threshold p* 
needs to be adjusted so that the probability J-'tot of false alarm over all the templates is still acceptable. Under 
the assumption that all the inner products (n, u,-) of the templates with noise alone are statistically independent 
variables [this hypothesis entails (ui,Uj) ~ 0], J- tot is just 1 — (1 — J 7 )^ ~ Ni? '. If the templates are not statistically 
independent, this number is an upper limit on the false alarm rate. However, we first need to note that, for any 
template Ui, there are a few obvious ways (parametrized by the so-called extrinsic parameters) of changing the signal 
shape that do not warrant the inclusion of the modified signals as separate templates [7(j 

The extrinsic parameters are the signal amplitude, phase and time of arrival. Any true signal h can be written in 
all generality as 

h(t) = A h a h [t - t h ] cos[$ h (t - t h ) + <f> h ], (6) 

where ah(t) = for t < 0, where ^(O) = 0, and where ah(t) is normalized so that (h, h) = A\. While the template 
bank {u.i} must contain signal shapes that represent all the physically possible functional forms a(t) and $(£), it is 
possible to modify our search strategy so that the variability in Ah, 4>h and th is automatically taken into account 
without creating additional templates. 

The signal amplitude is the simplest extrinsic parameter. It is expedient to normalize the templates it, so that 
(ui,Ui) = 1, and p(wj) = (s,Ui). Indeed, throughout the rest of this paper we shall always assume normalized 
templates. If s contains a scaled version hi = Aui of a template Ui (here A is known as the signal strength), then 
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p(ui) = A. However, the statistical distribution of p is the same in the absence of the signal. Then the problem of 
detection signals of known shape and unknown amplitude is easily solved by using a single normalized template and 
the same threshold as used for the detection of completely known signals 32]. Quite simply, the stronger an actual 
signal, the easier it will be to reach the threshold. 

We now look at phase, and we try to match h with a continuous one-parameter subfamily of templates u(4> t ] t) — 
ah(t) cos[$/,.(£) + (f>t\. It turns out that for each time signal shape {a(t), we need to keep in our template bank 

only two copies of the corresponding Ui, for <pt — and <pt — tt/2, and that the signal to noise of the detector output 
s against u^, for the best possible value of (f>t, is automatically found as |32J 



P4> 



max (s, u 



(&)> = ^/\( S ,u t (o))\ 2 + \{s,M*/m 2 , 



(7) 



where iti(O) and Ui(ir/2) have been orthonormalized. The statistical distribution of the phase-maximized statistic 
for the case of (normal) noise alone, is the Raleigh distribution |32| 



Po(P4>) = P4> e P * /2 > 
and the false-alarm probability for a threshold p<p* is just 



(8) 



(9) 



Throughout this paper, we will find it useful to consider inner products that are maximized (or minimized) with 
respect to the phases of both templates and reference signals. In particular, we shall follow Damour, Iyer and 
Sathyaprakash in making a distinction between the best match or maxmax match 



maxmax(/i, Ui) = maxmax (/i(0^), Ui(<pt)), 



(10) 



which represents the most favorable combination of phases between the signals h and itj, and the minmax match 

minmax ( ft,, iti) = minmax (h(<ph), Wj(<^t)), (11) 



which represents the safest estimate in the realistic situation, where we cannot choose the phase of the physical 
measured signal, but only of the template used to match the signal. Damour, Iyer and Sathyaprakash [see Appendix 
B of Ref. [2q] show that both quantities are easily computed as 



where 



maxmax 
minmax 



.4 
B 
C 



A + B 



± 



A-B 



C 1 



1/2' 



1/2 



(/i(0), Ul (0)) 2 + (M0), Ul (V2)) 2 , 
(h(n/2), Ul (0)) 2 + (&(tt/2), u t (n/2)y 

(ft(0),tii(0))(M7r/2),ii<(0)) + 
(/i(0) )Ui (7r/2))(ft(7r/2), Ui (7r/2)). 



(12) 



(13) 
(14) 
(15) 



In these formulas we have assumed that the two bases {h(0), /i(7r/2)} and {ui(0), Ui(Tr/2)} have been orthonormalized. 

The time of arrival th is an extrinsic parameter because the signal to noise for the normalized, time-shifted template 
u(t — to) against the signal s is just 



(s,u(t )) = 4 Re 



+ 00 ~* 



CfKf) 

Sn(f) ' 



df, 



(16) 



where we have used a well-known property of the Fourier transform of time-shifted signals. These integrals can be 
computed at the same time for all the time of arrivals {to}, using a fast Fourier transform technique that requires 
~ N s log N s operations (where N s is the number of the samples that describe the signals) as opposed to ~ N% required 
to compute all the integrals separately |37j . Then we can look for the optimal to that yields the maximum signal to 
noise. 
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We now go back to adjusting the threshold for a search over a vast template bank, using the estimate 10 for 
the false-alarm probability. Assuming that the statistics p^ for each signal shape and starting time are independent, 
we require that 

e-"^ 2 - — , (17) 

1 v timcs i v shapes 

or 



/J* ~ J 2(l0g iVtimes + log ^shapes ~ log Ftot) ■ (18) 

It is generally assumed that iVtimes ~3x 10 10 (equivalent to templates displaced by 0.01 s over one year \vL, 38]) and 
that the false-alarm probability Ttot ~ 10~ 3 . Using these values, we find that an increase of /?* by about ~ 3% is 
needed each time we increase iV s h ape s by one order of magnitude. So there is a tradeoff between the improvement in 
signal-to-noisc ratio obtained by using more signal shapes and the corresponding increase in the detection threshold 
for a fixed false-alarm probability. 



C. Imperfect detection and discrete families of templates 



There are two distinct reasons why the detection of a physical signal h by matched filtering with a template bank 
{ui} might result in signal-to-noise ratios lower than the optimal signal-to-noise ratio, 



Popt 



First, the templates, understood as a continuous family {u(X )} of functional shapes indexed by one or more intrinsic 
parameters X A (such as the masses, spins, etc.), might give an unfaithful representation of h, introducing errors in the 
representation of the phasing or the amplitude. The loss of signal to noise due to unfaithful templates is quantified 
by the fitting factor FF, introduced by Apostolatos [3|j, and defined by 

FF(h,u(X A )) = maXXA ^ XA)) . (20) 
V(h,h) 

In general, we will be interested in the FF of the continuous template bank in representing a family of physical 
signals {h(8 A )}, dependent upon one or more physical parameters 9 A : so we shall write FF(9 A ) — FF(h(9 A ), u(X A )). 
Although it is convenient to index the template family by the same physical parameters 9 A that characterize h(9 A ), 
this is by no means necessary; the template parameters X A might be a different number than the physical parameters 
(indeed, this is desirable when the 9 A get to be very many), and they might not carry any direct physical meaning. 
Notice also that the value of the FF will depend on the parameter range chosen to maximize the X A . 

The second reason why the signal-to- noise will be degraded with respect to its optimal value is that, even if our 
templates are perfect representations of the physical signals, in practice we will not adopt a continuous family of 
templates, but we will be limited to using a discrete bank hu = u(X A )}. This loss of signal to noise depends 
on how finely templates are laid down over parameter space ^3 0, |43 > a n °ti° n of metric in template space (the 
mismatch metric [271 1281 143| ) can be used to guide the disposition of templates so that the loss (in the perfect-template 
abstraction) is limited to a fixed, predetermined value, the minimum match MM, introduced in Refs. [2H Ecj. and 
defined by 

MM = min max {u(X A ), u(X A )) = minmax (u(X A ), u(X A + AX A )), (21) 

where AX A = X A — A . The mismatch metric <7bc(^ ) f°r the template space {u(X A )} is obtained by expanding the 
inner product (or match) (u(A ),u(A + AX A )) about its maximum of 1 at AX A = 0: 



1 rP-M 

(u(X A ), u(X A + AX A )) = M(X A , X A + AX A ) = 1 + 



2 dAX B dAX c 



AX B AX C + • • • , (22) 



so the mismatch 1 — M between u(X ) and the nearby template u(A + AX A ) can be seen as the square of the proper 
distance in a differential manifold indexed by the coordinates X A pflj . 

1 — M(X A , X A + AX A ) — gscAX B AX C : (23) 
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where 



1 d 2 M 
9BC ~ ~2 dAX B dAX c 



(24) 



If, for simplicity, we lay down the n-dimensional discrete template bank {u(X A )} along a hypercubical grid of cellsize 
dl in the metric qab (a grid in which all the templates on nearby corners have a mismatch of dl with each other), 
the minimum match occurs when X A lies exactly at the center of one of the hypercubes: then 1 — MM = n(dl/2) 2 . 
Conversely, given MM, the volume of the corresponding hypercubes is given by Vmm = (2y(l — MM)/n) n . The 
number of templates required to achieve a certain MM is obtained by integrating the proper volume of parameter 
space within the region of physical interest, and then dividing by V^mm: 

Af[g, MM] = f _ (25) 

2 v / [l-MM]/( 

In practice, if the metric is not constant over parameter space it will not be possible to lay down the templates on an 
exact hypercubical grid of cellsize dl, so A/" will be somewhat higher than predicted by Eq. I|25(l . However, we estimate 
that this number should be correct within a factor of two, which is adequate for our purposes. 

In the worst possible case, the combined effect of unfaithful modeling (FF < 1) and discrete template family 
(MM < 1) will degrade the optimal signal to noise by a factor of about FF + MM — 1. This estimate for the total 
signal-to-noise loss is exact when, in the space of signals, the two segments that join h(9 A ) to its projection u(X A ) 
and u(X A ) to the nearest discrete template u(Xf) can be considered orthogonal: 

(h(9 A )-u(X A ),uCX A )-u(X A ))cO. (26) 

This assumption is generally very accurate if FF and MM are small enough, as in this p aper ; so we will adopt this 
estimate. However, it is possible to be more precise, by defining an external metric g\ B [27l |44| that characterizes 
directly the mismatch between h(9 A ) and a template u(X A + AX A ) that is displaced with respect to the template 
u(X A ) that is yields the maximum match with h(9 A ). 

Since the strength of gravity- wave signals scales as the inverse of the distance jj^l ; the matched-filtering scheme, 
with a chosen sig nal-to- noise threshold p», will allow the reliable detection of a signal h, characterized by the signal 
strength Ad = y (h, h) at the distance do, out to a maximum distance 

^max _ -Ada (27) 

d p* 

If we assume that the measured GW events happen with a homogeneous event rate throughout the accessible portion 
of the universe, then the detection rate will scale as e?^ ax . It follows that the use of unfaithful, discrete templates 
{ui} to detect the signal h will effectively reduce the signal strength, and therefore d max , by a factor FF + MM — 1. 
This loss in the signal-to-noise ratio can also be seen as an increase in the detection threshold p* necessary to achieve 
the required false-alarm rate, because the imperfect templates introduce an element of uncertainty. In either case, 
the detection rate will be reduced by a factor (FF + MM — l) 3 . 



D. Approximations for detector noise spectrum and gravitational- wave signal 



For LIGO-I we use the analytic fit to the noise power spectral density given in Ref. [TJ, and plotted in Fig.^ 
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4.49 



0.16 



0.52 + 0.32 



(28) 



where fa — 150 Hz. The first term in the square brackets represents seismic noise, the second and third, thermal 
noise, and the fourth, photon shot noise. 

Throughout this paper, we shall compute BBH waveforms in the quadrupole approximation (we shall compute 
the phase evolution of the GWs with the highest possible accuracy, but we shall omit all harmonics higher than the 
quadrupole, and we shall omit post-Newtonian corrections to the amplitude; this is a standard approach in the field, 
see, e.g., Hfl). The signal received at the interferometer can then be written as pi l3l| 



e 



h(t) = —M V (TrMf GW ) 2 / 3 



COS V?GW, 



(29) 



8 




/(Hz) 



FIG. 1: Square root of the noise spectral density \/S n (f) versus frequency /, for LIGO-I [Eq. 11280 ]. and VIRGO (from Tab. 
IV of Ref. 111). 



model 


shorthand 


evolution equation 


section 


adiabatic model with 
Taylor-expanded energy £(v) and 
flux F{v) 


T(nPN, raPN; 8) 


energy-balance equation 


Sec.|TTTA| 


adiabatic model with Pade-expanded 
energy £ (v) and flux J-(v) 


P(nPN,mPN;6») 


energy-balance equation 


Sec.|TTTBl 


adiabatic model with 
Taylor-expanded energy £ (v) and 
flux J-(v) in the stationary-phase 
approximation 


SPA(nPN = mPN) 


energy-balance equation in the freq. 
domain 


Sec.lVIFI 


nonadiabatic Hamiltonian model 
with Taylor-expanded GW flux 


HT(nPN, mPN; 8) 


Hamilton equations 


Sec.|lVA| 


nonadiabatic Hamiltonian model 
with Pade-expanded GW flux 


HP(nPN,mPN;6>) 


Hamilton equations 


Sec.|iV^l 


nonadiabatic Lagrangian model 


L(nPN, mPN) 


F = ma 


Sec.HVBI 


nonadiabatic effective-one-body 
model with Taylor-expanded GW flux 


ET(nPN, mPN; 8; zi,z 2 ) 


eff. Hamilton equations 


Sec.|TVTl 


nonadiabatic effective-one-body 
model with Pade-expanded GW flux 


EP(nPN,raPN;0;zi,2 2 ) 


eff. Hamilton equations 


Sec.|TvT!l 



TABLE I: Post-Newtonian models of two-body dynamics defined in this paper. The notation X(nPN, mPN; 9) denotes the 
model X, with terms up to order nPN for the conservative dynamics, and with terms up to order mPN for radiation-reaction 
effects; for m > 3 we also need to specify the arbitrary flux parameter 8 (see Sec. 1111 A"l ; for n > 3, the effective-one-body 
models need also two additional parameters z\ and Zi (see Sec. II V CI . 

where / and (/?GW are the instantaneous GW frequency and phase at the time t, g?l is the luminosity distance, M 
and Tj are respectively the BBH total mass mi + m-i and the dimcnsionless mass ratio m 1 m2/M 2 , and where we have 
taken G = c = 1. The coefficient G depends on the inclination of the BBH orbit with respect to the plane of the sky, 
and on the polarization and direction of propagation of the GWs with respect to the orientation of the interferometer. 
Finn and Chernoff [3l| examine the distribution of O, and show that max = 4, while rmsG = 8/5. We shall use this 
last value when we compute optimal signal-to-noise ratios. The waveform given by Eq. (|29J) . after dropping the factor 
QMrj/dh, is known as restricted waveform. 

III. ADIABATIC MODELS 

We turn, now, to a discussion of the currently available mathematical models for the inspiral of BBHs. Table [I] 
shows a list of the models that we shall consider in this paper, together with the shorthands that we shall use to 



9 



denote them. We begin in this section with adiabatic models. BBH adiabatic models treat the orbital inspiral as a 
quasistationary sequence of circular orbits, indexed by the invariantly defined velocity 

v = {Mp) 1 ' 3 = (^Af/ GW ) 1/3 . (30) 

The evolution of the inspiral (and in particular of the orbital phase ip) is completely determined by the energy-balance 
equation 

= (31) 

This equation relates the time derivative of the energy function £(v) (which is given in terms of the total relativistic 
energy £ to t by £ = £tot — m\ — m^, and which is conserved in absence of radiation reaction) to the gravitational flux 
(or luminosity) function T(v). Both functions are known for quasicircular orbits as a PN expansion in v. It is easily 
shown that Eq. i|3*T|) is equivalent to the system (see, e.g., Ref. 26] ) 

alt M ' dt Md£(v)/dv' 1 ' 

In accord with the discussion around Eq. 1)29(1 . we shall only consider the restricted waveform h(t) = v 2 cos <pgw(^)> 
where the GW phase (pew is twice the orbital phase ip. 



A. Adiabatic PN expanded models 

The equations of motion for two compact bodies at 2.5PN order were first derived in Refs. The 3PN equations 
of motion have been obtained by two separate groups of researchers: Damour, Jaranowski and Schafer |46| used the 
Arnowitt-Deser-Misner (ADM) canonical approach, while Blanchet, Faye and de Andrade 1471 worked with the PN 
iteration of the Einstein equations in the harmonic gauge. Recently Damour and colleagues |48|. working in the ADM 
formalism and applying dimensional regularization, determined uniquely the static parameter that enters the 3PN 
equations of motion [461 l47j and that was until then unknown. In this paper we shall adopt their value for the static 
parameter. Thus at present the energy function £ is known up to 3PN order. 

The gravitational flux emitted by compact binaries was first computed at 1PN order in Ref. 49]. It was subse- 
quently determined at 2PN order with a formalism based on multipolar and post-Minkowskian approximations, and, 
independently, with the direct integration of the relaxed Einstein equations |50|. Nonlinear effects of tails at 2.5PN 
and 3.5PN orders were computed in Refs. More recently, Blanchet and colleagues derived the gravitational-flux 
function for quasicircular orbits up to 3.5PN order [52l l53|. However, at 3PN order |53j the gravitational-flux 
function depends on an arbitrary parameter 9 that could not be fixed in the regularization scheme used by these 
authors. 



PN energy and flux 

Denoting by £t n and J-t n the iV th -order Taylor approximants (T-approximants) to the energy and the flux func- 
tions, we have 

N 

£t 2N (v) = £Ncwt(v) £k(v) v 2k , (33) 

N 

T Tn (v) = J- Ncwt (v) Fk (r?) v k , (34) 

where "Newt" stands for Newtonian order, and the subscripts 2N and N stand for post 2Ar -Newtonian and post^- 
Newtonian order. The quantities in these equations are 

1 32 

^Newt(w) = ~2 Tll>2 ' :F Newt(w) = — if V W , (35) 

£o( V ) = l, £i(rj) = -|~ Y2 > £2 ^ = -Y + Y V ~i' (36) 
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FIG. 2: Normalized flux function Ft n /J~Newt versus v, at different PN orders for equal-mass binaries, rj = 0.25. Note that 
the 1.5PN and 2PN flux, and the 3PN and 3.5PN flux, are so close that they cannot be distinguished in these plots. The two 
long-dashed vertical lines correspond to v ~ 0.18 and v ~ 0.53; they show the velocity range that corresponds to the LIGO 
frequency band 40 < faw < 240 Hz for BBHs with total mass in the range 1O-4OM0. 
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Here rj = m-\ mil(m\ 



- n 7.J my i yi 1 1\ -i- TO2) 2 , 1e is Euler's gamma, and is the arbitrary 3PN flux parameter |52l l53| . From Tab. 
I of Ref. [531 we read that the extra number of GW cycles accumulated by the PN terms of a given order decreases 
(roughly) by an order of magnitude when we increase the PN order by one. Hence, we find it reasonable to expect 
that at 3PN order the parameter 9 should be of order unity, and we choose as typical values 9 — 0,±2. (Note for v3 
of this paper on g r-qc: Eqs. (39) and (41) are now revised as per Ref. |6jj; the parameter 9 has been determined to 
be 1039/4620 

In Fig. [2] we plot the normalized flux J~t n / ^Ncwt as a function of v at various PN orders for the equal mass case 
rj = 0.25. To convert v to a GW frequency we can use 



/GW - 3.2 x 10 4 {^\ , 



(42) 



The two long-dashed vertical lines in Fig. |21 correspond to v ~ 0.18 and v ~ 0.53; they show the velocity range that 
corresponds to the LIGO frequency band 40 < /gw < 240 Hz for BBHs with total mass in the range 10-40M Q . 
At the LIGO-I peak-sensitivity frequency, which is 153 Hz according to our noise curve, and for a (10+10)Af Q 
BBH, we have v ~ 0.362; and the percentage difference between subsequent PN orders is Newt — > 1PN : —58%; 
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FIG. 3: In the left panel, we plot the energy function £t n versus v, at different PN orders, for r\ — 0.25. The two long-dashed 
vertical lines in the left figure correspond to v ~ 0.18 and v ~ 0.53; they show the velocity range that corresponds to the 
LIGO frequency band 40 < /gw < 240 Hz, for BBHs with total mass in the range 1O-4OM . In the right panel, we plot the 
percentage difference 8£t n = 100 \(£t n+1 — £t n )/£t n \ versus the total mass M, for N = 1,2, at the LIGOT peak-sensitivity 
GW frequency, / pcak = 153 Hz [note: i; pcak = {nM f pcak ) 1/3 ]. 



M 




/gw 


(Hz) at MECO 






/gw(Hz) 


at ISCO 






T (1PN) 


T (2PN) 


T (3PN) 


P (2PN) 


P (3PN) 


H (1PN) 


E (1PN) 


E (2PN) 


E (3PN) 


(5 + 5)M 


3376 


886 


832 


572 


866 


183 


446 


473 


570 


(10 + 1O)M 


1688 


442 


416 


286 


433 


92 


223 


236 


285 


(15 + 15)M 


1125 


295 


277 


191 


289 


61 


149 


158 


190 


(20 + 2O)M 


844 


221 


208 


143 


216 


46 


112 


118 


143 



TABLE II: Location of the MECO/ISCO. The first six columns show the GW frequency at the Maximum binding Energy for 
Circular Orbits (MECO), computed using the T- and P-approximants to the energy function; the remaining columns show the 
GW frequency at the Innermost Stable Circular Orbit (ISCO), computed using the H-approximant to the energy, and using 
the EOB improved Hamiltonian l|9()|l with z\ = 22 = 0. For the H-approximant the ISCO exists only at 1PN order. 



1PN -> 1.5PN : +142%; 1.5PN -> 2PN : -0.2%; 2PN -> 2.5PN : -34%; 2.5PN -> 3PN(# = 0) : +43%; 3PN -> 
3.5PN(# = 0) : +0.04%. The percentage difference between the 3PN fluxes with 6 = ±2 is ~ 7%. It is interesting 
to notice that while there is a big difference between the 1PN and 1.5PN orders, and between the 2PN and 2.5PN 
orders, the 3PN and 3.5PN fluxes are rather close. Of course this observation is insufficient to conclude that the PN 
sequence is converging at 3.5PN order. 

In the left panel of Fig. [31 we plot the T-approximants for the energy function versus v, at different PN orders, 
while in the right panel we plot (as a function of the total mass M, and at the LIGO-I peak-sensitivity GW frequency 
/peak = 153 Hz) the percentage difference of the energy function between T-approximants to the energy function of 
successive PN orders. We note that the 1PN and 2PN energies are distant, but the 2PN and 3PN energies are quite 
close. 



Definition of the models 

The evolution equations 132fl for the adiabatic inspiral lose validity (the inspiral ceases to be adiabatic) a little 
before v reaches UmecO' wnere MECO stands for Maximum-binding-Energy Circular Orbit [54I l6o| . This ujy[ECO ^ s 
computed as the value of v at which d£r N {v) / 'dv = 0. In building our adiabatic models we evolve Eqs. (|32l) right up 
to Vmego and stop there. We shall refer to the frequency computed by setting v = wmeco i n Eq. (|42|) as the ending 
frequency for these waveforms, and in Tab. [H] we show this frequency for some BH masses. However, for certain 
binaries, the 1PN and 2.5PN flux functions can go to zero before v = u^ECO t see Fig- 12- In those cases we choose as 
the ending frequency the value of / = v 3 /(nM) where T{v) becomes 10% of J-Newt(^)- [When using the 2.5PN flux, 
our choice of the ending frequency differs from the one used in Ref. where the authors stopped the evolution at 
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TV 


(5 + 20) Mq 


(Tat, Tjv+i) 
(10+ 10) M© 


(15 + 15)M© 





0.432 


0.553 (0.861, 19.1, 0.241) 


0.617 


1 


0.528 [0.638] 


0.550 (0.884, 22.0, 0.237) 


0.645 [0.712] 


2 = +2) 


0.482 [0.952] 


0.547 (0.841, 18.5, 0.25) 


0.563 [0.917] 


20= -2) 


0.457 [0.975] 


0.509 (0.821, 18.7, 0.241) 


0.524 [0.986] 



TABLE III: Test for the Cauchy convergence of the T-approximants. The values quoted are maxmax matches obtained by 
maximizing with respect to the extrinsic parameters, but not to the intrinsic parameters (i.e., the matches are computed for 
T waveforms with the same masses, but different PN orders). Here we define T = T(0,0), Ti = T(l, 1.5), T 2 = T(2,2.5), 
T3 = T(3, 3.5,6). In the Newtonian case, To = (0,0), the MECO does not exist and we stop the integration of the balance 
equation at v = 1. The values in brackets, "[...]," are obtained by setting T2 = T(2, 2) instead of T(2, 2.5); the values in 
parentheses, "(...)," are obtained by maximizing with respect to the extrinsic and intrinsic parameters, and they are shown 
together with the Tjv + i parameters M and n where the maxima are achieved. In all cases the integration of the equations is 
started at a GW frequency of 20 Hz. 




/(Hz) 



FIG. 4: Frequency-domain amplitude versus frequency for the T-approximated (restricted) waveforms, at different PN orders, 
for a (15 + 15)M© BBH. The T(3, 3.5, 6 = 0) curve, not plotted, is almost identical to the T(3, 3,9 = 0) curve. 



the GW frequency corresponding to the Schwarzschild innermost stable circular orbit. For this reason there are some 
differences between our overlaps and theirs.] 

We shall refer to the models discussed in this section as T(nPN, mPN), where nPN (mPN) denotes the maximum 
PN order of the terms included for the energy (the flux). We shall consider (nPN, mPN) = (1, 1.5), (2, 2), (2, 2.5) and 
(3, 3.5, 9) [at 3PN order we need to indicate also a choice of the arbitrary flux parameter 9]. 

Waveforms and matches 

In Tab. IIIII for three typical choices of BBH masses, we perform a convergence test using Cauchy 's criterion (2||, 
namely, the sequence Tjv converges if and only if for each k, (Tjv,Tjv+fc) — > 1 as TV — > 00. One requirement of this 
criterion is that (TV, Tjv+i) 1 as TV -> 00, and this is what we test in Tab. EH setting T N = T(TV,TV + 0.5). The 
values quoted assume maximization on the extrinsic parameters but not on the intrinsic parameters. [For the case 
(10 + 10)TW©, we show in parentheses the maxmax matches obtained by maximizing with respect to the intrinsic and 
extrinsic parameters, together with the intrinsic parameters M and r\ of T/v+i where the maxima are attained.] These 
results suggest that the PN expansion is far from converging. However, the very low matches between TV = 1 and 
TV = 2, and between TV = 2 and TV = 3, are due to the fact that the 2.5PN flux goes to zero before the MECO can be 
reached. If we redefine T 2 as T(2, 2) instead of T(2, 2.5), we obtain the higher values shown in brackets is Tab. IIIII 

In Fig. 0] we plot the frequency-domain amplitude of the T-approximated waveforms, at different PN orders, for a 
(15 + 15) Mq BBH. The Newtonian amplitude, A^ e wt(f) = f~ 7 ^ 6 , is also shown for comparison. In the T(l, 1) and 
T(2, 2.5) cases, the flux function goes to zero before v = ^ M ecO' this means that the radiation-reaction effects become 



13 



negligible during the last phase of evolution, so the binary is able to spend many cycles at those final frequencies, 
skewing the amplitude with respect to the Newtonian result. For T(2,2), T(3,3) and T(3, 3.5), the evolution is 
stopped at v = ^mecO' anc ^ although /meco — 270-300 Hz (see Tab. Hill the amplitude starts to deviate from J -7 / 6 
around 100 Hz. This is a consequence of the abrupt termination of the signal in the time domain. 

The effect of the arbitrary parameter 9 on the T waveforms can be seen in Tab. IXIIII in the intersection between 
the rows and columns labeled T(3, 3.5, +2) and T(3, 3.5,— 2). For three choices of BBH masses, this table shows 
the maxmax matches between the search models at the top of the columns and the target models at the left end of 
the rows, maximized over the mass parameters of the search models in the columns. These matches are rather high, 
suggesting that for the range of BBH masses we are concerned, the effect of changing 9 is just a remapping of the 
BBH mass parameters. Therefore, in the following we shall consider only the case of 9 = 0. 

A quantitative measure of the difference between the T(2, 2), T(2, 2.5) and T(3, 3.5) waveforms can be seen in Tab. 
IXII in the intersection between the rows and columns labeled T(. . .). For four choices of BBH masses, this table shows 
the maxmax matches between the search models in the columns and the target models in the rows, maximized over 
the search-model parameters M and rj; in the search, r\ is restricted to its physical range < rj < 1/4, where 
corresponds to the test-mass limit, while 1/4 is obtained in the equal-mass case. These matches can be interpreted as 
the fitting factors [see Eq. 1|20|) ] for the projection of the target models onto the search models. For the case T(2, 2.5) 
the values are quite low: if the T(3, 3.5) waveforms turned out to give the true physical signals and if we used the 
T(2, 2.5) waveforms to detect them, we would lose ~ 32-49% of the events. The model T(2, 2) would do match better, 
although it would still not be very faithful. Once more, the difference between T(2, 2) and T(2, 2.5) is due to the fact 
that the 2.5PN flux goes to zero before the BHs reach the MECO. 



B. Adiabatic PN resummed methods: Pade approximants 



The PN approximation outlined above can be used quite generally to compute the shape of the GWs emitted by 
BNSs or BBHs, but it cannot be trusted in the case of binaries with comparable masses in the range M ~ 10-40M Q , 
because for these sources LIGO and VIRGO will detect the GWs emitted when the motion is strongly relativistic, 
and the convergence of the PN series is very slow. To cope with this problem, Damour, Iyer and Sathyaprakash |2^| 
proposed a new class of models based on the systematic application of Pade resummation to the PN expansions of 
£(v) and f(v). This is a standard mathematical technique used to accelerate the convergence of poorly converging 
or even divergent power series. 

If we know the function g(v) only through its Taylor approximant Gn(v) = go + giv + ■ ■ ■ + v N = Ti^[g(v)], 
the central idea of Pade resummation [55| is the replacement of the power series Gn(v) by the sequence of rational 
functions 

with M + K = N and Tm+k[Pk ( v )] = Gn( v ) (without loss of generality, we can set by = 1). We expect that for 
M, K — > +oo, Pjf [g(v)] will converge to g(v) more rapidly than T^[g{v)] converges to g(v) for N — > +oo. 



PN energy and flux 



Damour, Iyer and Sathyaprakash |2t| . and then Damour, Schafer and Jaranowski |T(| . proposed the following 
Pade-approximated (P-approximated) £p N (v) and Tp N {v) (for AT = 2,3): 

£p n = V /l + 27/ v /lT^>)-l-l, (44) 
32 1 

Fp n = ^V 2 v W - -p^/i^M, (45) 

5 1 - «A£L 



where 



e»(v) = - y + b-C-ln+W)/ t (46) 



i+i^-(3-i^) 



2 1 - (1 + ±7? + waft)) v 2 - (3 - ffry - (l + ±ry) w 3 ( V )) v 4 



e Pa (v) = * ' J ' ' * +y * • , (47) 

1 — w^yrj) v z 
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v 

FIG. 5: Normalized flux function J-p N / Tkewt versus v, at different PN orders. The two long-dashed vertical lines give v ~ 0.18 
and v ~ 0.53; they show the velocity range that corresponds to the LIGO frequency band 40 < /gw < 240 Hz for BBHs with 
total mass in the range 10-40Mq . Compare with Fig. 
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Here the dimensionless coefficients Ci depend only on r]. The c^'s are explicit functions of the coefficients fk (k — 1, ...5), 



Cl 



f'4 



C5 



-A, 



C2 = /l 



A _ /l /3 ~ f l 

h ' 



C3 = 



/l (/l ' / 2 ) ' 



/l (/I + fj + /l /4 ~ h (2 /l / 3 + h)) 
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where 



fk — J~k — 



fc-1 



pole 



(51) 
(52) 
(53) 

(54) 



Here J-k is given by Eqs. (|38(I - H41|I [for k = 6 and k — 7, the term —856/105 log 16v 2 should be replaced by 
—856/105 log 16(w^ 2 ECO ) 2 ]. The coefficients cj and c§ are straightforward to compute, but we do not show them 
because they involve rather long expressions. The quantity Wmeco ^ s ^ ne MECO of the energy function ep 2 [defined 
by dep 2 (v)/dv = 0]. The quantity v^^, given by 



"'pole 




35„ 
36 



(55) 



is the pole of ep 2 , which plays an important role in the scheme proposed by Damour, Iyer and Sathyaprakash [26j |. It 
is used to augment the Pade resummation of the PN expanded energy and flux with information taken from the test- 
mass case, where the flux (known analytically up to 5.5PN order) has a pole at the light ring. Under the hypothesis 
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FIG. 6: In the left panel, we plot the energy function £p N versus v, at different PN orders. In the right panel, we plot the 
percentage difference between 2PN and 3PN P-approximants, <5£p(v pea k) = 100 | [£p 3 (u pea k) — £p 2 {vpeak)]/£p 2 (v P oa,k)\ versus 
the total mass M, again evaluated at the LIGO-I peak-sensitivity GW frequency / pca k = 153 Hz [note: t> pca k = (7rAf/ pca k) 1//3 ]. 



N 




(Pn,Pn+i) 






(20 + 5)M 


(10 + 1O)M 


(15 + 15)M Q 


2(9 = +2) 


0.902 


0.915 (0.973, 20.5, 0.242) 


0.868 


2 = -2) 


0.931 


0.955 (0.982, 20.7, 0.236) 


0.923 



TABLE IV: Test for the Cauchy convergence of the P-approximants. The values quoted are maxmax matches obtained by 
maximizing with respect to the extrinsic parameters, but not to the intrinsic parameters (i.e., the matches are computed for 
P waveforms with the same masses, but different PN orders). Here we define P2 = P(2,2.5), P3 = P(3, 3.5). The values in 
parentheses are the maxmax matches obtained by maximizing with respect to the extrinsic and intrinsic parameters, shown 
together with the Pjv+i parameters M and n where the maxima are attained. In all cases the integration of the equations is 
started at a GW frequency of 20 Hz. 

of structural stability |26j| . the flux should have a pole at the light ring also in the comparable-mass case. In the 
test- mass limit, the light ring corresponds to the pole of the energy, so the analytic structure of the flux is modified 
in the comparable-mass case to include v ^(rj). At 3PN order, where the energy has no pole, we choose (somewhat 

arbitrarily) to keep using the value v^ B (rf); the resulting 3PN approximation to the test-mass flux is still very good. 

In Fig. |21 we plot the P-approximants for the flux function Fp N (v), at different PN orders. Note that at 1PN order 
the P-approximant has a pole. At the LIGO-I peak-sensitivity frequency, 153 Hz, for a (10+10)A/ Q BBH, the value 
of v is ~ 0.362, and the percentage difference in Tp N (0.362), between successive PN orders is 1.5PN — > 2PN : —8%; 
2PN -> 2.5PN : +2.2%; 2.5PN -> 3PN = -2) : +3.6%; 3PN -> 3.5PN((9 = -2) : +0.58%. So the percentage 
difference decreases as we increase the PN order. While in the test-mass limit it is known that the P-approximants 
converge quite well to the known exact flux function (see Fig. 3 of Ref. |2(j), in the equal- mass case we cannot be 
sure that the same is happening, because the exact flux function is unknown. (If we assume that the equal-mass flux 
function is a smooth deformation of the test-mass flux function, with 77 the deformation parameter, then we could 
expect that the P-approximants are converging.) In the left panel of Fig. El we plot the P-approximants to the energy 
function as a function of v, at 2PN and 3PN orders; in the right panel, we plot the percentage difference between 
2PN and 3PN P-approximants to the energy function, as a function of the total mass M, evaluated at the LIGO-I 
peak-sensitivity GW frequency /p Ca k = 153 Hz. 

Definition of the models 

When computing the waveforms for P-approximant adiabatic models, the integration of the Eqs. I|32[) is stopped 
at v = i>mecO' which is the solution of the equation d£p N (v)/dv — 0. The corresponding GW frequency will be the 
ending frequency for these waveforms, and in Tab. ITU we show this frequency for typical BBH masses. Henceforth, we 
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FIG. 7: Frequency-domain amplitude versus frequency for the P-approximated (restricted) waveform, at different PN orders, 
for a (15 + 15)M BBH. 



shall refer to the P-approximant models as P(nPN,mPN), and we shall consider (rtPN,mPN) = (2, 2.5), (3, 3.5, 9). 
[Recall that nPN and mPN are the maximum post-Newtonian order of the terms included, respectively, in the energy 
and flux functions £(v) and !F{v)\ at 3PN order we need to indicate also a choice of the arbitrary flux parameter 9.] 



Waveforms and matches 

In Tab. II VI for three typical choices of BBH masses, we perform a convergence test using Cauchy's criterion |2fij . 
The values are quite high, especially if compared to the same test for the T-approximants when the 2.5PN flux is 
used, see Tab. IIIII However, as we already remarked, we do not have a way of testing whether they are converging 
to the true limit. In Fig. we plot the frequency-domain amplitude of the P-approximated (restricted) waveform, 
at different PN orders, for a (15 + 15) Mq BBH. The Newtonian amplitude, -4Newt(/) = f~ 7 ^ 6 , is also shown for 
comparison. At 2.5PN and 3.5PN orders, the evolution is stopped at v = v^ co ; although /meco — 190 — 290 Hz 
(see Tab. ITT)) . the amplitude starts to deviate from /~ 7 / 6 around 100 Hz, well inside the LIGO frequency band. Again, 
this is a consequence of the abrupt termination of the signal in the time domain. 

A quantitative measure of the difference between the P(2, 2.5) and P(3,3.5) waveforms can be seen in Tab. IXII in 
the intersection between the rows and columns labeled P(. . .). For three choices of BBH masses, this table shows the 
maxmax matches between the search models in the columns and the target models in the rows, maximized over the 
search-model parameters M and r), with the restriction < 77 < 1/4. These matches are quite high, but the models 
are not very faithful to each other. The same table shows also the maximized matches (i.e., fitting factors) between T 
and P models. These matches are low between T(2, 2.5) and P(2, 2.5) (and viceversa), between T(2, 2.5) and P(3, 3.5) 
(and viceversa), but they are high between T(2, 2), T(3,3.5) and 3PN P-approximants (although the estimation of 
mass parameters is imprecise). Why this happens can be understood from Fig. [8] by noticing that at 3PN order the 
percentage difference between the T-approximated and P-approximated binding energies is rather small (< 0.5%), 
and that the percentage difference between the T-approximated and P-approximated fluxes at 3PN order (although 
still ~ 10%) is much smaller than at 2PN order. 



IV. NONADIABATIC MODELS 

By contrast with the models discussed in Sec. IIIII in nonadiabatic models we solve equations of motions that involve 
(almost) all the degrees of freedom of the BBH systems. Once again, all waveforms are computed in the restricted 
approximation of Eq. I|29[) . taking the GW phase (y5GW as twice the orbital phase if. 
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FIG. 8: In the left panel, we plot the percentage difference <5£pT(i>peak) = 100 [£p N (« pcak ) - £t n (v pe ak)]/£pjv (%>eak)| versus 
the total mass M, for N = 2, 3, at the LIGO-I peak-sensitivity GW frequency / poa k = 153 Hz [note: n pca k = (7rAf/ pca k) 1//3 ]. 
In the right panel, we plot the percentage difference between 2PN and 3PN P-approximants, 5J-p(v pea k) = 100 | [Tp 3 (w pca k) — 
J r p 2 (vpea,k)]/ '3 r p 2 {v- P esk)\ versus the total mass M, again evaluated at the LIGO-I peak-sensitivity GW frequency / pca k = 153 
Hz. 



A. Nonadiabatic PN expanded methods: Hamiltonian formalism 



Working in the ADM gauge, Damour, Jaranowski and G. Schafer have derived a PN expanded Hamiltonian for the 
general- relativistic two-body dynamics [l(| |4g, |48j : 



H(q, p) = iJ Newt (q, p) + i?iPN(q, p) + #2PN(q, p) + -ff3PN(q, p) , 



(56) 



where 

-ffNcwt (q, p) 
#ipn (q, p) 
#2pn (q, p) 



P_ _ I 

|(3r? - l)(p 2 ) 2 - \ [(3 + r,)p 2 + 77(11 • p) 2 ] J + ± , 

i (1 - 5 V + 5t7 2 ) (p 2 ) 3 + i [(5 - 20 ?? - St? 2 ) (p 2 ) 2 - 2?7 2 (n • p) 2 p 2 - Sr^n • p) 4 ] 1 
16 8 <7 

+\ [(5 + 8?7)p 2 + 3rj(n • p) 2 ] 1 - ~(1 + 3r?)i , 

1 



3\ ^„2n4 



(57) 
(58) 

(59) 



#3PN(q,p) = T ^(-5 + 3 57 7 - 70r; 2 + 3 57 ? a )(p^) 

+1 [(-7 + 427? - 53r? 2 - 5T? 3 ) (p 2 ) 3 + (2 - S^n • p) 2 (p 2 ) 2 + 3(1 - ^(n • p) 4 p 2 - 5rj 3 (n • p) 6 ] - 
16 q 

^ (-27 + 136?? + 10977 2 ) (p 2 ) 2 + 1(17 + 3077)77(11 • p) 2 p 2 + 1(5 + 4377)77(11 ■ p) 4 



+ 



25 
~8~ 
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85 3 2 7 . 

7T 77 Tlln • pf 

16 64 4 ' ' ,v H7 



1 



7T 2 I J] 



(60) 



Here the reduced Non-Relativistic Hamiltonian in the center-of-mass frame, H = iJ NR //x, is written as a function 
of the reduced canonical variables p = pi//i = — P2//J, and q = (xi — X2VM, where Xi and X2 are the positions 
of the BH centers of mass in quasi-Cartesian ADM coordinates (see Refs. 0,EIEI|); the scalars q and p are the 
(coordinate) lengths of the two vectors; and the vector n is just q/q. 
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Equations of motion 

We now restrict the motion to a plane, and we introduce radiation-reaction (RR) effects as in Ref. 15). The 
equations of motion then read (using polar coordinates r and <p obtained from the q with the usual Cartesian-to-polar 
transformation) 

dr dH dip _^ dH 

-~= ^— {r,Pr,P<p), —~:=u) = - — {r,p r ,p v ), (61) 

dt Op r V dt Op v 

~7F = -^r-{r,Pr,P v ) + F r {r,p r ,p v ) , = F v [£D(r,p r ,p v )] , (62) 

dt or dt 

where t — t/M, Q — ujM; and where F v = F v j fi and F r = F r / ' [i are the reduced angular and radial components of 
the RR force. Assuming F r <C F v [l^, averaging over an orbit, and using the balance equation l|31f) . we can express 
the angular component of the radiation-reaction force in terms of the GW flux at infinity [l5| . More explicitly, if we 
use the P-approximated flux, we have 

pv = f Pn K\ = 3 % K = - if »7 u w , p 2 , : , 63) 

while if we use the T-approximated flux we have 

= F Tn [v u ] ^ T Tn [v u ] , (64) 

where v w = a) 1 / 3 = (dip/dt) 1 / 3 . This w w is used in Eq. (|29|l to compute the restricted waveform. Note that at each PN 
order, say nPN, we define our Hamiltonian model by evolving the Eqs. (|61|l and (|62|) without truncating the partial 
derivatives at the nPN order (differentiation with respect to the canonical variables can introduce terms of order 
higher than nPN). Because of this choice, and because of the approximation used to incorporate radiation- reaction 
effects, these nonadiabatic models are not, strictly speaking, purely post-Newtonian. 



Innermost stable circular orbit 



Circular orbits are defined by setting r = constant while neglecting radiation- reaction effects. In our PN Hamiltonian 
models, this implies dH/dp r = through Eq. (|61(l : because at all PN orders the Hamiltonian H [Eqs. H56JI - H60J1 ] ^ s 
quadratic in p r , this condition is satisfied for p r = 0; in turn, this implies also dH /dr = [through Eq. H62J) ]. which 
can be solved for p v . The orbital frequency is then given by Cj = dH/dp v . 

The stability of circular orbits under radial perturbations depends on the second derivative of the Hamiltonian: 

> <^4> stable orbit ; < unstable orbit . (65) 



dr 2 dr 2 

For a test particle in Schwarzschild geometry (the i] — > of a BBH), an Innermost Stable Circular Orbit (ISCO) 
always exists, and it is defined by 

nrrSchw n2 rrSchw 

ISCO (Schwarzschild) : — = —— 2 — = 0, (66) 

® r \p r =a or \p T =0 

where 7J Schw (r, p r ,P v ) is the (reduced) nonrelativistic test-particle Hamiltonian in the Schwarzschild geometry. Sim- 
ilarly, if such an ISCO exists for the (reduced) nonrelativistic PN Hamiltonian H [Eq. it is defined by 

f)M f) 2 M 

ISCO (Hamiltonian) : — = — - = 0. (67) 

ot |j> r =o or \ Pr =o 

Any inspiral built as an adiabatic sequence of quasicircular orbits cannot be extended to orbital separations smaller 
than the ISCO. In our model, we integrate the Hamiltonian equations (|61|l and l|62|l including terms up to a given PN 
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order, without re-truncating the equations to exclude terms of higher order that have been generated by differentiation 
with respect to the canonical variables. Consistently, the value of the 1SCO that is relevant to our model should be 
derived by solving Eq. 1)67(1 without any further PN truncation. 

How is the ISCO related to the Maximum binding Energy for Circular Orbit (MECO), used above for nonadiabatic 
models such as T? The PN expanded energy for circular orbits £t„ (p) at order nPN can be recovered by solving the 
equations 



dH(r,p r = 0,p v ) 



0. 



dH(r,p r = 0,p v ) 



dr dp v 
for r and p v as functions of Co, and by using the solutions to define 

H(Co)=H[r(Co),p r = 0,p v (Co)}. 



(68) 



(69) 



Then H(Co 



= „3\ - 



£t„(w) as given by Eq. (133(1 . if and only if in this procedure we are careful to eliminate all terms 



of order higher than nPN (see, e.g., Ref. 54j). 

In the context of nonadiabatic models, the MECO is then defined by 



MECO 



dH 

duo 



0: 



(70) 



and it also characterizes the end of adiabatic sequences of circular orbits. Computing the variation of Eq. I|69|) between 
nearby circular orbits, and setting p r — 0, dp r = 0, we get 



duj = TT^ — dr 
drdp v 

and combining these two equations we get 

dp v 
duo 



d 2 H 
dpi 



d 2 H 
dr 2 



dp 



v > 



d 2 H 
dr 2 



dr 



d 2 H 
drdp v 



dp v = : 



(71) 



d 2 H 
drdp v 



d 2 H d 2 H 
dp 2 dr 2 



(72) 



So finally we can write 



dH 

duj 



dH dp v _ d 2 H dH 



dp v du) 



dr 2 dp^ 



d 2 H 
drdp v 



d 2 H d 2 H 



dpi 



dr 2 



(73) 



Not surprisingly, Eqs. (|73|) and (|69fl together are formally equivalent to the definition of the ISCO, Eq. (|67|) [note that 
the second and third terms on the right-hand side of Eq. I|73|) are never zero.] Therefore, if we knew the Hamiltonian 
H exactly, we would find that the MECO defined by Eq. J7UJ), is numerically the same as the ISCO defined by 
Eq. 1671) . Unfortunately, we are working only up to a finite PN order (say nPN); thus, to recover the MECO as given 
by Eq. 1)33(1 . all three terms on the right-hand side of Eq. ((73(1 must be written in terms of lj, truncated at nPN order, 
then combined and truncated again at nPN order. This value of the MECO, however, will no longer be the same as 
the ISCO obtained by solving Eq. ((67|) exactly without truncation. 

If the PN expansion was converging rapidly, then the difference between the ISCO and the MECO would be mild; 
but for the range of BH masses that we consider the PN convergence is bad, and the discrepancy is rather important. 
The ISCO is present only at 1PN order, with nsco = 9.907 and tDisco = 0.02833. The corresponding GW frequencies 



are given in Tab. HTlfor a few BBHs with equal masses. At 3PN order we find the formal solution 



„isco _ 



1.033 and 



piSCO _ 0.355, but since we do not trust the PN expanded Hamiltonian when the radial coordinate gets so small, we 
conclude that there is no ISCO at 3PN order. 



Definition of the models 



In order to build a quasicircular orbit with initial GW frequency /o, our initial conditions (nmtjPrini^JVinit) are 
set by imposing cpinit = 7r/o, p r init = and dr- ln it/dt = — F / (r/dH / dr) c i IC , as in Ref. |56j. The initial orbital phase 
(^init remains a free parameter. For these models, the criterion used to stop the integration of Eqs. I(61|) . I(62|) is rather 
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N 


(5 + 2O)M 


{HTjv,HT]v+i) 
(10 + 1O)M 


(15 + 15)M 


(5 + 2O)M 


(HPjv,HPjv+i) 
(10 + 10)M o 


(15 + 15)M 





0.118 


0.191 (0.553, 13.7, 0.243) 


0.206 


0.253 


0.431 (0.586, 16.7, 0.242) 


0.316 


1 


0.102 


0.174 (0.643, 61.0, 0.240) 


0.170 


0.096 


0.161 (0.623, 17.4, 0.239) 


0.151 


2 (0 = +2) 


0.292 


0.476 (0.656, 18.6, 0.241) 


0.377 


0.266 


0.369 (0.618, 17.6, 0.240) 


0.325 


2 (0 = -2) 


0.287 


0.431 (0.671, 19.0, 0.241) 


0.377 


0.252 


0.354 (0.622, 17.4, 0.239) 


0.312 



TABLE V: Test for the Cauchy convergence of the HT- and HP-approximants. The values quoted are maxmax matches 
obtained by maximizing with respect to the extrinsic parameters, but not to the intrinsic parameters (i.e., the matches are 
computed for H waveforms with the same masses, but different PN orders). Here we define HTo = HT(0, 0), HTi = HT(1, 1.5), 
HT 2 = HT(2, 2) [because the 2.5PN flux goes to zero before the MECO is reached, so we use the 2PN flux], HT 3 = HT(3, 3.5, 0); 
we also define HP = HP(0,0), HPi = HP(1, 1.5), HP 2 = HP(2,2.5), and HP 3 = HP(3, 3.5,6). The values in parentheses are 
the maxmax matches obtained by maximizing with respect to the extrinsic and intrinsic parameters, shown together with the 
Hjv+i parameters M and n where the maxima are attained. In all cases the integration of the equations is started at a GW 
frequency of 20 Hz. 
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FIG. 9: Inspiraling orbits in (x, y)-plane when r\ — 0.25 for HT(1, 1.5) (in the left panel) and HT(3, 3.5, 0) (in the right panel). 
For a (15 + 15)M BBH the evolution starts at /gw = 34 Hz and ends at /gw = 97 Hz for HT(1, 1.5) panel and at /gw = 447 
Hz for HT(3, 3.5, 0). The dynamical evolution is rather different because at 1PN order there is an ISCO (nsco — 9.9M), while 
at 3PN order it does not exist. 



arbitrary. We decided to push the integration of the dynamical equations up to the time when we begin to observe 
unphysical effects due to the failure of the PN expansion, or when the assumptions that underlie Eqs. (|62[) [such 
as F r <C F 9 ], cease to be valid. When the 2.5PN flux is used, we stop the integration when Tt n equals 10% of 
•^Newt, and we define the ending frequency for these waveforms as the instantaneous GW frequency at that time. To 
be consistent with the assumption of quasicircular motion, we require also that the radial velocity be always much 
smaller than the orbital velocity, and we stop the integration when |r| > 0.3(r<^), if this occurs before Tt n equals 
10% of JFNcwt- In some cases, during the last stages of inspiral Q reaches a maximum and then drops quickly to zero 
[see discussion in Sec.|V]. When this happens, we stop the evolution at u; = 0. 

We shall refer to these models as HT(nPN, toPN) (when the T-approximant is used for the flux) or HP(nPN, mPN) 
(when the P-approximant is used for the flux), where nPN (mPN) denotes the maximum PN order of the terms 
included in the Hamiltonian (the flux). We shall consider (nPN, mPN) = (1, 1.5), (2, 2), (2, 2.5), and (3,3.5,0) [at 
3PN order we need to indicate also a choice of the arbitrary flux parameter 9} . 
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FIG. 10: Frequency-domain amplitude versus frequency for the HT and HP (restricted) waveforms, at different PN orders, for 
a (15 + 15) Af© BBH. The HT(3, 3.5, = 0) curve, not plotted, is almost identical to the HT(3, 3,9 = 0) curve. 



Waveforms and matches 

In Tab. Ivl for three typical choices of BBH masses, we perform a convergence test using Cauchy's criterion |2fij . 
The values are very low. For N = and N = 1, the low values are explained by the fact that at 1PN order there 
is an ISCO [see the discussion below Eq. (|73J) ] . while at Newtonian and 2PN, 3PN order there is not. Because of 
the ISCO, the stopping criterion [|f| > 0.3(rtp) or Q = 0] is satisfied at a much lower frequency, hence at 1PN order 
the evolution ends much earlier than in the Newtonian and 2PN order cases. In Fig. |5J we show the inspiraling 
orbits in the (x,y) plane for equal-mass BBHs, computed using the HT(1,1.5) model (in the left panel) and the 
HT(3, 3.5, 0) model (in the right panel). For N = 2, the low values are due mainly to differences in the conservative 
dynamics, that is, to differences between the 2PN and 3PN Hamiltonians. Indeed, for a (10 + 10)M Q BBH we find 
(HT(2,2),HT(3,2)) = 0.396, still low, while (HT(2, 2), HT(2, 3.5)) = 0.662, considerably higher than the values in 
Tab. E| 

In Fig. 1101 we plot the frequency-domain amplitude of the HT-approximated (restricted) waveforms, at different 
PN orders, for a (15+ 15)M Q BBH. The Newtonian amplitude, ^4.Newt(/) = /~ 7 ^ 6 , is also shown for comparison. For 
HT(1, 1.5), because the ISCO is at r ~ 9.9M, the stopping criterion |r| > 0.3 ip r is reached at a very low frequency 
and the amplitude deviates from the Newtonian prediction already at / ~ 50Hz. For HT(2,2.5), the integration of 
the dynamical equation is stopped as the flux function goes to zero; just before this happens, the RR effects become 
weaker and weaker, and in the absence of an ISCO the two BHs do not plunge, but continue on a quasicircular orbit 
until Ft{v) equals 10% of .T-Ncwt- So the binary spends many cycles at high frequencies, skewing the amplitude with 
respect to the Newtonian result, and producing the oscillations seen in Fig. 1101 We consider this behaviour rather 
unphysical, and in the following we shall no longer take into account the HT(2, 2.5) model, but at 2PN order we shall 
useHT(2,2). 

The situation is similar for the HP models. Except at 1PN order, the HT and HP models do not end their evolution 
with a plunge. As a result, the frequency-domain amplitude of the HT and HP waveforms does not decrease markedly 
at high frequencies, as seen in Fig. 1101 and in fact it does not deviate much from the Newtonian result (especially at 
3PN order). 

Quantitative measures of the difference between HT and HP models at 2PN and 3PN orders, and of the difference 
between the Hamiltonian models and the adiabatic models, can be seen in Tables EH EH For some choices of BBH 
masses, these tables show the maxmax matches between the search models in the columns and the target models in 
the rows, maximized over the search-model parameters M and 77, with the restriction < 77 < 1/4. The matches 
between the H(2, 2) and the H(3, 3.5) waveforms are surprisingly low. More generally, the H(2, 2) models have low 
matches with all the other PN models. We consider these facts as an indication of the unreliability of the H models. 
In the following we shall not give much credit to the H(2, 2) models, and when we discuss the construction of detection 
template families we shall consider only the H(3, 3.5) models. [We will however comment on the projection of the 
H(2, 2) models onto the detection template space.] 

As for the H(3,3.5) models, their matches with the 2PN adiabatic models are low; but their matches with the 3PN 
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adiabatic models are high, at least for M < 30M©. For M — 40M© (as shown in Tables IXll and |XTT|) . the matches 
can be quite low, as the differences in the late dynamical evolution become significant. 



B. Nonadiabatic PN expanded methods: Lagrangian formalism 

Equations of motion 



In the harmonic gauge, the equations of motion for the general-relativistic two-body dynamics in the Lagrangian 
formalism read I45U57L 581: 



where 



X — aN + apN + a 2 PN + &2.5RR + a 3.5RR ; 



M 



a N = --^-n, 



(74) 
(75) 
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For the sake of convenience, in this section we are using same symbols of Sec. IIV Al to denote different physical 
quantities (such as coordinates in different gauges). Here the vector x = xi — X2 is the difference, in pseudo- 
Cartesian harmonic coordinates ^5|, between the positions of the BH centers of mass; the vector v = dx/dt is the 
corresponding velocity; the scalar r is the (coordinate) length of x; the vector n = x/r; and overdots denote time 
derivatives with respect to the post-Newtonian time. We have included neither the 3PN order corrections 03pn 
derived in Ref. nor the 4.5PN order term a4.5PN for the radiation-reaction force computed in Ref. [5^|. Unlike 
the Hamiltonian models, where the radiation-reaction effects were averaged over circular orbits but were present up 
to 3PN order, here radiation-reaction effects are instantaneous, and can be used to compute generic orbits, but are 
given only up to 1PN order beyond the leading quadrupole term. 

We compute waveforms in the quadrupole approximation of Eq. I|29|l . defining the orbital phase ip as the angle 
between x and a fixed direction in the orbital plane, and the invariantly defined velocity v as (Mip) 1 / 3 . 



Definition of the models 



For these models, just as for the HT and HP models, the choice of the endpoint of evolution is rather arbitrary. 
We decided to stop the integration of the dynamical equations when we begin to observe unphysical effects due 
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FIG. 11: Frequency-domain amplitude versus frequency for the L-approximated (restricted) waveforms, at different PN orders, 
for a (15 + 15)M BBH. 





T(2, 2) T(3, 3.5, 0) P(2, 2.5) P(3, 3.5, 0) EP(2, 2.5) EP(3, 3.5, 0) HT(3, 3.5, 0) 
mm M 7] mm M rj mm M i] mm M rj mm M rj mm M r} mm M rj 


L(2, 1) (20+20)M o 
(15+15)M 
(15+5)M Q 
(5+5)M Q 


0.994 78.83 0.05 
0.991 55.16 0.06 
0.981 35.51 0.05 
0.956 10.68 0.20 


0.998 61.24 0.09 
0.995 44.50 0.10 
0.991 29.03 0.08 
0.965 11.49 0.18 


0.999 52.76 0.13 
0.999 39.96 0.13 
0.995 26.02 0.10 
0.971 11.33 0.19 


0.998 57.96 0.11 
0.998 43.57 0.11 
0.994 27.99 0.09 
0.964 11.89 0.17 


0.935 70.76 0.05 
0.912 46.67 0.09 
0.942 27.46 0.09 
0.964 11.03 0.19 


0.944672.04 0.06 
0.916 50.90 0.07 
0.941 28.85 0.08 
0.960 11.69 0.17 


0.994 49.53 0.14 
0.994 37.08 0.15 
0.994 22.89 0.13 
0.966 11.32 0.18 



TABLE VI: Fitting factors [see Eq. H2U1 1 for the projection of the L(2, 1) (target) waveforms onto the T, P, EP and HP (search) 
models at 2PN and 3PN order. The values quoted are obtained by maximizing the maxmax (mm) match over the search-model 
parameters M and n. 



to the failure of the PN expansion. For many (if not all) configurations, the PN-expanded center-of-mass binding 
energy (given by Eqs. (2.7a)-(2.7e) of Ref. [l!|) begins to increase during the late inspiral, instead of continuing to 
decrease. When this happens, we stop the integration. The instantaneous GW frequency at that time will then be the 
ending frequency for these waveforms. We shall refer to these models as L(nPN,mPN), where nPN (nPN) denotes 
the maximum PN order of the terms included in the Hamiltonian (the radiation-reaction force). We shall consider 
(nPN,mPN) = (2,0), (2,1). 



Waveforms and matches 

In Fig. 1111 we plot the frequency-domain amplitude versus frequency for the L-approximated (restricted) waveforms, 
at different PN orders, for a (15 + 15)A/0 BBH. The amplitude deviates from the Newtonian prediction slightly before 
100 Hz. Indeed, the GW ending frequencies are 116 Hz and 107 Hz for the L(2, 0) and L(2,l) models, respectively. 
These frequencies are quite low, because the unphysical behavior of the PN-expanded center-of-mass binding energy 
appears quite early [at r onc i = 6.6 and r cn( j = 7.0 for the L(2, 0) and L(2, 1) models, respectively]. So the L models do 
not provide waveforms for the last stage of inspirals and plunge. 

Tablc lVll shows the maxmax matches between the L-approximants and a few other selected PN models. The overlaps 
are quite high, except with the EP(2,2.5) and EP(3, 3.5, 0) at high masses, but extremely unfaithful. Moreover, we 
could expect the L(2, 0) and L(2, 1) models to have high fitting factors with the adiabatic models T(2, 0) and T(2, 1). 
However, this is not the case. As Table IV 111 shows, the T models are neither effectual nor faithful in matching the 
L models, and vice versa. This might be due to one of the following factors: (i) the PN-expanded conservative 
dynamics in the adiabatic limit (T models) and in the nonadiabatic case (L models) are rather different; (ii) there is 
an important effect due to the different criteria used to end the evolution in the two models, which make the ending 
frequencies rather different. All in all, the L models do not seem very reliable, so we shall not give them much credit 
when we discuss detection template families. However, we shall investigate where they lie in the detection template 
space. 
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L(2,0) T(2,0) L(2,l) T(2, 1) 

mm M r\ mm M r\ mm M r\ mm M r\ 



(15+15)M 
L(2,0) (15+5)M 
(5+5)M Q 




0.884 42.02 0.237 
0.769 24.71 0.201 
0.996 21.70 0.068 






(15+15)M 
T(2,0) (15+5)M 
(5+5)M 


0.834 23.44 0.247 
0.823 14.90 0.250 
0.745 9.11 0.250 








(15+15)M 
L(2,l) (15+5)M 
(5+5)M 








0.837 60.52 0.236 
0.844 55.70 0.052 
0.626 11.47 0.238 


(15+15)M 
T(2, 1) (15+5)M 
(5+5)M 






0.663 19.38 0.250 
0.672 13.56 0.250 
0.631 9.22 0.243 





TABLE VII: Fitting factors [see Eq. Pty] for the projection of the L(2, 1) and L(2, 0) (target) waveforms onto the T(2, 0) 
and T(2, 1) (search) models. The values quoted are obtained by maximizing the maxmax (mm) match over the search-model 
parameters M and rj. 



C. Nonadiabatic PN resummed methods: the Effective-One-Body approach 



The basic idea of the effective-one-body (EOB) approach 0] is to map the real two-body conservative dynamics, 
generated by the Hamiltonian l|56|) and specified up to 3PN order, onto an effective one-body problem where a test 
particle of mass /i = mim2/M (with m\ and TO2 the BH masses, and M = mi + TO2) moves in an effective background 
metric g e ^, given by 



ds'i s = gf„ dx" dx v = —A(R) c l dt 



D{R) 



dR 2 +R 2 {d0 2 + sin 2 6 dip 2 



where 



A(R) 

, GM ( GM \ 2 ( GM\ 3 (GM\ 
A{R) = 1 + a ^ +a2 {-^R) +03 Ui?J +ai {^R) 



CM 

D(R) = l + d 1 — + d 2 



(GM\ (GM\ 
<*R ' \^R) +d3 \^R) 



The motion of the particle is described by the action 



S e fi 



fj,c / ds c s ■ 



(80) 

(81) 
(82) 

(83) 



For the sake of convenience, in this section we shall use the same symbols of Sees. II V Al and HVBI to denote different 
physical quantities (such as coordinates in different gauges). The mapping between the real and the effective dynamics 
is worked out within the Hamilton- Jacobi formalism, by imposing that the action variables of the real and effective 
description coincide (i.e., J loa i = J c ff, T Tea ,i = ^effj where J denotes the total angular momentum, and X the radial 
action variable 14j), while allowing the energy to change, 



cNR 
c eff 

flC 2 



cNR. 
'-'real 

/iC 2 
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„ c rcal 1 „ 

\xc z 



cNR 
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jJLC 2 



OLZ 



cNR. 
c rcal 

lie 2 



(84) 



here is the Non-Relativistic effective energy, while is related to the relativistic effective energy £ c h by the equation 
£^p- — £ c ff — fic 2 ; £ c ff is itself defined uniquely by the action (|83|l . The Non-Relativistic real energy £^^ = i?(q, p), 
where H (q, p) is given by Eq. (j56(l with H(q, p) = fiH(q,p). From now on, we shall relax our notation and set 
G = c = 1. 



Equations of motion 



Damour, Jaranowski and Schafer [l6l | found that, at 3PN order, this matching procedure contains more equations 
to satisfy than free parameters to solve for (ai, 0,2, 03, di, d-2, efa, and a±, a 2l ota). These authors suggested the following 
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two solutions to this conundrum. At the price of modifying the energy map and the coefficients of the effective metric 
at the fPN and 2PN levels, it is still possible at 3PN order to map uniquely the real two-body dynamics onto the 
dynamics of a test mass moving on a geodesic (for details, see App. A of Ref. HE])- However, this solution appears 
very complicated; more importantly, it seems awkward to have to compute the 3PN Hamiltonian as a foundation for 
deriving the matching at the 1PN and 2PN levels. The second solution is to abandon the hypothesis that the effective 
test mass moves along a geodesic, and to augment the Hamilton-Jacobi equation with (arbitrary) higher-derivative 
terms that provide enough coefficients to complete the matching. With this procedure, the Hamilton-Jacobi equation 
reads 



= + (x) fv p v + A mpa (x) ^ Pri p p p (j + 



(85) 



Because of the quartic terms A a/3l,s , the effective 3PN relativistic Hamiltonian is not uniquely fixed by the matching 
rules defined above; the general expression is |l6j |: 



cNR _ 
C eS = 



#eff(q,p) - \ A(q) 



1 



(i|- 1 ) (n ' p)2 + ? <2 ' ,p3)2 



z 2 p 2 (n • p) 2 + z 3 (n • p) 



(86) 



here we use the reduced relativistic effective Hamiltonian H e g = H g/fj,, and q and p are the reduced canonical 

variables, obtained by rescaling the canonical variables by M and /i,, respectively. The coefficients 21,22 and 23 are 
arbitrary, subject to the constraint 

8z!+4z2 + 3z 3 = 6(4-3??) 77. (87) 

Moreover, we slightly modify the EOB model at 3PN order of Ref. [HI by requiring that in the test mass limit the 3PN 
EOB Hamiltonian equal the Schwarzschild Hamiltonian. Indeed, one of the original rationales of the PN resummation 
methods was to recover known exact results in the test-mass limit. To achieve this, Z\, z 2 and 23 must go to zero as 
77 — > 0. A simple way to enforce this limit is to set z\ = rjzi, z 2 = i]z 2 and 23 = 7723. With this choice the coefficients 
A(r) and D(r) in Eq. (jHSJ) read: 



A(r) 
D(r) 



1 
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2rj 
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41 , 
32' 



V 
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(89) 



where we set r — |q|. The authors of Ref. ^(| restricted themselves to the case z\ = z 2 = (21 = 22 = 0). Indeed, 
they observed that for quasicircular orbits the terms proportional to 22 and 23 in Eq. (|86|l are very small, while for 
circular orbits the term proportional to z\ contributes to the coefficient A(r), as seen in Eq. (|88|l . So, if the coefficient 
z\ = r\z\ + 1 0, its value could be chosen such as to cancel the 3PN contribution in A(r). To avoid this fact, which can 
be also thought as a gauge effect due to the choice of the coordinate system in the effective description, the authors of 
Ref. ^(| decided to pose z\ = (21 = 0). By contrast, in this paper we prefer to explore the effect of having z\, 2 + 1 0. 
So we shall depart from the general philosophy followed by the authors in Ref. , pushing (or expanding) the EOB 
approach to more extreme regimes. 

Now, the reduction to the one-body dynamics fixes the arbitrary coefficients in Eq. I|84|l uniquely to a± = 77/2, 
a 2 = 0, and 03 = 0, and provides the resummed (improved) Hamiltonian [obtained by solving for m -^q- l|84|) 
and imposing H im P TOVcd = £NRj . 



H 



improved 



M 4/1 + 277 



-Heff - H 



(90) 



Including radiation-reaction effects, we can then write the Hamilton equations in terms of the reduced quantities 

^improved = ^improved/^ f = t j M ^ ^ = [3, 



dr ^improved 

dt 

dip 
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(91) 

(92) 

(93) 
(94) 
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where for the if component of the radiation-reaction force we use the T- and P-approximants to the flux function [see 
Eqs. (E3J]. Note that at each PN order, say nPN, we integrate the Eqs. (|!JT)l - lj$4l without further truncating 

the partial derivatives of the Hamiltonian at nPN order (differentiation with respect to the canonical variables can 
introduce terms of order higher than nPN). 

Following the discussion around Eq. I|67[l . the ISCO of these models is determined by setting dHQ llpiovcd /dr — 
d*H™ pmved /dr 2 = 0, where #j roproved (r,Pr,2V) = Approved (r, 0,p v ). If we define 

H 2 s (r,0,p v ) = W Pv = A(r) U + ^+tjzt^J , (95) 

we extract the ISCO by imposing dW p (r)/dr = = d 2 W Pv (r)/d 2 r. Damour, Jaranowski and Schafer ^(| noticed 
that at 3PN order, for Z\ — z 2 — 0, and using the PN expanded form for A{r) given by Eq. Q88JI. there is no ISCO. 
To improve the behavior of the PN expansion of A(r) and introduce an ISCO, they proposed replacing A(r) with the 
Pade approximants 

ApJr) = - , (96) 

2W 2r 2 + 2rj + rrj ' y ' 
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(98) 



In Table [HJ we show the GW frequency at the ISCO for some typical choices of BBH masses, computed using the 
above expressions for A(r) in the improved Hamiltonian Q90|l with z\ — z 2 = 0. 

We use the Pade resummation for A(r) of Ref. 0] also for the general case z\ ^ 0, because for the PN expanded 
form of A(r) the ISCO does not exist for a wide range of values of z\. [However, when we discuss Fourier-domain 
detection template families in Sec. I VII we shall investigate also EOB models with PN-expanded ^4(r).] 

In Fig. El we plot the binding energy as evaluated using the improved Hamiltonian 190J1 . at different PN orders, 
for equal-mass BBHs. At 3PN order, we use as typical values z\ = 0, ±4. [For z\ > 4 the location of the ISCO is 
no longer a monotonic function of z\. So we set z.\ < 4.] In the right panel of Fig. 1121 we show the variation in the 
GW frequency at the ISCO as a function of Si for a (15+15)M Q BBH. Finally, in Fig.^J we compare the binding 
energy for a few selected PN models, where for the E models we fix z\ = z 2 — [see the left panel of Fig. ^| for the 
dependence of the binding energy on the coefficient z%\. Notice, in the left panel, that the 2PN and 3PN T energies are 
much closer to each other than the 2PN and 3PN P energies are, and than the 2PN and 3PN E energies are; notice also 
that the 3PN T and P energies are very close. The closeness of the binding energies (and of the MECOs and ISCOs) 
predicted by PN expanded and resummed models at 3PN order (with z\ =0), and of the bindin g en ergy predicted 
by the numerical quasiequilibrium BBH models of Ref. [25| was recently pointed out in Refs. |54l l65|. However, the 
EOB results are very close to the numerical results of Ref. only if the range of variation of z\ is restricted. 



Definition of the models 

For these models, we use the initial conditions laid down in Ref. an d also adopted in this paper for the HT 
and HP models (see Sec. IIV A"|) . At 2PN order, we stop the integration of the Hamilton equations at the light ring 
given by the solution of the equation r 3 — 3r 2 + 5r/ = 15]. At 3PN order, the light ring is defined by the solution of 

±[u 2 A P3 (u)} =0, (99) 

with u = 1/r and Ap 3 is given by Eq. I|97|) . For some configurations, the orbital frequency and the binding energy 
start to decrease before the binary can reach the 3PN light ring, so we stop the evolution when lj — [see discussion 
in Sec. IIV Dl below]. For other configurations, it happens that the radial velocity becomes comparable to the angular 
velocity before the binary reaches the light ring; in this case, the approximation used to introduce the RR effects into 
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FIG. 12: In the left panel, we plot the binding energy evaluated using the improved Hamiltonian l|9()^ . as a function of the 
velocity parameter v, for equal-mass BBHs, n — 0.25. We plot different PN orders for the E-model varying also the parameter 
Zi. In the right panel, we plot the GW frequency at the ISCO at 3PN order as a function of the parameter zi for (15+15)M0 
BBH. 
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FIG. 13: Binding energy as a function of the velocity parameter v, for equal-mass BBHs. We plot different PN orders for 
selected PN models. For the E model at 3PN order we fix z\ — = Z2- 



the conservative dynamics is no longer valid, and we stop the integration of the Hamilton equations when \f/(rip)\ 
reaches 0.3. For some models, usually those with 7^ 0, the quantity \r/(np)\ reaches a maximum during the 
last stages of evolution, then it starts decreasing, and f becomes positive. In such cases, we choose to stop at the 
maximum of \f/(rip)\. 

In any of these cases, the instantaneous GW frequency at the time when the integration is stopped defines the 
ending frequency for these waveforms. 

We shall refer to the EOB models (E-approximants) as ET(nPN, toPN) (when the T-approximant is used for the 
flux) or EP(nPN,mPN) (when the P-approximant is used for the flux), where rtPN (mPN) denotes the maximum 
PN order of the terms included in the Hamiltonian (flux). We shall consider (nPN, mPN) = (1,1.5), (2,2.5), and 
(3, 3.5, 0) [at 3PN order we need to indicate also a choice of the arbitrary flux parameter 9]. 
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Waveforms and matches 

In Table IXIVI we investigate the dependence of the E waveforms on the values of the unknown parameters z.\ 
and Z2 that appear in the EOB Hamiltonian at 3PN order. The coefficients Z\ and z-i are in principle completely 
arbitrary. When Z\ ^ 0, the location of the ISCO changes, as shown in Fig. ^] Moreover, because in Eq. (|86|1 z\ 
multiplies a term that is not zero on circular orbits, the motion tends to become noncircular much earlier, and the 
criteria for ending the integration of the Hamilton equations are satisfied earlier. [See the discussion of the ending 
frequency in the previous section.] This effect is much stronger in equal-mass BBHs with high M . For example, for 
(15+15)Ajf BBHs and for £2 = 0, the fitting factor (the maxmax match, maximized over M and 77) between an EP 
target waveform with z± = and EP search waveforms with —40 < z± < —4 can well be < 0.9. However, if we restrict 
Z\ to the range [—4,4], we get very high fitting factors, as shown in Table IXTVl 

In Eq. (|86J) . the coefficients 2,2 and £3 multiply terms that are zero on circular orbits. [The coefficient £2 appears 
also in D(r), given by Eq. (|89|l .] So their effect on the dynamics is not very important, as confirmed by the very high 
matches obtained in Table IXTVl between EP waveforms with S2 = and EP waveforms with £2 = ±4. It seems that 
the effect of changing z 2 is nearly the same as a remapping of the BBH mass parameters. 

We investigated also the case in which we use the PN expanded form for A(r) given by Eq. (|88|) . For example, 
for (15+15)M BBHs and £2 — 0, the fitting factors between EP target waveforms with Z\ = —40,-4,4,40 and EP 
search waveforms with z x = are (maxmax, M, 77) = (0.767,39.55,0.240), (0.993,30.83,0.241), (0.970,30.03,0.241), 
and (0.915,28.23,0.242), respectively. So the overlaps can be quite low. 

In Table IVIIII for three typical choices of BBH masses, we perform a convergence test using Cauchy's criterion. 
The values are quite high. However, as for the P-approximants, we have no way to test whether the E-approximants 
are converging to the true limit. In Fig. 1141 we plot the frequency-domain amplitude of the EP-approximated 
(restricted) waveforms, at different PN orders, for a (15 +15)Af BBH. The evolution of the EOB models contains a 
plunge characterized by quasicircular motion [15j. This plunge causes the amplitude to deviate from the Newtonian 
amplitude, y^Nowt = /~ 7 ^ 6 around 200 Hz, which is a higher frequency than we found for the adiabatic models [see 

Figs. mm. 

In Table IXTTTl for some typical choices of the masses, we evaluate the fitting factors between the ET(2,2.5) and 
ET(3, 3.5) waveforms (with z\ — £2 — 0) and the T(2, 2.5) and T(3, 3.5) waveforms. This comparison should emphasize 
the effect of moving from the adiabatic orbital evolution, ruled by the energy-balance equation, to the (almost) full 
Hamiltonian dynamics, ruled by the Hamilton equations. More specifically, we see the effect of the differences in the 
conservative dynamics between the PN expanded T-model and the PN resummed E-model (the radiation-reaction 
effects are introduced in the same way in both models). While the matches are quite low at 2PN order, they are high 
(> 0.95) at 3PN order, at least for M < 3OAf , but the estimation of mi and 1112 is poor. This result suggests that, 
for the purpose of signal detection as opposed to parameter estimation, the conservative dynamics predicted by the 
EOB resummation and by the PN expansion are very close at 3PN order, at least for M < 3OM . Moreover, the 
results of Table IXIIII suggest also that the effect of the unknown parameter 9 is rather small, at least if 9 is of order 
unity, so in the following we shall always set 9 — 0. 

In Tables IXT1 and IXTTl we study the difference between the EP(2,2.5) and EP(3,3.5) models (with £\ — Z2=0), and 
all the other adiabatic and nonadiabatic models. For some choices of BBH masses, these tables show the maxmax 
matches between the search models in the columns and the target models in the rows, maximized over the search- 
model parameters M and 77, with the restriction < 77 < 1/4. At 2PN order, the matches with the T(2, 2.5), HT(2, 2) 
and HP(2,2.5) models are low, while with the matches with the T(2,2) and P(2,2.5) models are high, at least for 
M < 3OM (but the estimation of the BH masses is poor). At 3PN order, the matches with T(3, 3.5,9), P(3, 3.5,9), 
HP(3, 3.5, 9) and HT(3, 3.5, 9) are quite high if M < 30M Q . However, for M = 4OM , the matches can be quite low. 
We expect that this happens because in this latter case the differences in the late dynamical evolution become crucial. 

D. Features of the late dynamical evolution in nonadiabatic models 

While studying the numerical evolution of nonadiabatic models, we encounter two kinds of dynamical behavior that 
are inconsistent with the assumption of quasicircular motion used to include the radiation-reaction effects, so when 
one of these two behaviors occurs, we immediately stop the integration of the equations of motion. First, in the late 
stage of evolution Q can reach a maximum, and then drop quickly to zero; so we stop the integration if Q = 0. Second, 
the radial velocity r can become a significant portion of the total speed, so we stop the integration if r = 0.3(rQ). 

The first behavior is found mainly in the H models at 3PN order, when 77 is relatively small (< 0.21). As we shall see 
below, it is not characteristic of either the Schwarzschild Hamiltonian or the EOB Hamiltonian. In the left panel of 
Fig. El we plot the binding energy evaluated from H(r,p r = 0,p v ) [given by Eq. JSHJ] as a function of r at 77 = 0.16, 
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N 


(5 + 20) M Q 


{EP]v,EPjv+i} 
(10 + 1O)M 


(15 + 15)M 





0.677 


0.584 (0.769, 17.4, 0.246) 


0.811 


1 


0.766 


0.771 (0.999, 21.8, 0.218) 


0.871 


2 = +2) 


0.862 


0.858 (0.999, 21.3, 0.222) 


0.898 


2 0= -2) 


0.912 


0.928 (0.999, 21.9, 0.211) 


0.949 



TABLE VIII: Test for the Cauchy convergence of the EP-approximants. The values quoted assume optimization on the 
extrinsic parameters but the same intrinsic parameters (i.e., they assume the same masses). Here we define EPo = EP(0,0), 
EPi = EP(1, 1.5), EP2 = EP(2, 2.5), and EP3 = EP(3, 3.5, 8, zi = z% = 0). The values in parentheses are the maxmax matches 
obtained by maximizing with respect to the extrinsic and intrinsic parameters, shown together with the EPjv+i parameters M 
and r\ where the maxima are attained. In all cases the integration of the equations is started at a GW frequency of 20 Hz. 




FIG. 14: Frequency-domain amplitude versus frequency for the EP-approximated (restricted) waveform, at different PN orders, 
for a (15 + 15)M BBH. 



for various values of the (reduced) angular momentum p v . As this plot shows, there exists a critical radius, r cr it, 
below which no circular orbits exist. This r cr i t can be derived as follows. From Fig. ^| (left), we deduce that 
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Combining these equations we obtain two conditions that define r cr i t : 
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FIG. 15: Ending points of the H models at 3PN order for low values of n. In the left panel, we plot as a function of r 
the Hamiltonian H(r,p r = 0,y v ) [given by Eq. I|56jl ] . evaluated at n = 0.16 for a (5+20)A/q BBH, for various values of the 
(reduced) angular momentum p v . The circular-orbit solutions are found at the values of r and H joined by the dashed line. At 
»*crit — 4.524 there is a critical radius, below which there is no circular orbit. In the right panel, we plot as a function of n the 
orbital angular frequency <2> C rit(??) corresponding to the critical radius, for 0.1 < n < 0.21 (solid line). This curve agrees well 
with the ending frequencies of the HT and HP models at 3PN order, which are shown as dotted and dashed lines in the figure. 



In the right panel of Fig.^J we plot the critical orbital frequency a} cr it as a function of r\ in the range [0.1, 0.21]. In the 
same figure, we show also the ending frequencies for the HT(3, 3.5, ±2) and HP(3, 3.5, ±2) models. For 0.1 < r\ < 0.21, 
these ending frequencies are in good agreement with the critical frequencies d> cr it; for r\ > 0.21, the ending condition 
r = 0.3(ra)) is satisfied before w = 0. For 0.1 < n < 0.21, this good agreement can be explained as follows: for the 
H models at 3PN order with rj < 0.21, the orbital evolution is almost quasicircular (i.e., r remains small and Q keeps 
increasing) until the critical point is reached; beyond this point, there is no way to keep the orbit quasicircular, as the 
angular motion is converted significantly into radial motion, and Q begins to decrease. This behavior (D — ► 0) is also 
present in the E model in the vicinity of the light ring, because the light ring is also a minimal radius for circular orbits 
[the conditions (|100fl are satisfied also in this case]. However, the behavior of the energy is qualitatively different for 
the H and E models: in the E models (just as for a test particle in Schwarzchild spacetime) the circular-orbit energy 
goes to infinity, while this is not the case for the H models. 

The second behavior is usually caused by radiation-reaction effects, and accelerated by the presence of an ISCO (and 
therefore of a plunge). However, it is worth to mention another interesting way in which the criterion r — 0.3(ru>) can 
be satisfied for some E evolutions at 3PN order. During the late stages of evolution, r sometimes increases suddenly 
and drastically, and the equations of motion become singular. This behavior is quite different from a plunge due to 
the presence of an ISCO (in that case the equations of motion do not become singular). The cause of this behavior 
is that at 3PN order the coefficient D{r) [see Eq. I|89|l ] can go to zero and become negative for a sufficiently small r. 
For z\ — Z2 — 0, this occurs at the radius rjj given by 



6nr D + 2(3?7 - 26)77 = 0; 



(105) 



ro can fall outside the light ring. For example, for 77 = 0.25 we have td = 2.54, while the light rings sits at r = 2.31. 
On the transition from D{r) > to D(r) < 0, the effective EOB metric unphysical, and the E model then becomes 
invalid. Using the Hamiltonian equation of motion (|91|l . it is straightforward to prove that a negative D(r) causes 
the radial velocity to blow up: 



Pr 



dH 



dp r D(r) 



as 



(106) 



SIGNAL-TO-NOISE RATIO FOR THE TWO-BODY MODELS 



In Fig. E| we plot the optimal signal-to- noise ratio p op t for a few selected PN models. The value of p op t is computed 
using Eqs. and l|19|) with the waveform given by Eq. (|29|l . for a luminosity distance of 100 Mpc and the rms = 8/5 
[see discussion around Eq. Il29fl ]: for the EP model we set z\ = z% = 0. Notice that, because the E models have a 
plunge, their signal-to-noise ratios are much higher (at least for M > 30M Q ) than those for the adiabatic models, 
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FIG. 16: Signal-to-noise ratio at 100 Mpc versus total mass M, for selected PN models. The S/N is computed for equal-mass 
BBHs using the LIGO-I noise curve 1281 and the waveform expression 1291 with the rms Q = 8/5; for the E model at 3PN we 
set Z\ — 5,2 = 0. 




FIG. 17: Effect of the plunge on the signal-to-noise ratio. The S/N is computed at 100 Mpc for equal- mass BBHs, as a function 
of the total mass, for the T(2, 2) adiabatic model (for comparison), and for the EP(2, 2.5) model with ending frequency at the 
ISCO, and at the light ring (in this latter case the signal includes a plunge). Here we use the LIGO-I noise curve 1281 and the 
waveform expression 1291 with the rms Q = 8/5. 



which we cut off at the MECO. See also Fig. El which compares the S/N for EP(2, 2.5) waveforms with and without 
the plunge; for M = 2OM , excluding the plunge decreases the S/N by ~ 4% (which corresponds to a decrease in 
detection rate of 12% for a fixed detection threshold); while for M — 30M Q , excluding the plunge decreases the S/N 
by ~ 22% (which corresponds to a decrease in detection rate of 54%). This result confirms the similar conclusion 
drawn in Ref. 12] . 

Because at 2PN and 3PN order the H models do not have a plunge, but the two BHs continue to move on 
quasicircular orbits even at close separations, the number of total GW cycles is increased, and so is the signal-to-noise 
ratio, as shown in the right panel of Fig. 1161 However, we do not trust the H models much, because they show a very 
different behavior at different PN orders, as already emphasized in Sec. IIV Al 
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VI. PERFORMANCE OF FOURIER-DOMAIN DETECTION TEMPLATES, AND CONSTRUCTION 
OF A FOURIER-DOMAIN DETECTION-TEMPLATE BANK 

In the previous sections we have seen [for instance, in Table IXl] that the overlaps between the various PN waveforms 
are not very high, and that there could be an important loss in event rate if, for the purpose of detection, we restricted 
ourselves to only one of the two-body models [see Figs. 1161 II 7j . To cope with this problem we propose the following 
strategy. We guess that the conjunction of the waveforms from all the PN models spans a region in signal space 
that includes (or almost includes) the true signals, and we build a detection template family that embeds all the PN 
models in a higher-dimensional space. The PN models that we have considered (expanded and resummed, adiabatic 
and nonadiabatic) rely on a wide variety of very different dynamical equations, so the task of consolidating them under 
a single set of generic equations seems arduous. On the other hand, we have reason to suspect, from the values of 
the matches, and from direct investigations, that the frequency-domain amplitude and phasing (the very ingredients 
that enter the determination of the matches) are, qualitatively, rather similar functions for all the PN models. We 
shall therefore create a family of templates that model directly the Fourier transform of the GW signals, by writing 
the amplitude and phasing as simple polynomials in the GW frequency Jew- We shall build these polynomials with 
the specific powers of few that appear in the Fourier transform of PN expanded adiabatic waveforms, as computed 
in the stationary-phase approximation. However, we shall not constrain the coefficients of these powers to have the 
same functional dependence on the physical parameters that they have in that scheme. More specifically, we define 
our generic family of Fourier-domain effective templates as 



where 



hes(f) = Aff(/)e"^\ (107) 
Aff (/) = r 7/6 (l - a / 2 / 3 ) 0(/ cut - /) , (108) 

Veff(/) = 2^/to + 00 + r 5/3 (Vo + V>l/2 / V3 + V>1 f /3 + V>3/2 / + V^2 f 1 * + •••), (109) 

where to and </>o are the time of arrival and the frequency-domain phase offset, and where 0{. . .) is the Heaviside step 
function. This detection template family is similar in some respects to the template banks implicitly used in Fast 
Chirp Transform techniques 63]. However, because we consider BBHs with masses 10-40M Q , the physical GW signal 
can end within the LIGO frequency band; and the predictions for the ending frequency given by different PN models 
can be quite different. Thus, we modify also the Newtonian formula for the amplitude, by introducing the cutoff 
frequency / cut and the shape parameter a. 

The significance of / cu t with respect to true physical signals deserves some discussion. If the best match for the 
physical signal g is the template hf oat , which ends at the instantaneous GW frequency / cut (so that hf cnt (f) ~ g(f) 
for / < / cut and hf cnt (f) — for / > / cu t), then we can be certain to lose a fraction of the optimal p that is given 
approximately by 



r/cut \~g(f)\ 2 , f roo \g(f)\ 2 M 

pent V JO S n (f) aJ „, 1 1 -//cut SZm aJ / 11fV . 



On the other hand, if we try to match g with the same template family without cuts (and if indeed the /i's are 
completely inadequate at modeling the amplitude and phasing of g above / cu t), then even the best-match template 
/inocut (defined by h nocut (f) ~ g(f) for / < / cut , and by zero correlation, h nocut (f)g*(f) ~ for / > / cut ) will yield 
an additional loss in p caused by the fact that we are spreading the power of the template beyond the range where it 
can successfully match g. Mathematically, this loss comes from the different normalization factor for the templates 
/i/ cot and ft-nocut, and it is given by 



r/c„t \HfWjf 1 r °° \hU)f Jf 

Pnocut < V-Jo S n (f) J _ 1 _ 1 J /cut S n (f) 1 ( U1 \ 



Pent t /r°° MfWJ f ' ' 2 f 00 Wf , )l2 df 



If we assume that g and h nocu t have roughly the same amplitude distribution, the two losses are similar. 

In the end, we might be better off cutting templates if we cannot be sure that their amplitude and phasing, beyond 
a certain frequency, are faithful representations of the true signal. Doing so, we approximately halve the worst-case 
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1 r°° \ h (f)\ 2 r1 f 1 r°° ig(/)i r°° ig(/)i % 7f 



loss of p, because instead of losing a factor 

PnOCUt Pcut ^ 1 * " Jcut ^n\J I • ^ " Jcut ^n\J I * ^ ^ " ■) cut "»UJ ' (112) 

we lose only the factor /? C ut/Popt- On the other hand, we do not want to lose the signal-to-noise ratio that is accumulated 
at high frequencies if our templates have a fighting chance of matching the true signal there; so it makes sense to 
include in the detection bank the same template with several different values of / cu t- 

It turns out that using only the two parameters ipo and V'3/2 m the phasing (and setting all other -0 coefficients to 
zero) and the two amplitude parameters, f cu t and a, we obtain a family that can already match all the PN models of 
Secs. lllllllVl with high fitting factors FF. This is possible largely because we restrict our focus to BBHs with relatively 
high masses, where the number of GW cycles in the LIGO range (and thus the total range of the phasing ip(f) that 
we need to consider) is small. 

In Tab. IX Vl we list the minmax (see Sec.|nJ fitting factor for the projection of the PN models onto our frequency- 
domain effective templates, for a set of BBH masses ranging from (5 + 5)M Q to (20 + 20) Mq. In computing the 
fitting factors, we used the simplicial search algorithm amoeba |6l| to search for the optimal set of parameters 
("0o> "03/27 /cut, ot) (as always, the time of arrival and initial phase of the templates were automatically optimized as 
described in Sec.[H]). From Tab. IXVl we draw the following conclusions: 

1. All the adiabatic models (T and P) are matched with fitting factors FF > 0.97. Lower-mass BBHs are matched 
better than higher-mass BBHs, presumably because for the latter the inspiral ends at lower frequencies within 
the LIGO band, producing stronger edge effects, which the effective templates cannot capture fully. 3PN models 
are matched better than 2PN models. 

2. The Effective-One-Body models (ET and EP) are matched even better than the adiabatic models, presumably 
because they have longer inspirals and less severe edge effects at the end of inspiral. Unlike the adiabatic 
models, however, ET and EP are matched better for higher-mass BBHs. In fact, all the FFs are > 0.99 except 
for (5 + 5)Mq BBHs, where FF > 0.979. The reason for this is probably that this low-mass BBH has more GW 
cycles in the LIGO frequency band than any other one, and the two phasing parameters of our effective templates 
cannot quite model the evolution of the phasing. [In the adiabatic models, these effects may be overshadowed 
by the loss in signal to noise ratio due to the edge effects at high frequencies.] When the parameters Zi t 2 are 
allowed to be nonzero, the matches get worse, but not by much. For all the plausible values of z\, the worst 
situation seems to happen at z\ = —40, where the overlaps are still higher than ~ 0.95 [with minimum 0.947.] 

3. The Hamiltonian models (HT and HP) at 3PN order are not matched as precisely, but the detection template 
family still works reasonably well. We usually have FF > 0.96, but there are several exceptions, with FF as 
low as 0.948. For these models, the overlaps are lower in the equal-mass cases, where the ending frequencies of 
the waveforms are much higher than for the other models; it seems that the effective templates are not able to 
reproduce this late portion of the waveforms (this might not be so bad, because it does not seem likely that this 
part of the signal reflects the true behavior of BBH waveforms). 

4. The Lagrangian models (L) are matched a bit worse than the Hamiltonian models (HT and HP) at 3PN, but 
they still have FF higher than 0.95 in most cases, with several exceptions [at either (20 + 20)M Q or (5 + 5)M Q ], 
which can be as low as 0.93. 

5. HT and HP models at 2PN are matched the worst, with typical values lower than 0.95 and higher than 0.85. 

Finally, we note that our amplitude function *4 G ff(/) is a linear combination of two terms, so we can search 
automatically over the correction coefficient a, in essentially the same way as discussed in Sec. [H] for the orbital 
phase. In other words, a is an extrinsic parameter. [Although we do search over a, it is only to show the required 
range, which will be a useful piece of information when one is deciding how to lay down a mesh of discrete templates 
on the continuous detection-template space.] 



A. Internal match and metric 



To understand the matches between the Fourier-domain templates and the PN models, and to prepare to compute 
the number of templates needed to achieve a given (internal) MM, we need to derive an expression for the match 
between two Fourier-domain effective templates. 
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We shall first restrict our consideration to effective templates with the same amplitude function (i.e., the same a 
and /cutoff)- The overlap (h(ipo, ^3/2) > h(ipo + A0O: "03/2 + A0> 3 / 2 )) between templates with close values of ipo and ^3/2 
can be described (to second order in Aip and Aip 3 / 2 ) by the mismatch metric gij psj : 



(h(ip , ip3/2),h(i/j + A^ , V'3/2 + A-03/2)) = 1 - ^2 3ij A^A^j 

»,j=0,3/2 

The metric coefficients can be evaluated analytically from the overlap 

(h(i>o, tp3/2),m>o + A0o, -03/2 + A^ 3/2 )> ~ 



max 

A0 o ,Ai o 



M W)\ 3 

dj - , „, cos 



Sh(f) 



max 

A0 o ,At o 



4f 



\Af)\ 2 
s h (f) 



+ + 2irfAt 



df 



\Af)\ 2 



Sh(f) 
\A(f)\ 2 



df- 



S h (f) 



where no = 5/3 and 713/2 = 2/3. Comparison with Eq. I|113|) then gives 
V 5 ij M>iN>i = \ A mm { ( At/;o A^ 3/2 ) M (1 



A0 O 

A^3/2 



2(A0 O 2^Ato)M (2) ( A A ^ 2 



A0„ 2 7 rAt )M (3) 



A0 O 



where the M(i),..(3) are the matrices 



M (1) = 
M (2) = 



M (3) = 



J(2n ) J(n Q + n 3 / 2 ) 
J(n + n 3/2 ) J(2n 3/2 ) 

J (no) J(n 3/2 ) 
J(n - 1) J(n 3/2 - 1) 

J(0) J(-l)" 
J(-l) J(-2) 



and where 



J(n) = 



4f 



l^(/)l 2 1 
Sh(f) f n 



df 



\A(f)\'< 



Sh(f) 



(113) 



(114) 



(115) 

(116) 
(117) 
(118) 

(119) 



Since M( 3 ) describes the mismatch caused by {A(j)o, Ato), it must be positive definite; because the right-hand side 
of 111511 reaches its minimum with respect to variations of A4>o and Ato when 



2M 



(2) 



A^ 

A0 3 / 2 



2M r 



A0 O 
W I 2vrAt 



= 0, 



we obtain 



to = 2 



M (1) - Mf 2) M^M (2) 



(120) 



(121) 



We note also that the mismatch (hfyo, 03/2), M"0o + A"0Oj 03/2 + A-0 3 / 2 )) is translationally invariant in the (ipo, ^3/2) 
plane, so the metric is constant everywhere. In the left panel of Fig. ^0 we plot the iso-match contours (at 
matches of 0.99, 0.975 and 0.95) in the (At/> , A?/> 3 / 2 ) plane, as given by the metric (|121|l [solid ellipses], compared 
with the actual values obtained from the numerical computation of the matches [dashed lines]. For our purposes, 
the second-order approximation given by the metric is quite acceptable. In this computation we use a Newtonian 
amplitude function A(f) = f~ 7 ^ 6 [i.e., we set a — and we set our cutoff frequency at 400 Hz]. 

We move now to the mismatch induced by different cutoff frequencies /cut- Unlike the case of the Vo, 03/2 
parameters, this mismatch is first order in A/ cut , so it cannot be described by a metric. Suppose that we have two 



35 




FIG. 18: In the left panel, we plot the iso-match contours for the function {h(ipo, V'3/2), h(%j)Q + Aipo, V>3/2 + Aip 3 / 2 )); contours 
are given at matches of 0.99, 0.975 and 0.95. Solid lines give the indications of the mismatch metric; dashed lines give actual 
values. Here we use a Newtonian amplitude function A(f) = / ~ 7 ^ 6 [we set a = and we do not cut the template in the 
frequency domain. In fact / cut = 400 Hz]. In the right panel we plot the values of A/ cut (versus /cut) required to obtain 
matches (h(f cut ), h(f cut + A/ Cu t)) of 0.95 (uppermost curve), 0.975 and 0.99 (lowermost). In the region below each contour the 
match is larger than the value quoted for the contour. Again, here we use a Newtonian amplitude function A(f) = f~ 7 ^ 6 [we 
set a = 0]. 



effective templates h(f cut ) and h(f cut + A/ cut ) with the same phasing and amplitude Af > 0, but different cutoff 
frequencies. The match is then given by 



(M/cut),M/cut+A/ cut )) = 



\Mf)\ 
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1/2 



Jo a J s h u) J |yo 
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This result depends strongly on / cut . In the right panel of Fig. ^] we plot the values of A/ cut that correspond to 
matches of 0.95, 0.975 and 0.99, according to the first order approximation [solid lines], and to the exact numerical 
calculations [dashed lines] , both of which are given in the second line of Eq. i|123|) ■ In the region below each contour the 
match is larger than the value that characterizes the contour. As we can see from the graph, the linear approximation 
is not very accurate, thus in the following we shall use the exact formula. 



B. Construction of the effective template bank: parameter range 

All the PN target models are parametrized by two independent numbers (e.g., the two masses or the total mass 
and the mass ratio); if we select a range of interest for these parameters, the resulting set of PN signals can be seen 
as a two-dimensional region in the (mi, 7712) or (M, 77) plane. Under the mapping that takes each PN signal into the 
Fourier-domain effective template that matches it best, this two-dimensional region is projected into a two-dimensional 
surface in the ("0o> ^3/2' /cut) parameter space (with the fourth parameter a — 0). As an example, we show in Fig. 1191 
the projection of the ET(2, 2.5) waveforms with (single-BH) masses 5-20 M Q . The 26 models tested in Sees. IHll II VI 
would be projected into 26 similar surfaces. In constructing the detection template families, we shall first focus on 
17 of the 26 models, namely the adiabatic T and P models at 2PN and 3PN, the E models at 2PN and at 3PN but 
with zx 2 = 0, and the H models at 3PN. We will comment on the E models with 7^ 0, on the L models, and on 
the HT and HP models at 2PN order at the end of this section. 

It is hard to visualize all three parameters at once, so we shall start with the phasing parameters ipo an d V'3/2- 
In Fig. 1201 we plot the (V'0,'03/2) section of the PN-model projections into the (ipo, 1^3/2, /cutoff) space, with solid 
diamonds showing the projected points corresponding to BBHs with the same set of ten mass pairs as in Tab. IXVI 
Each PN model is projected to a curved-triangular region, with boundaries given by the sequences of BBHs with 
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FIG. 19: Projection of the ET(2, 2.5) waveforms onto the frequency-domain effective template space. For a we choose the 
optimal value found by the search. The (ipo, V'3/2, /cut) surface is interpolated from the then mass pairs shown in Tab. |W| 
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TABLE IX: End-to-end matches and ending frequencies along the BH mass lines of Fig. 1201 The first three columns show 
the end-to-end matches and the corresponding number of templates (for MM ~ 0.98) along the BH mass lines; the remaining 
columns show the minimum ending frequencies of PN waveforms along the BH mass lines, the match between the two effective 
templates at the ends of the range, and the number of templates needed to step along the range while always maintaining 
a match ~ 0.98 between neighboring templates. When computing these matches, we use a Newtonian amplitude function 
-4(/) = /" 7/6 I we set a = 0], and we maximize over the parameters ipo and 1P3/2 (which is equivalent to assuming perfect 
phasing synchronization) . 

masses (m + to) (equal mass), (20 + to) and (to + 5). In Fig.[3U| these boundaries are plotted using thin dashed lines, 
for the models T(2,2.5) (the uppermost in the plot), HT(3, 3.5,0 = 2) (in the middle), and P(2,2.5) (lowest). 

As we can see, different PN models can occupy regions with very different areas, and thus require a very different 
number of effective templates to match them with a given MMt- Among these three models, T(2,2.5) requires the 
least number of templates, P(2,2.5) requires a few times more, and HT(3, 3.5,0 = 2) requires many more. This is 
consistent with the result by Porter 60] who found that, for the same range of physical parameters, T waveforms are 
more closely spaced than P waveforms, so fewer are needed to achieve a certain MM. In this plot we have also linked 
the points that correspond to the same BBH parameters in different PN models. In Fig. E01 these lines (we shall call 
them BH mass lines) lie all roughly along one direction. 

A simple way to characterize the difference between the PN target models is to evaluate the maxmax end-to-end 
match between effective templates at the two ends of the BH mass lines (i.e., the match between the effective templates 
with the largest and smallest "03/2 among the projections of PN waveforms with the same mass parameters toi, 7712)5 
we wish to focus first on the effects of the phasing parameters, so we do not cut the templates in the frequency domain 
and we set a = 0. We compute also a naive end-to-end number of templates, A/" en d to end , by counting the templates 
required to step all along the BH mass line while maintaining at each step a match ~ 0.98 between neighboring 
templates. A simple computation yields A/" cn dtocnd = log(end-to-end match)/ log(0. 98). The results of this procedure 
are listed in Table ITXl Notice that, as opposed to the fitting factors between template families computed elsewhere in 
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FIG. 20: Projection of the PN waveforms onto the (^0,^3/2) plane, for BBHs with masses (5 + 5)Mq, (10 + 5)Mq, . . . , 
(20 + 20)M Q (see Tab. lXVt . The projection was computed by maximizing the maxmax match over the parameters ipo, ^3/2 
and / cut ; the correction coefficient a was set to zero. The thin dotted and dashed lines show the boundaries of the projected 
images for the models (from the top) T(2, 2.5), HT(3, 3.5, 9 — 2) and P(2, 2.5). Solid lines (the BH mass lines) link the images 
of the same BBH for different PN models. The ends of the BH mass lines are marked with the BBH masses and with the 
minimum value min{/ en d, /cut} across all the PN models. The thick dashed lines delimit the region that will be covered by 
the effective template bank; the (^ojV'3/2) coordinates are marked on the vertices. The region is further subdivided into four 
subregions I-IV that group the BH mass lines with very similar ending frequencies /cndmm- 



this paper (which arc maximized over the BBH mass parameters of one of the families), these matches give a measure 
of the dissimilarity between different PN models for the same values of the BBH parameters; thus, they provide a 
crude estimate of how much the effective template bank must be enlarged to embed all the various PN models. 

We expect that the projection of a true BBH waveform onto the (^0,^3/2) plane will lie near the BH mass line 
with the true BBH parameters, or perhaps near the extension of the BH mass line in either direction. For this reason 
we shall lay down our effective templates in the region traced out by the thick dashed lines in Fig. 1201 which was 
determined by extending the BH mass lines in both directions by half of their length. 

We move on to specifying the required range of / cut for each (-00, ^3/2) ■ F° r a given PN model and BBH mass 
parameters, we have defined the ending frequency / on( j as the instantaneous GW frequency at which we stop the 
integration of the PN orbital equations. We find that usually the / cu t of the optimally-matched projection of a PN 
template is larger than the / on( j of the PN template. This is because the abrupt termination of the PN waveforms 
in the time domain creates a tail in the spectrum for frequencies higher than / en d- With / cut > / on d and a > 0, the 
effective templates can mimic this tail and gain a higher match with the PN models. In some cases, however, the 
optimal / cu t can be smaller than / Gn d [for example, P(2, 2.5) with (10 + 5)M©, (15 + 5)M Q and (10+10)Af©] suggesting 
that the match of the phasing in the entire frequency band up to / en d is not very good and we have to shorten the 
Fourier-domain template. Now, since we do not know the details of the plunge for true BBH inspiral, it is hard to 
estimate where the optimal / cu t might lie, except perhaps imposing that it should be larger than nnn(/ cnc ] , / cu t ) ■ A 
possibility is to set the range of / cu t to be above /cutmin = min{/ C ut) /end}> with the minimum evaluated among all 
the PN models. 
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In Table IIXI we show the / cu t m in found across the PN models for given BBH mass parameters. We have also 
marked this minimum frequency in Fig. 1201 under the corresponding BH mass lines. In the table we also show the 
match of the two detection templates h(f cut = /cut min) and h(f cut = +00), and the number N^ ss \ inc of intermediate 
templates with different / cu t needed to move from /i(/ C utmin) to h(+oo) while maintaining at each step a match ~ 0.98 
between neighboring templates. It is easy to see that this number is A/^ssline = l°g(M/cut min), h(+oo)) / log(0.98). 
The match was computing using a Newtonian amplitude function A(f) = f~ 7 ^ 6 [we set a = 0], and maximizing 
over the parameters "00 and ^3/2- Under our previous hypothesis that the projection of a true BBH waveform would 
lie near the corresponding BH mass line, we can use the numbers in Table IIXI to provide a rough estimate of the 
range of / cut that should be taken at each point (ipo, ^3/2) within the dashed contour of Fig. |3U] We trace out four 
subregions I, II, III, IV, such that the BH mass lines of each subregion have approximately the same values of / cutmin ; 
we then use these minimum ending frequencies to set a lower limit for the values of / cut required in each subregion: 

/cutmin(I) — 143, /cutmin(H) — 192, /cutmin 

(III) = 232, /cutmin(IV) = 346. The maximum / cu t is effectively set by 
the detector noise curve, which limits the highest frequency at which signal-to-noise can be still accumulated. 

Moving on to the last parameter, a, we note that it is probably only meaningful to have a/^/f < 1, so that -4 e ff(/) 
cannot become negative for / < / cut . [A negative amplitude in the detection template will usually give a negative 
contribution to the overlap, unless the phasing mismatch is larger than tt/2, which does not seem plausible in our 
cases.] Indeed, the optimized values found for a in Tab. IX Vl seem to follow this rule, except for a few slight violations 
that are probably due to numerical error (since we had performed a search to find the optimal value of a). For the 
17 models considered here, the optimal a is always positive [Tab. IXV) which means that, due to cutoff effects, the 
amplitude at high frequencies becomes always lower than the / _7 ^ 6 power law. So for the 17 models considered in 
this section < af^it < 1- [Note that this range will have to be extended to include negative a's if we want to 
incorporate the models discussed in Sec. I VI El ] 



C. Construction of the effective templates bank: parameter density 

At this stage, we have completed the specification of the region in the (ipo, V , 3/2i /cut, ct) parameter space where we 
shall lay down our bank of templates. We expect that the FF for the projection of the true physical signals (emitted 
by nonspinning BBHs with M — 10-A0Mq) onto this template bank should be very good. We now wish to evaluate 
the total number of templates N needed to achieve a certain MM. 

We shall find it convenient to separate the mismatch due to the phasing from the mismatch due to the frequency cuts 
by introducing two minimum match parameters MM^ and MM cu t, with MM = MM,/, • MM cut MM,/, + MM cut - 1. 
As mentioned at the beginning of this section, the correction coefficient a is essentially an extrinsic parameter [see 
Sec. Ill B) : we do not need to discretize the template bank with respect to a, and there is no corresponding MM 
parameter. 

We evaluate M in three refinement steps: 

1. We start by considering only the phasing parameters, and we compute the parameter area Si [in the (V'OiV'3/2) 
plane] for each of the subregions i = I, II, III, IV of Fig. [2U1 We then multiply by the determinant yfg of the 
constant metric, and divide by 2(1— MM,/,), according to Eq. 1)25(1 . to get 

M = y S iy /g 

^2(1 -MM,/,) v ; 

This expression is for the moment only formal, because we cannot compute y/g without considering the amplitude 
parameters a and / cu t- 

2. Next, we include the effect of / cu t- In the previous section, we have set / m incut for each of the subregions by 
considering the range swept by / en d along the mass lines. Recalling our discussion of A/^ssline' we approximate 
the number of distinct values of / cu t that we need to include for each parameter pair (1^0,1^3/2) as 

cut/, , \ -1 fog( / «(V'0,'/'3/2,a,/mincut),/i(V'0,l/'3/2,a,nOCUt)) 

rii (i>o,fa/2,a)- 1 + i — ™ • ( 125 ) 

log MMcut 

2 /3 

For a in the physical range < a < / cut this match is minimized for a = 0, so this is the value that we use 
to evaluate the n£ ut 's. Note that the choice of cutoff frequencies does not depend on the values of the phasing 
parameters. This allows us to have a single set of cutoff frequencies for all points in one subregion. For subregion 
i, we denote this set by Fj. 
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3. The final step is to include the effect of a and / cu t on the computation of ^fg. For simplicity, we shoot for an 
upper limit by maximizing ^fg with respect to a. [Because a is essentially an extrinsic parameter, we do not 
multiply Af by the number of its discrete values: the matches are automatically maximized on the continuous 
range < a < f c ^J 3 ■] Our final estimate for the total number of templates is 

^ = 2(l-U) £Si E " aX ^ (126) 

v VJ i /outers 



We have evaluated this Af numerically. We find that the contributions to the total number of templates from the 
four subregions, for MM = 0.96 (taking MM + = MM cut = 0.98 ), are W(I) ~ 6, 410, Af{U) ~ 2, 170, M{lll) ~ 1, 380, 
Af(TV) ~ 1,230, for a total of N = 11,190. This number scales approximately as [0.04/(1 - MM)] 2 . Notice that 
subregion I, which contains all the BBHs with total mass above 25M Q , requires by far the largest number of templates. 
This is mostly because these waveforms end in the LIGO band, and many values of / cu t are needed to match different 
ending frequencies. 

Remember that the optimal signal-to-noise ratio p for filtering the true GW signals by a template bank is approxi- 
mately degraded (in the worst case) by the factor MM T = FF + MM - 1 [7^. 

While MM depends on the geometry of the template bank, we can only guess at the fitting factor FF for the 
projection of the true signal onto the template space. In this section we have seen that all PN models can be projected 
onto the effective frequency-domain templates with a good FF: for a vast majority of the waveforms FF > 0.96 (and 
the few exceptions can be explained). It is therefore reasonable to hope that the FF for the true GW signals is ~ 0.96, 
so the total degradation from the optimal p will be MMt ^ 0.92, corresponding to a loss of % 22% in event rate. This 
number can be improved by scaling up the number of templates, but of course the actual FF represents an upper 
limit for MMt- For instance, about 47,600 templates should get us MMt ^ 0.94, corresponding to a loss of < 17% 
in event rate. 



D. Parameter estimation with the detection template family 

Although our family of effective templates was built for the main purpose of detecting BBHs, we can still use 
it (once a detection is made) to extract partial information about the BH masses. It is obvious from Fig. [2U| that 
the masses cannot in general be determined unambiguously from the best- match parameters [i.e., the projection of 
the true waveform onto the (tpo,ipa/2) plane], because the images of different PN models in the plane have overlaps. 
Therefore different PN models will have different ideas, as it were, about the true masses. Another way of saying this 
is that the BH mass lines can cross. 

However, it still seems possible to extract at least one mass parameter, the chirp mass A4 = Mr] 3 / 5 , with some 
accuracy. Since the phasing is dominated by the term i/iof -5 ^ 3 at low frequencies, we can use the leading Newtonian 
term ip^if) — jfg (irMf)~~ 5 / 3 obtained for a PN expanded adiabatic model in the stationary-phase approximation to 
infer 

o / 1 \V3 1 / Q \3/5 

~ _ M appiox = - — — . (127) 

If this correspondence was exact, the BH mass lines in Fig. 1201 would all be vertical. They are not, so this estimation 
has an error that gets larger for smaller ip (i.e., for binaries with higher masses). In Table Ixl we show the range of 
chirp-mass estimates obtained from Eq. (|127|) for the values of ipo at the projections of the PN models in Fig. 1201 
together with their percentage error e = (A4^££°* — A4^^ ox )/M. In this table, .M max and A4 m i n correspond to the 
endpoints of the BH mass lines. If we take into account the extension of the BH mass lines by a factor of two in the 
effective template bank, we should double the e of the table. 

It seems quite possible that a more detailed investigations of the geometry of the projections into the effective 
template space (and especially of the BH mass lines) could produce better algorithms to estimate binary parameters. 
But again, probably only one parameter can be estimated with certain accuracy. 



E. Extension of the two-dimensional Fourier-domain detection template 



In our construction of the effective template bank, we have been focusing until now on a subset of 17 models. The 
models we left out are: E models at 3PN with z\$ nonzero, HT and HP models at 2PN, and L models. 
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(5 + 5)M 


4.35 


4.16 


4.27 


2.6 


(10 + 5)M 


6.08 


5.75 


6.00 


4.2 


(15 + 5)M 


7.33 


6.85 


7.28 


5.9 


(10 + 1O)M 


8.71 


8.10 


8.72 


7.1 


(20 + 5)M 


8.33 


7.55 


8.31 


9.1 


(15 + 1UJA2 


10.62 


9.76 


10.96 


11.3 


(20 + 1O)M 


12.17 


10.92 


12.50 


13.0 


(15 + 15)M 


13.06 


11.69 


14.88 


24.4 


(20 + 15)M 


15.05 


13.15 


17.74 


30.6 


(20 + 2O)M 


17.41 


14.91 


21.52 


38.0 



TABLE X: Estimation of the chirp masses M from the projections of the PN target models onto the Fourier-domain effective 
template space. The numbers in the second column (labeled ".M") give the values of the chirp mass corresponding to the 
BH masses to their left; the numbers in the third and fourth columns give the range of estimates obtained from Eq. 1127H 
for the values of ipo at the projections of the target models shown in Fig. I2UI The last column shows the percentage error 
e={M^ii°*-M^)/M. ' 
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FIG. 21: Projection of the E models with nonzero z\ into the (ipo, ^3/2) plane (shown in black dots.) The new points sit quite 
well along the BH mass lines of the 17 models investigated in Secs. lVIBllVI Cl and lVI Dl We use the notation EP(3, 3.5, 9, Zi, Z2) 
and denote by EP(T3, ...) the two-body model in which the coefficient A(r) is PN expanded [see Eq. 1881 1. 



As we can see from Fig. E models with z\ t 2 nonzero have a very similar behavior to the 17 models investigated 
above. Indeed: (i) the projection of the PN waveforms from the same model occupy regions that are triangular, and 
(ii) the projections of PN waveforms of a given mass lies on the BH mass line spanned by the previous 17 models. In 
addition, their projections lie roughly in the region we have already defined in Sees. IVTBl IVI (Jl and IVTDl However, 
the ending frequencies of these models can be much lower than the values we have set for the detection templates: 
the detection templates (in all four subregions) should be extended to lower cutoff frequencies if we decide to match 
these models, up to FF ~ 0.95. A rough estimate shows that this increases the number of templates to about twice 
the original value. 

In Fig. we plot the projections of the L(2,0), L(2, 1), HT(2,2) and HP(2,2.5) waveforms into the (^0,^3/2) 
plane. As we already know, these models are not matched by the detection templates as well as the other 17 models. 
Here we can see that their projections onto the ("0Oj -! / , 3/2) plane are also quite dissimilar from those models. For L 
models, although different masses project into a triangular region, the projection of each mass configuration does 
not align along the corresponding BH mass line generated by the 17 models. In order to cover the L models up 
to FF~ 0.93, we need to expand the (ipo, 1^3/2) region only slightly. However, as we read from Tab. IX Vl the cutoff 
frequencies need to be extended to even lower values than for the E models with nonzero iii,2. Luckily, this expansion 
will not cost much. In the end the total number of templates needed should be about three times the original value. 

For HT and HP models at 2PN, the projections almost lie along the BH mass lines, but the regions occupied by 
these projections have weird shapes. We have to extend the (i/^o, ^3/2) region by a factor ~ 2 in order to cover the 
phasings. [The ending/cutoff frequencies for these models are higher than for the previous two types of models.] An 
additional subtlety in this case is that, as we can read from Tab. IXVI the optimal values of a are often negative, 
since the amplitude becomes higher than the /~ 7 ^ 6 power law at higher frequencies. This expansion of the range of 
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FIG. 22: Projections of HT and HP models at 2PN and L models into the (0o, 03/2) plane (shown in black dots.) The 
projections of the previous 17 models are shown in gray dots. 



a affects both the choice of the discrete cutoff frequencies and the placement of (-00, "03/2) lattices. This effect is yet 
to be estimated. 

Finally, we notice that if these extensions are made, then the estimation of the chirp mass from the coefficient 0o 
becomes less accurate than the one given in Table Ixl 

F. Extension of the Fourier-domain detection template family to more than two phasing parameters 




50 



FIG. 23: Projection of the models P(2,2.5), ET(2, 2.5), ET(3, 3.5,0), and SPA(1.5) onto the three-parameter Fourier-domain 
detection template, for many BBH masses that lie within the same ranges taken in Fig. 1201 The variables (X, Y, Z) are related 
to (0o, 01, 03/2) by a linear transformation, constructed so that the mismatch metric is just Sij and that the (00,0,03/2) plane 
is mapped to the (X, Y, 0) plane. The dots show the value of the parameters (X, Y, Z) where the match with one of the PN 
waveforms is maximum. 

It might seem an accident that by using only two phasing parameters, 0o and 03/2, we are able to match very 
precisely the wide variety of PN waveforms that we have considered. Indeed, since the waveforms predicted by each 
PN model span a two-dimensional manifold (generated by varying the two BH masses m\ and mi or equivalently 
the mass parameters M and ?y), we could naturally expect that a third parameter is required to incorporate all the 
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FIG. 24: (X,Z) section of Fig. 1231 Comparison with Fig. 1231 shows that all the projections lie near the (X,0,Z) plane. 



PN models in a more general family, and to add even more signal shapes that extrapolate beyond the phasings and 
amplitudes seen in the PN models. 

In particular, because the accumulation of signal-to-noise ratio is more sensitive to how well we can match the 
phasing (rather than the amplitude) of PN templates, such a third parameter should probably interpolate between 
phasings predicted by different PN models. As a consequence, the amplitude parameters / cut and A do not generate 
a real dimensional extension of our detection template family. In this section, we present a qualitative study of the 
extension of our detection template family obtained by adding one phasing parameter, the parameter ipi of Eq. (|109fl . 

We use the (^Oj ipi, 1P3/2) Fourier-domain detection templates to match the PN waveforms from the models P(2, 2.5), 
ET(2, 2.5), and ET(3, 3.5,0); these models were chosen because their projections onto the (V'Oj V'3/2) detection tem- 
plates were rather distant in the (V'o 1^3/2) parameter space. Throughout this section (and unlike the rest of this 
paper) , we use an approximated search procedure whereby we essentially replace the amplitude of the target models 
with the Newtonian amplitude A(f) = f~ 7 ^ 6 with a cutoff frequency / cut [we always assumed A = and / cut = 400 
Hz]. As expected, the matches increase, and indeed they are almost perfect: always higher than 0.994 (it should be 
remembered however that these should be considered as matches of the PN phasings rather than as matches of the 
PN waveforms; especially for high masses, the frequency dependence of the amplitude is likely to change these values). 

If we plot the projections of the PN waveforms in the (ipo, ^3/2) space, we find that the clusters of points corre- 
sponding to each PN target model look quite different from the projections [onto the ('0cbV'3/2) template space] 
shown in Fig. I2UI but this is just an artifact of the parametrization. We can perform a linear transformation 
(ipo, ipi, "03/2) — > (A, Y, Z), defined in such a way that (i) in the (A, Y, Z) parameters, the mismatch metric is just Sij, 
and that (ii) the ("0o> 0, 1P3/2) plane is mapped to the (A, Y, 0) plane. These conditions define the linear transformation 
up to a translation and a rotation along the Z axis; to specify the transformation completely we require also that all 
the projections of the PN models lie near the origin, and be concentrated around the A axis. Figure ESI shows the 
projection of the PN models P(2, 2.5), ET(2, 2.5), and ET(3, 3.5, 0) onto the (V>o, ^1^3/2) detection template family, 
as parametrized by the (A, Y, Z) coordinate system, for many BBH masses that lie within the same ranges of Fig. |2"U1 
Each dot marks the parameters (A, Y, Z) that best match the phasing of one of the PN waveforms. We include also 
the projection of a further PN model, SPA(1.5), obtained by solving the frequency-domain version of the balance 
equation, obtained in the stationary-phase approximation from our T model. The expression of the SPA(1.5) phasing 
as a function of / coincides with our Eq. I|109[) . but the coefficients that correspond to (ipo, tpi, "03/2) o^e functions of 
the two mass parameters M and 7/. 

By construction, the match between nearby detection templates is related to their Euclidian distance in the (A, Y, Z) 
by 

1 - overlap = AA 2 + AY 2 + AZ 2 . (128) 

We see immediately that all the PN models are not very distant from the (A, Y, 0) plane [also shown in the figure], 
which coincides with the (ipo, V'3/2) plane. The farthest model is P(2, 2.5), with a maximum distance ~ 0.18. It 
is important to notice that, since this number is obtained by assuming / cut = 400 Hz and A = 0, it tends to 
underestimate the true overlaps for models that end below 400 Hz, such as the P models at higher masses. See also 
Fig.Elfor an (A, Z) section of Fig. EU 

We can study the relation between this three-dimensional family of templates and the two-dimensional family 
considered earlier by projecting the points of Fig. EHonto the (A, Y, 0) plane [which corresponds to the (tpo, 0, V'3/2) 
plane]. The resulting images resemble closely the projections of the PN models onto the (ipo, V'3/2) parameter space of 



43 




'0.0 



0.5 



1.0 1.5 

¥o/10 s 



2.0 




0.0 



0.5 



1.0 1.5 

Vo/10 s 



2.0 



FIG. 25: In this figure, we compare the projection of the PN models onto the three-dimensional (ipo, ipi, ^3/2) Fourier- domain 
detection template family [shown by the dots as a two-dimensional section in the (i/to, ^3/2) submanifold] with the projection 
of the PN models in the two-dimensional (t/'o, ^3/2) template family [shown by the lines]. In the left panel, we use .4 = and 
/cut ~ 400 Hz to maximize the matches; in the right panel we use A = and / cu t = 200 Hz. 



the two-dimensional family, as seen in the left panel of Fig. [23 However, the agreement is poor for P(2, 2.5) because 
of the relatively high cut frequency / cu t = 400 Hz. The right panel of Fig. 1251 was obtained by taking / cut = 200 Hz. 
The agreement is much better. This result goes some way toward explaining why using only two phasing parameters 
was enough to match most PN models in a satisfactory way. 

As stated at the beginning of this section, the parameter Z can indeed be used to expand the dimensionality of 
our detection template family, because it appears to interpolate between different PN models. It is possible that the 
number of Z values needed when laying down a discrete template family might not be too large, because the PN 
models do not seem to lie very far from the Z = plane [remember that distances in the (A, Y, Z) parameter space 
are approximately mismatch distances]. 

The good performance that we find for the two- and three-dimensional Fourier-domain families confirms the results 
obtained in Refs. , Q and ■ In Ref- , the authors point out that the waveforms obtained from the stationary 
phase approximation at 2PN and 2.5PN order are able to approximate the E models, throughout most of the LIGO 
band, by maximizing over the mass parameters [see Ref. [12j, in particular the discussion of their model "Tf2," and 
the discussion around their Fig. 2] . 

In Ref. 0], Chronopolous and Apostolatos show that what would be in our notation the SPA(2) model (where the 
phasing is described by a fourth-order polynomial in the variable Z 1 ^ 3 ) can be approximated very well, at least for the 
purpose of signal detection, by the SPA(1.5) model, with the advantage of having a much lower number of templates. 
In Ref. |§4|, the authors go even further, investigating the possibility of approximating the SPA(2) phasing with a 
polynomial of third, second and even first degree obtained using Chebyshev approximants. 

It is important to underline that in all of these analyses the coefficients that appear in the expression of the phasing 
[corresponding to our ipo,ipi, ... in Eq. i|109l) ] depend on only two BBH mass parameters, either directly 0, |44| . 
or indirectly [64J through specific PN relations at each PN order. As a consequence, the phasings assumed in these 
references are confined to a two-dimensional submanifold analog to the surface labeled "SPA(1.5)" in Fig. 1231 

In this paper we follow a more general approach, because the phasing coefficients i/'i ar e initially left completely 
arbitrary. Only after studying systematically the projection of the PN models onto the template bank we have 
determined the region where a possible detection template bank would be laid down. The high matches that we find 
between detection templates and the various PN models depend crucially on this assumption. As a consequence, our 
parameters ipi do not have a direct physical meaning, and they cannot easily be traced back to specific functions of 
the BBH masses, except for the chirp mass, as seen in Sec. I VI Dl This is natural, because our detection templates are 
built to interpolate between different PN models, each of which has, as it were, a different idea of what the waveform 
for a BBH of given masses should be. 



VII. PERFORMANCE OF THE TIME-DOMAIN DETECTION TEMPLATES AND CONSTRUCTION 

OF THE DETECTION BANK IN TIME DOMAIN 

Another possibility of building a detection template family is to adopt one or more of the physical models discussed 
in Sees. II VI as the effective template bank used for detection. Under the general hypothesis that underlies this work 
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(that is, that the target models span the region in signal space where the true physical signals reside), if we find that 
one of the target models matches all the others very well, we can use it as the effective model; and we can estimate 
its effectualness in matching the true physical signal from its cffectualness in matching all the other models. 

As shown in Tables fxTl IXIII and discussed in Sec.[V] the fitting factors FF for the projection of the PN models onto 
each other are low (at least for PN order n < 2.5 or for high masses); in other words, the models appear to be quite 
distant in signal space. This conclusion is overturned, however, if we let the dimensionless mass ratio rj move beyond 
its physical range < r\ < 1/4. For instance, the P(2, 2.5) and EP(3, 3.5,0) models can be extended formally to the 
range < r\ < 1. Beyond those ranges, either the equations (of energy-balance, or motion) become singular, or the 
determination of the MECO or light ring (the evolutionary endpoint of the inspiral for the P(2, 2.5) model and the 
EP(3, 3.5,0) model, respectively) fails. 

When the models are extended to < rj < 1, they appear to lie much closer to each other in signal space. In 
particular, the P(2, 2.5) and EP(3, 3.5, 0) models are able to match all the other models, with minmax FF > 0.95, for 
almost all the masses in our range, and in any case with much improved FF for most masses; see Tables IXVII and 
IX VIII Apparently, part of the effect of the different resummation and approximation schemes is just to modulate 
the strength of the PN effects in a way that can be simulated by changing 77 to nonphysical values in any one model. 
This fact can be appreciated by looking at Figs. [23 and [2H1 I2H1 which show the projection of several models 
onto the P(2, 2.5) and EP(3, 3.5, 0) effective template spaces, respectively. For instance, in comparison with T(2, 2.5), 
the model P(2,2.5) seems to underestimate systematically the effect of 77, so a satisfactory FF for t\t — 0.25 can be 
obtained only if we let r\p > 0.25 (quite consistently, in the comparison of Tables IXll IXIII where 77 was confined to its 
physical range, T(2, 2.5) could match P(2, 2.5) effectively, but the reverse was not true). 

The other (and perhaps crucial) effect of raising 77 is to change the location of the MECO for the P-approximant 
model (or the light ring, for the EP model), where orbital evolution ends. (Remember that one of the differences 
between the Pade and the EOB models is that the latter includes a plunge part between the ISCO and the light ring.) 
More specifically, for P(2, 2.5) [EP(3, 3.5, 0)] the position of the MECO [light ring] is pushed to smaller radii as 77 is 
increased. This effect can increase the FF for target models that have very different ending frequencies from those of 
P(2, 2.5) and EP(3, 3.5) at comparable 77's. 

Because for the EP model the frequency at the light ring is already quite high, we cannot simply operate on 77 
to improve the match between the EP model and other models that end at much lower frequencies [see the values 
of minmax matches in Table |XV11| . Thus, we shall enhance the effectualness of EP by adding an arbitrary cut 
parameter that modifies the radius r (usually the light-ring radius) at which we stop the integration of the Hamilton 
equations (|91|l - l|94[l : the effect is to modify the final instantaneous GW frequency of the waveform. This is therefore 
a time-domain cut, as opposed to the frequency-domain cuts of the frequency-domain effective templates examined 
in the previous section. 

We can then compute the FF by searching over / cut in addition to M and 77, and we shall correspondingly account 
for the required number of distinct / cut when we estimate the number of templates required to give a certain MM to t ■ 
Even so, if we are unsure whether we can model successfully a given source over a certain range of frequencies that 
falls within LIGO range (as it is the case for the heavy BBHs with MECOs at frequencies < 200 Hz), the correct 
way to estimate the optimal p (and therefore the expected detection rate) is to include only the signal power in the 
frequency range that we know well. 

The best matches shown in Tables IXVII and IXVIII and in Figs. I26H29I were obtained by searching over the target 
model parameter space with the simplicial amoeba algorithm [61|. We found (empirically) that it was expedient to 
conduct the searches on the parameters (3 = Mr/ 2 / 5 and r\ rather than on M and 77. This is because iso-match surfaces 
tend to look like thin ellipses clustered around the best match parameter pair, with principal axes along the (3 and 77 
directions. As shown in Table IX VII the values of the maxmax and minmax FFs are very close to each other for the 
P(2,2.5) model; the same is true for the EP(3,3.5) model (so in Table IXVTTl we do not show both). For EP(3,3.5), 
the search over the three parameters (/3, 77, f r .nt) was performed as a refinement step after a first search on (/3, 77). 

We have evaluated the mismatch metric [28j gij (see Sec. [HJ with respect to the parameters (/3, rf) for the models 
P(2, 2.5) and EP(3, 3.5, 0) (while evaluating </y, the EP waveforms were not cut). The metric components at the point 
(/3o,%) were obtained by first determining the ranges (/3 min , /3 m ax), (Vn, w) for which 

(u(f3o,Vo),u(Pmin,Vo)) = («(A)> »7o), u(/3 max , rjo)) = 1 - 0.05 (129) 
(u(f3 , ryo), u(fl , f] min )) = (u(J3 , 770), u(J3 , r7 max )) = 1 - 0.05; (130) 

then a quadratic form was least-squares-fit to 16 values of the match along the ellipse Ti with axes given by (/3 m i n , /3 m ax) 
and (77 m i n , 77 max ). The first quadratic form was used only to determine the principal axes of two further ellipses T2 
and Ts, at projected matches of 1 — 0.025 and 1 — 0.0125. Another quadratic form (giving the final result for the 
metric) was then fit at the same time to 16 points along T2 and to 16 points along T3, but the two ellipses were given 
different fitting weights to cancel the quartic correction terms in the Taylor expansion of the match around (/3o,?7o) 
[the cubic terms were canceled automatically by taking symmetric points along the ellipses] . The rms error of the fit 
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FIG. 26: Projection of 2PN waveforms onto the P(2, 2.5) effective template space. Dots are shown for the same BBH masses 
of Tab. ESI and f or PN models T(2, 2.5), P(2, 2.5), ET(2, 2.5), and EP(2, 2.5). The thin solid lines show the BH mass lines 
(introduced in Sec. lVIBjl . while the dashed and dotted lines show the contours of the projections of selected PN models. 
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FIG. 27: Projection of 3PN waveforms onto the P(2, 2.5) effective template space. Dots are shown for the same BBH masses of 
Tab.|XVl and for PN models T(3,3.5,+2), P(3,3.5,+2), ET(3, 3.5, +2), EP(3, 3.5, +2), HT(3, 3.5, +2), and HP(3, 3.5,0). The 
dots for — —2 are only slightly displaced, and they are not shown. The thin solid lines show the BH mass lines (introduced 
in Sec. I VI 51 . while the dashed and dotted lines show the contours of the projections of selected PN models. 



was in all cases very good, establishing that the quadratic approximation held in the close vicinity (matches ~ 0.95) 
of each point . 

We estimate that the numerical error ~ 20% is in any case less than the error associated with using Eq. I|25|) to 
evaluate the required number of templates, instead of laying down a lattice of templates more accurately. 

The resulting \/\g\ for P(2,2.5) and EP(3, 3.5,0) is shown in Fig. 1201 It is evident that most of the mismatch 
volume is concentrated near the smallest /3's and 77's in parameter space. This is encouraging, because it means that 
the extension of the effective template family to high masses and high 7?'s (necessary, as we have seen, to match 
several target models with very high FF) will be relatively cheap with respect to the size of the template bank (this 
picture, however, changes when we introduce frequency-domain cuts for the EP models). With the ^/jg|'s we then 
computed the number of P and EP templates necessary to cover the parameter ranges (3 : (4, 24), 77 : (0.15, 1.00), and 
(3 : (4,24), 77 : (0.1, 1.00) which span comfortably all the projected images of the target spaces onto the P and EP 
template spaces, respectively. [Note the ranges include also BBHs where one of the BH has a mass less than 5M Q .] 
We obtained 
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FIG. 28: Projection of 2PN waveforms onto the EP(3, 3.5) effective template space. This projection includes the effect of the 
frequency cut. Dots are shown for the same BBH masses of Tab. IXYI and for PN models T(2, 2.5), P(2, 2.5), ET(2, 2.5), and 
EP(2, 2.5). The thin solid lines show the BH mass lines (introduced in Sec. I VI 51 . while the dashed and dotted lines show the 
contours of the projections of selected PN models. 
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FIG. 29: Projection of 3PN waveforms onto the EP(3, 3.5) effective template space. This projection includes the effect of 
the frequency cut. Dots are shown for the same BBH masses of Tab. IXV1 and for PN models T(3,3.5, +2), P(3,3.5, +2), 
ET(3,3.5, +2), EP(3,3.5,+2), HT(3, 3.5, +2), and HP(3, 3.5, +2). The dots fo r 9 = -2 are only slightly displaced, and they 
are not shown. The thin solid lines show the BH mass lines (introduced in Sec. IV! B"t . while the dashed and dotted lines show 
the contours of the projections of selected PN models. 



( 0.02 \ , / 0.02 \ , x 

Mp * 3260 (l-MM ) ' 6700 (iTMM ) ' < 131 > 

where MM is the required minimum match (analog to the parameter MM^ of the preceding section). By comparison, 
these numbers are reduced to respectively 1230 and 3415 if we restrict 77 to the physical range. 

The number Me does not include the effect of multiple ending frequencies (cuts). We estimate the number of 
distinct / cu t needed for each (3 by an argument similar to the one used for the Fourier-domain effective templates 
(see Sec. IVIII : it turns out that more cuts are required for higher masses. The resulting number of templates is 
Nec — 51,000 for MM = 0.98, which is comparable to the result for the effective Fourier-domain templates. 

If we assume that the distance between the time-domain templates and the target models is representative of the 
distance to the true physical signal, we can guess that FF > 0.95 for P and FF > 0.97 for EP with cuts. Under these 
hypotheses, 6,500 P templates can buy us a (worst-case) MMt — 0.94, corresponding to a loss in event rate of ~ 17%. 
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FIG. 30: Determinant of the mismatch metric for the P(2, 2.5) models [left panel], and for the EP(3, 3.5, 0) models [right panel]. 
The determinant \/\g\ is shown as a function of rj and j3 = Mr/ 2 ' 5 . 



For 51,000 EP templates, we get MMt — 0.95, corresponding to a loss in event rate of ~ 14%. 

Before ending this section we would like to point out another time-domain detection-template family which can be 
consider kindred of the Fourier-domain detection-template family introduced in Sec. EH see Eq. lfTU7J) . We can use. 
for example, the following expression suggested by PN calculations [see, e.g., Ref. 62]] 



h cS (t) = Aj s (t) e^«W 



where 



•45f(*) = (* 



{t c - t) 



0(t cu t — t) , 



(132) 



(133) 



iM*) = 0c + (i c -t) 5/8 [V»0 +il>y 2 (tc-t)- l l % +il[(t c -t)- 1 l i 



(134) 



where <t> c ,t c , ct 1 > V'o > V'f) ^3/2 anc ^ V'J are arbitrary parameters whose range of values are determined maximizing the 
matches with the target two-body models. 



VIII. SUMMARY 



This paper deals with the problem of detecting GWs from the most promising sources for ground-based GW 
interferometers: comparable-mass BBHs with total mass M = 10-40M Q moving on quasicircular orbits. The detection 
of these sources poses a delicate problem, because their transition from the adiabatic phase to the plunge, at least in 
the nonspinning case, is expected to occur in the LIGO and VIRGO frequency bands. Of course, the true GW signals 
from these inspirals should be obtained from exact solutions of the Einstein equations for two bodies of comparable 
mass. However, the theoretical templates used to search for these signals will be, at best, finite-order approximations 
to the exact solutions, usually derived in the PN formalism. Because the perturbative PN approach begins to fail 
during the final stages of the inspiral, when stron g cu rvature and nonlinear effects can no longer be neglected, various 
PN resummation methods have been introduced |l4l ilaL Il6| to improve the convergence of the PN series. 

In the first part of this paper [see Sec. II I II IIVI and we studied and compared in detail all the PN models of 
the relativistic two-body dynamics currently available, including PN Taylor-expanded and resummed models both in 
the adiabatic approximation and in the nonadiabatic case. We noticed the following features [see Tables IXll IXII| . 
At least for PN orders n < 2.5, the target models T, P, and E have low cross matches, if the 2.5PN Taylor flux 
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is used. For example, for almost all the masses in our range, we found maxmaxFF < 0.9; the matches were much 
better only for P against E (and viceversa). However, if the 2PN Taylor flux is used the overlaps are rather high. 
At 3PN order we found much higher matches between T, P, and E, and also with the nonadiabatic model H, at 
least for masses M < 3OAf , and restricting to z% = = z%. These results make sense because at 3PN order the 
various approximations to the binding energy and the flux seem to be much closer to each other than at lower orders. 
This "closeness" of the different analytical approaches, which at 3PN order are also much closer to some examples 
of numerical quasiequilibrium BBH models |25j . was recently pointed out in Refs. |54l l65j|. On the other hand, the 
extraction of BBH parameters from a true measured signal, if done using the 3PN models, would still give a range 
of rather different estimates. However, we want to point out that for quite high masses, e.g., M = 40M Q , the 3PN 
models can have again lower overlaps, also from the point of view of detection. 

In addition, by studying the frequency-domain amplitude of the GW signals that end inside the LIGO frequency 
band [see Figs. 01 [7| El El i we understood that if high matches are required it is crucial to reproduce their deviations 
from the Newtonian amplitude evolution, f~ 7 ^ 6 (on the contrary, the Newtonian formula seems relatively adequate 
to model the PN amplitude for GW frequencies below the instantaneous GW frequency at the endpoint of orbital 
evolution). 

Finally, the introduction of the HT, HP and L models in Sees . IIV Al IIV Bl provided another example of two-body 
nonadiabatic dynamics, quite different from the E models. In the H models, the conservative dynamics does not have 
an ISCO [see the discussion below Eq. I|73l) ] at 2PN and 3PN orders. As a consequence, the transition to the plunge 
is due to secular radiation-reaction effects, and it is pushed to much higher frequencies. This means that, for the 
H models, the GW signals for BBHs of total mass M = 10-40Mq end outside the LIGO frequency band, and the 
frequency-domain amplitude does not deviate much from the Newtonian result, at least until very high frequencies 
[see Fig. EI]. The L models do not provide the waveforms during the late inspiral and plunge. This is due to the fact 
that because of the appearance of unphysical effects, e.g., the binding energy starts to increase with time instead of 
continuing decreasing, we are obliged to stop the evolution before the two BHs enter the last stages of inspiral. It is 
important to point out that differently from the nonadiabatic E models, the nonadiabatic H and L models give rather 
different predictions when used at various PN orders. So, from these point of view they are less reliable and robust 
than the E models. 

In the second part of this paper [Sees. IVII IVII| we pursued the following strategy. We assumed that the target 
models spanned a region in signal space that (almost) included the true GW signal. We were then able to provide a 
few detection template families (either chosen among the time-domain target models, or built directly from polynomial 
amplitude and phasings in the frequency domain) that approximate quite well all the targets [FF > 0.95 for almost all 
the masses in our range, with much better FFs for most masses]. We speculate that the effectualness of the detection 
model in approximating the targets is indicative of its effectualness in approximating the true signals. 

The Fourier-domain detection template family, discussed in Sec. IVII is simple and versatile. It uses a PN polynomial 
structure for the frequency-domain amplitude and phasing, but it does not constrain the coefficients to the PN 
functional dependencies on the physical parameters. In this sense this bank follows the basic idea that underlies the 
Fast Chirp Transform 63]. However, because for the masses that we consider the GW signal can end within the 
LIGO frequency band, we were forced to modify the Newtonian-order formula for the amplitude, introducing a cutoff 
frequency and a parameter to modify the shape of the amplitude curve (the parameter a) . As discussed at the end of 
Scc. lVIFl the good performance of the two and three-dimensional families confirms also results obtained in Refs. |12| . 
and 0. 

We showed that our Fourier-domain detection template space has a FF higher than 0.97 for the T, P and E models, 
and > 0.96 for most of the 3PN HT and HP models; we then speculate that it will match true BBH waveforms 
with FF ~ 0.96. We have computed the number of templates required to give MM ~ 0.96 (about 10 4 ). The total 
MM T should be larger than FF • MM ~ 0.92, which corresponds to a loss of event rate of I — MM| « 22%. This 
performance could be improved at the price of introducing a larger number of templates, with the rough scaling law 
of N = 1 4 [0 4/(0.96- MM)] 2 . 

In Sec. IVlEl we investigated where the less reliable 2PN H and L models, and the E models at 3PN order further 
expanded considering z\ ^= 0, lie in the detection template space. The Fourier-domain template family has FF in 
the range [0.85,0.95] with the 2PN H models, and FF mostly higher than 0.95, but with several exceptions which 
can be as low as 0.93 with the L models. The E models with z\ ^ are matched by the detection template family 
with FF almost always higher than 0.95. The E models with z\ ^ and the L models are (almost) covered by the 
region delimiting the adiabatic models and the E models with Z\ = 0. However, these models require lower cutoff 
frequencies, which will increase the number of templates up to a factor of 3. The 2PN H models sit outside this region 
and if we want to include them the number of templates should be doubled. 

The time-domain detection template families, discussed in Sec. IVIII followed a slightly different philosophy. The 
idea in this case was to provide a template bank that, for some choices of the parameters, could coincide with one 
of the approximate two-body models. Quite interestingly, this can be achieved by relaxing the physical hypothesis 
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that < i] < 0.25. However, the good performances of these banks are less systematic, and harder to generalize 
than the performance of the Fourier-domain effective bank. As suggested at the end of Sec. IV11I [see Eq. <|132[) ]. the 
time-domain bank could be improved by using a parametrization of the time-domain amplitude and phase similar to 
the one used for the Fourier-domain templates. The detection template families based on the extension of the P(2, 2.5) 
and EP(3, 3.5) to nonphysical values of r\ were shown to have FF respectively > 0.95 and > 0.97 for all the PN target 
models, and considerably higher for most models and masses. We have computed the number of P templates needed 
to obtain a MM = 0.99 (about 6,500) and of EP templates to obtain a MM = 0.98 (about 51,000). The expected 
total MMt is then respectively > 0.94 and > 0.95, corresponding to losses in event rate of < 17% and < 14%. The 
MMs scale roughly as [0.01/(1 - MM)] for P and [0.02/(1 - MM)] 2 for EP (because of the additional frequency-cut 
parameter) . 

We notice that the number of templates that we estimate for the Fourier- and time-domain detection template 
families is higher than the number of templates we would obtain using only one PN model. However, the number of 
independent shapes that enters the expression for the p„ threshold [see Eq. I18j does not coincide with the number 
of templates that are laid down within a discrete template bank to achieve a given MM; indeed, if MM is close to 
one, these are almost guaranteed to be to yield S/N statistics that are strongly correlated. A rough estimate of the 
number of independent shapes can be obtained taking a coarse-grained grid in template space. For example by setting 
MM=0 in Eq. l|25|l . the number of independent shapes would be given roughly by the volume of the template space. 
As explained at the end of Sec. Ill Bl if we wish to keep the same false-alarm probability, we have to increase the 
threshold by ~ 3% if we increase the number of independent shapes by one order of magnitude. This effect will cause 
a further loss in event rates 66]. 

Finally, in Sec. IV! Fl we extended the detection template family in the Fourier domain by requiring that it embeds 
the targets in a signal space of higher dimension (with more parameters) . We investigated the three dimensional case 
and we found, as expected, the maxmax matches increase. In particular, the match of the phasings are nearly perfect: 
always higher than 0.994 for the two-body models which are farthest apart in the detection template space. Moreover, 
by projecting the points in the three-dimensional space back to the two-dimensional space, we can get nearly the same 
projections we would have got from matching directly the PN waveforms with the two-parameter-phasing model. The 
analysis done in Sec. I VI Fl could suggest ways of systematically expand the Fourier-domain templates. Trying to guess 
the functional directions in which the true signals might lie with respect to the targets was the most delicate challenge 
of our investigation. However, our suggestions are not guaranteed to produce templates that will capture the true 
signal, and they should be considered as indications. When numerical relativity provides the first good examples 
of waveforms emitted in the last stages of the binary inspiral and plunge, it will be very interesting to investigate 
whether the matches with our detection template families are high and in which region of the detection template 
space do they sit. 
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T(2, 2) T(2, 2.5) T(3, 3.5, 0) P(2, 2.5) P(3, 3.5, 0) 

mm M rj mm M rj mm M rj mm M i] mm M i] 



T(2, 2) (2O+2O)Ai 
(15+15)M 
(15+5)M 
(5+5)M 




0.924 54.47 0.23 
0.873 39.46 0.24 
0.885 29.45 0.10 
0.988 21.28 0.06 


0.999 40.47 0.24 
0.999 30.35 0.24 
0.998 19.64 0.19 
0.998 10.61 0.22 


0.977 39.13 0.25 
0.980 29.69 0.25 
0.992 18.07 0.22 
0.994 10.54 0.22 


0.999 41.93 0.24 
0.998 31.54 0.23 
0.998 20.23 0.18 
0.999 11.16 0.20 


T(2, 2.5) (2O+2O)A/ 
(15+15)M 
(15+5)M© 
(5+5) Mq 


0.882 31.44 0.25 
0.845 24.85 0.25 
0.848 15.34 0.25 
0.801 9.41 0.25 




0.870 31.54 0.25 
0.835 25.21 0.25 
0.865 15.74 0.25 
0.823 9.51 0.25 


0.824 30.25 0.25 
0.796 25.35 0.25 
0.870 15.85 0.25 
0.826 9.51 0.25 


0.893 33.09 0.25 
0.863 26.20 0.25 
0.894 15.90 0.25 
0.849 9.61 0.25 


T(3, 3.5, 0) (2O+2O)M 
(15+15)M© 
(15+5)M© 
(5+5)Mr> 


0.999 39.57 0.24 
0.999 29.71 0.24 
0.999 20.98 0.16 
0.991 9.67 0.25 


0.916 54.63 0.23 
0.855 39.46 0.24 
0.877 29.20 0.10 
0.986 19.49 0.07 




0.989 39.03 0.24 
0.992 29.25 0.25 
0.997 18.82 0.20 
0.998 9.90 0.24 


0.997 41.56 0.23 
1.000 31.97 0.21 
1.000 20.81 0.17 
1.000 10.57 0.22 


P(2, 2.5) (2O+2O)M 
(15+15)M 
(15+5)M 
(5+5)M 


0.970 40.47 0.24 
0.967 30.15 0.24 
0.989 23.77 0.12 
0.989 9.67 0.25 


0.879 56.77 0.23 
0.816 39.66 0.24 
0.792 20.56 0.20 
0.882 13.04 0.15 


0.991 41.80 0.22 
0.998 32.66 0.20 
0.996 21.55 0.15 
0.998 10.08 0.24 




0.999 46.01 0.18 
0.999 34.02 0.19 
0.998 21.83 0.15 
0.997 10.75 0.21 


P(3, 3.5, 0) (2O+2O)M 
(15+15)M 
(15+5)M 

(5+5) Mn 


0.999 38.33 0.24 
0.997 28.47 0.25 
0.997 19.53 0.18 
0.949 9.80 0.24 


0.923 51.51 0.24 
0.979 51.01 0.10 
0.825 20.89 0.19 
0.988 17.70 0.09 


0.997 38.97 0.24 
0.997 28.96 0.25 
1.000 19.12 0.19 
0.993 9.75 0.25 


0.971 37.70 0.25 
0.961 28.88 0.25 
0.998 18.32 0.21 
0.991 9.75 0.25 




EP(2, 2.5) (2O+2O)M 
(15+15)M 
(15+5)M 
(5+5)M 


0.954 38.10 0.25 
0.965 29.34 0.25 
0.988 20.79 0.16 
0.996 9.70 0.25 


0.936 51.14 0.24 
0.895 37.45 0.25 
0.769 21.97 0.19 
0.980 20.46 0.07 


0.933 39.10 0.25 
0.960 29.60 0.25 
0.983 20.22 0.18 
0.997 10.29 0.23 


0.878 38.22 0.25 
0.903 29.56 0.25 
0.969 19.54 0.19 
0.995 10.22 0.23 


0.962 39.94 0.25 
0.975 30.15 0.25 
0.980 20.85 0.17 
0.997 10.83 0.21 


EP(3, 3.5, 0) (2O+2O)M 
(15+15)M Q 
(15+5)M© 
(5+5)M 


0.946 37.11 0.25 
0.955 28.78 0.24 
0.992 18.51 0.20 
0.968 9.65 0.25 


0.949 48.90 0.24 
0.913 35.38 0.24 
0.808 22.15 0.18 
0.985 18.41 0.08 


0.930 37.84 0.25 
0.948 28.89 0.25 
0.985 18.92 0.20 
0.994 9.76 0.25 


0.867 36.72 0.25 
0.893 28.82 0.25 
0.970 18.34 0.21 
0.992 9.77 0.25 


0.954 38.80 0.24 
0.968 29.50 0.25 
0.983 19.63 0.19 
0.998 10.16 0.23 


HT(2, 2) (2O+2O)M 
(15+15)M 
(15+5)M 
(5+5)M 


0.777 21.39 0.25 
0.674 20.20 0.24 
0.616 15.88 0.20 
0.796 9.62 0.25 


0.890 27.58 0.25 
0.780 21.83 0.25 
0.666 18.84 0.18 
0.935 10.00 0.25 


0.768 21.61 0.25 
0.673 21.02 0.25 
0.625 17.37 0.18 
0.833 9.73 0.25 


0.732 21.63 0.25 
0.657 21.03 0.25 
0.645 16.10 0.22 
0.834 9.74 0.25 


0.789 22.57 0.25 
0.687 21.07 0.25 
0.631 17.14 0.18 
0.856 9.75 0.25 


HT(3, 3.5, 0) (2O+2O)M 
(15+15)Af 
(15+5)M 
(5+5)M 


0.812 32.35 0.25 
0.848 27.97 0.25 
0.998 23.08 0.13 
0.952 9.65 0.25 


0.925 44.91 0.24 
0.919 33.30 0.25 
0.788 21.15 0.20 
0.828 10.36 0.24 


0.795 34.76 0.25 
0.835 28.70 0.25 
0.999 21.25 0.16 
0.984 9.76 0.25 


0.737 32.98 0.25 
0.788 28.78 0.25 
0.994 19.77 0.18 
0.984 9.77 0.25 


0.812 37.10 0.24 
0.875 29.07 0.25 
0.999 21.81 0.15 
0.992 9.99 0.24 


HP(2, 2.5) (2O+2O)M 
(15+15)M 
(15+5)M 
(5+5)M 


0.756 18.71 0.25 
0.631 17.87 0.24 
0.582 14.33 0.25 
0.731 9.41 0.25 


0.853 23.74 0.24 
0.714 18.06 0.25 
0.631 16.88 0.20 
0.869 9.75 0.25 


0.752 18.96 0.25 
0.634 17.86 0.25 
0.587 14.54 0.25 
0.755 9.51 0.25 


0.725 19.09 0.25 
0.630 18.46 0.25 
0.600 16.40 0.18 
0.755 9.54 0.25 


0.769 19.70 0.25 
0.642 18.53 0.25 
0.589 17.88 0.15 
0.765 9.54 0.25 


HP(3, 3.5, 0) (2O+2O)M 
(15+15)Af 
(15+5)M 
(5+5)M 


0.748 32.36 0.25 
0.789 27.41 0.24 
0.998 21.75 0.15 
0.912 9.62 0.25 


0.879 42.53 0.25 
0.915 31.80 0.25 
0.792 20.41 0.21 
0.990 16.20 0.10 


0.733 32.51 0.25 
0.782 27.43 0.25 
1.000 20.57 0.17 
0.959 9.73 0.25 


0.679 30.72 0.25 
0.741 27.43 0.25 
0.995 19.29 0.19 
0.961 9.76 0.25 


0.756 34.48 0.25 
0.817 28.60 0.25 
0.999 21.17 0.16 
0.982 9.76 0.25 



TABLE XI: (Continued into Table IxTTl 'l Fitting factors between several PN models, at 2PN and 3PN orders. For three choices 
of BBH masses, this table shows the maxmax matches [see Eq. (llOI ll between the search models at the top of the columns and 
the target models at the left end of the rows, maximized over the intrinsic parameters of the search models in the columns. For 
each intersection, the three numbers mm, M = mi + 7712 and rj = mi 772,2 /M 2 denote the maximized match and the search-model 
mass parameters at which the maximum is attained. In computing these matches, the parameter r\ of the search models was 
restricted to its physical range < 77 < 1/4. The arbitrary flux parameter 9 was always set equal to zero. 

These matches represent the fitting factors [see Eq. 120H 1 for the projection of the target models onto the search models. The 
reader will notice that the values shown are not symmetric across the diagonal: for instance, the match for the search model 
T(2, 2.5) against the target model P(2, 2.5) is higher than the converse. This is because the matches represent the inner product 
Q between two different pairs of model parameters: in the first case, the target parameters (mi = 15M0,m2 = 15Mq)p = 
(M = 30Mo,77 = 0.25)p are mapped to the maximum- match search parameters (M = 39. 7M©, 77 = 0.24)t; in the second case, 
the target parameters (mi = 15M0,m.2 = 15Mq)t = (M = 30Mq,t] — 0.25)t are mapped to the maximum-match parameters 
(M = 25.37M0,r7 = 0.24)p [so the symmetry of the inner product Q is reflected by the fact that the search parameters 
(M = 25.3M©,7i = 0.24)p are mapped into the target parameters (M = 3OM0,77 = 0.25)t]- 
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EP 
mm 



(2, 2.5) 

M J7 



EP(3, 3.5, 
mm M 



0) 

V 



HT(2, 2) 
mm M ij 



HT(3, 3.5, 
mm M 



0) 
»7 



HP(2, 2.5) 
mm M rj 



HP(3, 3.5, 
mm M 



0) 



T(2, 2) 



T(2, 2.5) 



(20+20 
(15+15 
(15+5 
(5+5 



M 3 
Mq 
M 
Mo 



0.953 
0.962 
0.988 
0.997 



41.67 
30.41 
19.11 
10.33 



0.24 
0.24 
0.20 
0.23 



0.952 
0.991 
0.992 
0.998 



43.00 
35.32 
20.93 
11.09 



0.24 
0.17 
0.17 
0.20 



0.951 
0.899 
0.924 
0.788 



80.34 
58.93 
69.96 
9.93 



0.24 
0.24 
0.05 
0.25 



0.855 
0.997 
0.998 
0.998 



56.69 
33.03 
19.38 
10.92 



0.24 
0.20 
0.19 
0.21 



0.965 
0.922 
0.876 
0.727 



90.12 
67.38 
57.94 
10.19 



0.24 
0.24 
0.07 
0.25 



0.859 
0.998 
0.999 
0.999 



74.80 
33.67 
19.81 
11.19 



0.25 
0.20 
0.18 
0.20 



T(3, 3.5, 0) 



(20+20 
(15+15 
(15+5 
(5+5 



Mq 
M 
Mq 
W 



0.908 
0.861 
0.822 
0.814 



31.37 
24.52 
15.40 
9.52 



0.25 
0.25 
0.25 
0.25 



0.929 
0.893 
0.867 
0.839 



32.98 
25.58 
15.81 
9.59 



0.25 
0.25 
0.25 
0.25 



0.959 
0.932 
0.790 
0.941 



58.39 
53.46 
16.59 
9.63 



0.24 
0.17 
0.25 
0.25 



0.928 
0.926 
0.903 
0.838 



35.74 
26.82 
15.81 
9.52 



0.24 
0.25 
0.25 
0.25 



0.955 
0.920 
0.839 
0.872 



67.85 
51.38 
51.91 
9.80 



0.24 
0.24 
0.07 
0.25 



0.892 
0.921 
0.955 
0.866 



36.87 
27.99 
16.03 
9.61 



0.23 
0.24 
0.25 
0.25 



P(2, 2.5) 



(20+20 
(15+15 
(15+5 
(5+5 



Mq 
Mq 
Mq 

.\/. 



0.925 
0.955 
0.983 
0.992 



40.09 
29.98 
19.68 
9.99 



0.24 
0.24 
0.18 
0.24 



0.918 
0.937 
0.985 
0.997 



42.90 
30.78 
20.97 
10.40 



0.24 
0.24 
0.16 
0.22 



0.940 
0.887 
0.926 
0.826 



80.76 
58.83 
69.81 
9.83 



0.24 
0.24 
0.05 
0.25 



0.833 
0.996 
0.999 
0.993 



57.71 
32.67 
19.47 
10.48 



0.24 
0.20 
0.19 
0.22 



0.958 
0.914 
0.887 
0.749 



89.85 
66.56 
60.02 
10.07 



0.24 
0.24 
0.07 
0.25 



0.840 
0.758 
1.000 
0.995 



73.84 
31.32 
19.79 
10.81 



0.25 
0.24 
0.18 
0.21 



P(3, 3.5, 0) 



(20+20 
(15+15 
(15+5 
(5+5 



Mq 

Mq 

Mq 
Mq 



0.866 
0.898 
0.966 
0.995 



41.72 
30.06 
20.48 
9.79 



0.24 
0.24 
0.17 
0.25 



0.859 
0.963 
0.966 
0.994 



43.14 
38.21 
21.86 
10.43 



0.24 
0.14 
0.15 
0.22 



0.912 
0.857 
0.907 
0.825 



83.09 
62.07 
70.42 

9.81 



0.24 
0.24 
0.05 
0.25 



0.795 
0.992 
0.993 
0.990 



65.45 
33.28 
20.08 
10.51 



0.24 
0.19 
0.17 
0.22 



0.934 
0.890 
0.904 
0.748 



92.91 
69.31 
64.71 
10.05 



0.24 
0.24 
0.06 
0.25 



0.805 
0.709 
0.997 
0.992 



82.71 
59.88 
20.29 
10.83 



0.25 
0.25 
0.17 
0.21 



EP(2, 2.5) 



(20+20 
(15+15 
(15+5 
(5+5 



Mq 
Mq 
Mq 
Mq 



0.960 
0.965 
0.982 
0.973 



40.10 
29.33 
18.87 
9.74 



0.23 
0.24 
0.20 
0.25 



0.953 
0.966 
0.983 
0.998 



41.06 
30.14 
20.29 

9.85 



0.24 
0.24 
0.17 
0.25 



0.943 
0.893 
0.926 
0.849 



76.61 
56.29 
68.98 
9.81 



0.24 
0.24 
0.05 
0.25 



0.835 
0.993 
0.996 
0.992 



53.85 
31.83 
19.15 
10.02 



0.24 
0.20 
0.19 
0.24 



0.961 
0.920 
0.886 
0.761 



86.56 
63.91 
58.97 
10.04 



0.24 
0.24 
0.07 
0.25 



0.842 
0.996 
0.999 
0.993 



70.76 
32.41 
19.45 
10.46 



0.25 
0.20 
0.19 
0.22 



EP(3, 3.5, 0) 



(20+20 
(15+15 
(15+5 
(5+5 



Mq 
Mq 
Mq 



0.996 
0.999 
0.999 
0.999 



41.72 
32.66 
21.35 
10.75 



0.24 
0.21 
0.16 
0.21 



0.953 
0.908 
0.909 
0.807 



75.09 
56.68 
70.41 

9.84 



0.24 
0.24 
0.05 
0.25 



0.929 
0.889 
0.992 
0.997 



47.51 
32.89 
19.52 
10.69 



0.24 
0.24 
0.19 
0.21 



0.948 
0.915 
0.858 
0.733 



84.61 
64.87 
64.23 
10.08 



0.24 
0.24 
0.06 
0.25 



0.907 
0.997 
0.986 
0.998 



59.72 
33.00 
20.00 
10.99 



0.24 
0.20 
0.18 
0.20 



HT(2, 2) 



(20+20 
(15+15 
(15+5 
(5+5 



Mq 
Mq 
Mq 



0.995 
0.992 
0.999 
0.982 



38.25 
28.77 
18.53 
9.74 



0.25 
0.25 
0.20 
0.25 



0.958 
0.938 
0.905 
0.832 



72.99 
70.37 
69.04 
10.00 



0.24 
0.14 
0.05 
0.24 



0.918 
0.999 
0.998 
0.996 



45.74 
31.41 
18.97 
10.24 



0.24 
0.21 
0.20 
0.23 



0.956 
0.922 
0.858 
0.748 



81.66 
61.77 
61.43 
10.06 



0.24 
0.24 
0.06 
0.25 



0.896 
1.000 
0.994 
0.997 



59.30 
32.11 
19.26 
10.61 



0.25 
0.21 
0.19 
0.22 



HT(3, 3.5, 0) 



(20+20 
(15+15 
(15+5 
(5+5 



Mq 
Mq 
Mq 



0.794 
0.651 
0.624 
0.817 



21.34 
18.40 
14.96 
9.72 



0.25 
0.24 
0.25 
0.25 



0.815 
0.674 
0.632 
0.845 



22.35 
19.03 
15.15 
9.74 



0.25 
0.24 
0.25 
0.25 



0.840 
0.377 
0.608 
0.845 



24.31 
37.58 
17.70 
9.74 



0.25 
0.25 
0.17 
0.25 



0.968 
0.936 
0.965 
0.841 



46.75 
36.99 
17.85 
9.97 



0.25 
0.24 
0.22 
0.25 



0.835 
0.392 
0.612 
0.865 



25.77 
47.22 
17.35 
9.76 



0.25 
0.25 
0.18 
0.25 



HP(2, 2.5) 



(20+20 
(15+15 
(15+5 
(5+5 



Mq 
Mq 
Mq 
M. 



0.904 
0.891 
0.986 
0.964 



34.61 
27.49 
20.73 
9.75 



0.24 
0.25 
0.16 
0.25 



0.920 
0.926 
0.986 
0.993 



37.64 
28.59 
21.99 
9.79 



0.24 
0.25 
0.15 
0.25 



0.903 
0.883 
0.919 
0.834 



65.68 
49.56 
71.02 
9.83 



0.24 
0.24 
0.05 
0.25 



0.873 
0.867 
0.886 
0.749 



74.44 
59.23 
61.90 
10.07 



0.25 
0.24 
0.07 
0.25 



0.999 
1.000 
1.000 
1.000 



41.41 
31.02 
20.34 
10.35 



0.23 
0.23 
0.17 
0.23 



HP(3, 3.5, 0) 



(20+20 
(15+15 
(15+5 
(5+5 



Mq 
Mq 
Mq 



0.762 
0.595 
0.577 
0.741 



18.74 
16.37 
16.04 
9.50 



0.25 
0.24 
0.20 
0.25 



0.784 
0.617 
0.599 
0.754 



19.44 
16.40 
14.32 
9.53 



0.25 
0.24 
0.25 
0.25 



0.973 
0.931 
0.957 
0.975 



36.64 
27.84 
22.10 
11.46 



0.21 
0.21 
0.14 
0.18 



0.794 20.75 0.24 
0.329 40.09 0.25 
0.589 15.53 0.21 
0.755 9.52 0.25 



0.801 
0.343 
0.593 
0.770 



21.53 
48.60 
15.59 
9.61 



0.25 
0.25 
0.21 
0.25 



(20+20 
(15+15 
(15+5 
(5+5 



Mq 
Mq 
Mq 

M_ 



0.832 
0.831 
0.986 
0.933 



31.43 
26.96 
20.13 
9.72 



0.25 
0.25 
0.17 
0.25 



0.840 
0.860 
0.986 
0.971 



35.15 
28.03 
21.50 
9.75 



0.25 
0.25 
0.15 
0.25 



0.850 
0.852 
0.922 
0.857 



60.63 
46.65 
70.24 
9.80 



0.25 
0.24 
0.05 
0.25 



0.974 37.71 0.25 
0.975 28.95 0.25 
1.000 19.64 0.18 
0.991 9.75 0.25 



0.806 72.61 0.25 

0.842 55.71 0.24 

0.884 60.67 0.07 

0.758 10.03 0.25 



TABLE XII: (Continued from Table EH) 
caption to Tabic |X"TI 



Fitting factors between several PN models, at 2PN and 3PN orders. Please see the 



T(2, 2.5) ET(2,2.5) T(3, 3.5, +2) T(3,3.5, -2) ET(3,3.5,+2) ET(3, 3.5, -2) 

mm M rj mm M rj mm M rj mm M rj mm M rj mm M rj 



(15+15)M a 
T(2,2.5) (15+5)M 
(5+5)M 




0.914 27.58 0.248 
0.916 16.81 0.249 
0.900 10.13 0.241 










(15+15)M 
ET(2,2.5) (15+5)M 
(5+5)M 


0.922 33.93 0.241 
0.971 33.17 0.076 
0.984 13.57 0.147 












(16+15)Af s 
T(3,3.5,+2) (15+5)M 

(5 + 5) Mq 








0.995 29.83 0.243 
1.000 19.06 0.204 
0.981 9.96 0.250 


0.963 30.52 0.240 
0.984 20.03 0.186 
0.991 10.16 0.242 


0.974 30.32 0.240 
0.974 20.09 0.182 
0.972 9.94 0.250 


(15+15)M 
T(3,3.5,-2) (15+5)M 
(5+5)M 






0.998 30.94 0.242 
1.000 20.93 0.173 
0.999 10.61 0.226 




0.951 31.27 0.239 
0.985 20.89 0.173 
0.994 10.26 0.240 


0.960 30.59 0.241 
0.983 20.27 0.181 
0.993 10.19 0.241 


(15+15)M 
ET(3,3.5,+2) (15+5)M 
(5+5)M 






0.951 30.39 0.240 
0.981 20.16 0.186 
0.996 10.22 0.240 


0.931 29.76 0.241 
0.985 18.97 0.207 
0.985 9.96 0.250 




0.994 30.06 0.241 
1.000 19.23 0.201 
0.979 9.95 0.250 


(16+15)Af s 
ET(3,3.5,-2) (15+5)M 

(5 + 5) Mq 






0.963 30.94 0.240 
0.983 20.65 0.179 
0.987 10.27 0.240 


0.953 30.30 0.241 
0.980 20.32 0.182 
0.996 10.21 0.241 


0.999 31.07 0.238 
1.000 20.83 0.175 
1.000 10.51 0.230 





TABLE XIII: Fitting factors between T and ET models, at 2PN and 3PN orders, and for different choices of the arbitrary flux 
parameter 8. For three choices of BBH masses, this table shows the maxmax matches [see Eq. HUH ] between the search models 
at the top of the columns and the target models at the left end of the rows, maximized over the mass parameters of the models 
in the columns. For each intersection, the three numbers mm, M and n denote the maximized match and the search-model 
mass parameters at which the maximum is attained. The matches can be interpreted as the fitting factors for the projection 
of the target models onto the search models. See the caption to Table IXTT1 for further details. 
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EP(3,3.5, 2, -4, 0) EP(3, 3.5, 2, 0, -4) EP(3, 3.5, 2, 0, 0) EP(3, 3.5, 2, 0, 4) EP(3, 3.5, 2, 4, 0) 
mm Af ?y mm M rj mm M r/ mm M mm M rj 



(15+15)M 
EP(3,3.5,2, -4,0) (15+5)M 
(5+5)M 




0.995 30.93 0.238 
0.998 20.61 0.177 
0.999 10.22 0.240 


0.994 30.85 0.240 
0.999 20.71 0.176 
0.999 10.22 0.240 


0.995 30.87 0.239 
0.999 20.60 0.177 
0.999 10.22 0.240 


0.952 31.17 0.242 
0.993 21.59 0.162 
0.996 10.46 0.231 


(15+15)AT© 
EP(3, 3.5,2,0, -4) (15+5)M 
(5+5)M 


0.983 30.12 0.241 
0.999 19.28 0.201 
0.993 10.01 0.249 




0.999 30.47 0.240 
1.000 20.06 0.186 
0.996 10.19 0.241 


0.999 30.43 0.241 
1.000 20.03 0.187 
0.996 10.19 0.241 


0.987 30.88 0.240 
0.999 20.70 0.175 
0.998 10.22 0.240 


(15+15)M 
EP(3, 3.5, 2, 0, 0) (15+5)M 
(5+5)M 


0.983 30.12 0.241 
0.999 19.26 0.202 
0.993 9.99 0.250 


0.999 30.47 0.241 
1.000 20.06 0.186 
1.000 10.00 0.250 




0.999 30.42 0.241 
1.000 20.03 0.187 
0.996 10.19 0.241 


0.987 30.88 0.240 
0.999 20.70 0.175 
0.998 10.22 0.240 


(15+15)Af© 
EP(3, 3.5, 2, 0, 4) (15+5)M 
(5+5)M 


0.982 30.12 0.241 
0.999 19.35 0.200 
0.993 10.01 0.249 


0.999 30.54 0.240 
1.000 20.05 0.187 
1.000 10.00 0.250 


0.999 30.54 0.240 
1.000 19.98 0.188 
0.996 10.19 0.241 




0.987 30.88 0.240 
0.998 20.73 0.175 
0.998 10.22 0.240 


(15+15)M 
EP(3, 3.5,2,4,0) (15+5)M 
(5+5)M 


0.929 29.60 0.240 
0.992 18.42 0.219 
0.970 10.17 0.241 


0.968 30.11 0.242 
0.998 19.29 0.201 
0.993 9.99 0.250 


0.968 30.16 0.240 
0.998 19.36 0.199 
0.993 9.99 0.250 


0.967 30.15 0.240 
0.998 19.29 0.201 
0.993 9.99 0.250 





TABLE XIV: Fitting factors for the projection of EP(3, 3.5, 0) templates onto themselves, for various choices of the parameters 
Zi and Z2- The values quoted are obtained by maximizing the maxmax (mm) match over the mass parameters of the (search) 
models in the columns, while keeping the mass parameters of the (target) models in the rows fixed to their quoted values, 
(15 + 15)M@, (15 + 5)Mq (5 + 5)Mq. The three numbers shown at each intersection are the maximized match and the search 
parameters at which the maximum was attained. In labeling rows and columns we use the notation EP(3, 3.5, 8, Zi, z-i). See 
the caption to Table IxTTI for further details. 



FF for projection onto the Fourier-domain detection template families 



PN model 




/end mn l/>0 V'3/2 <*fc,{t /cut 




/e„d mn VO '03/2 a flit /cut 


T(2, 2) 


(2O+2O)M 
(2O+15)M 
(15+15)Af© 
(2O+1O)M 
(15+10)Af© 


221.4 0.983 23891. -554.63 0.949 240.7 
252.4 0.987 30200. -606.41 0.975 272.5 
295.2 0.989 38126. -653.61 0.968 313.5 
291.7 0.989 41735. —677.51 1.002 314.2 
352.7 0.991 52565. -713.54 0.968 387.1 


(2O+5)M 
(10+10)Af© 
(15+5)Af© 
(1O+5)M 
(5+5)M 


341.2 0.992 77508. -1041.30 0.897 347.0 
442.8 0.992 72639. -768.78 0.632 331.4 

431.3 0.993 96191. -1030.20 0.831 440.8 

583.4 0.993 130600. —1019.10 1.001 805.3 
885.6 0.989 225060. -1056.80 0.531 894.4 


T(2, 2.5) 


(2O+2O)M 
(2O+15)M 
(15+15)Af© 
(2O+1O)M 
(15+1O)M 


161.2 0.970 19807. 62.32 0.691 224.4 
185.9 0.975 25398. 57.59 0.347 220.3 
214.9 0.979 32787. 40.11 0.210 245.0 

222.3 0.980 36540. 28.23 0.160 255.5 
261.2 0.983 47008. 2.24 0.107 293.7 


(2O+5)M 
(10+10)Af© 
(15+5)M 
(lO+5)M 
(5+5)Af© 


281.6 0.987 71552. -188.92 0.227 312.7 
322.4 0.983 66783. -37.92 0.490 630.9 
345.6 0.988 89296. -166.70 0.107 373.8 
443.3 0.989 123100. —159.28 0.379 746.0 
643.9 0.994 217090. -194.81 0.253 1033.1 


T(3,3.5,+2) 


(2O+2O)M 
(20+15)Af© 
(15+15) AT© 
(2O+1O)M 
(15+10)Af© 


207.9 0.983 25219. -575.44 1.002 265.8 
234.5 0.984 31622. -623.54 1.005 268.5 

277.2 0.987 38891. -612.96 0.990 306.3 

259.3 0.986 43944. —729.80 0.979 301.6 
324.3 0.987 53869. -688.38 0.865 315.6 


(2O+5)M 
(10+10)Af© 
(15+5)M 
(1O+5)M 
(5+5)Af© 


276.1 0.986 79630. -1095.00 0.743 238.3 
415.9 0.988 73738. -701.48 0.923 437.8 
362.3 0.990 97371. -988.17 0.617 277.0 
518.5 0.990 131210. —899.96 0.642 392.3 
831.7 0.985 224370. -826.19 0.563 886.2 


T(3,3.5, -2) 


(2O+2O)M 
(20+15)Af© 
(15+15)Af© 
(20+10)Af© 
(15+10)Af© 


207.9 0.981 24857. -603.44 0.983 246.4 
234.5 0.985 31773. -681.75 0.983 252.8 

277.2 0.986 39565. -707.26 0.933 277.9 

259.3 0.985 44027. -787.96 0.900 251.9 
324.3 0.988 54194. -761.61 0.984 341.1 


(2O+5)M 
(10+10)M o 
(15+5)M 
(1O+5)M 
(5+5)M 


276.1 0.987 80359. -1188.90 0.825 257.0 
415.8 0.988 74637. -810.89 0.750 350.3 
362.3 0.989 97861. -1070.50 0.661 267.7 
518.5 0.988 131840. -992.35 0.901 553.3 
831.7 0.982 225550. -943.65 0.577 916.3 


P(2,2.5) 


(2O+2O)M 
(20+15)Af© 
(15+15) AT© 
(20+10)Af© 
(15+10)Af© 


142.9 0.972 27006. -743.88 0.991 208.5 

162.5 0.977 33307. -778.72 0.987 206.7 

190.6 0.980 40486. -752.07 0.991 237.0 
185.0 0.977 45403. -864.50 1.116 288.3 
226.3 0.981 54709. -771.73 0.867 232.9 


(2O+5)M 
(10+10)M o 
(15+5)M 
(10+5)Af© 
(5+5)M 


207.8 0.978 81397. -1244.40 0.698 192.4 

285.9 0.985 73970. -743.09 0.681 245.7 
267.5 0.984 98390. -1074.60 0.709 231.8 
370.0 0.985 131920. -961.15 0.758 346.4 
571.8 0.983 224810. -867.58 0.813 764.5 


P(3,3.5,+2) 


(20+20) M Q 
(2O+15)M 
(15+15) AT© 
(20+10)Af© 
(15+10)Af© 


216.4 0.984 24922. -523.74 0.995 265.2 

243.6 0.985 31204. -564.86 1.007 299.3 

288.5 0.987 38194. -541.27 0.971 328.2 

265.7 0.986 43280. -660.41 1.001 328.8 
336.2 0.987 52941. -605.52 0.902 356.7 


(2O+5)M 
(10+10)M o 
(15+5)M 
(1O+5)M 
(5+5)M 


265.0 0.985 79624. -1070.20 0.830 258.4 
432.8 0.990 72663. -617.31 0.896 488.2 

359.2 0.990 96933. -935.65 0.619 279.6 

531.3 0.991 130310. -827.00 0.843 588.6 
865.6 0.988 223830. -780.35 0.537 896.7 


P(3,3.5,-2) 


(20+20) AT© 
(2O+15)M 
(15+15)Af© 
(2O+1O)M 
(15+1O)M 


216.4 0.984 24830. -545.66 1.062 291.4 

243.6 0.984 31086. -583.34 0.988 269.5 

288.5 0.988 38426. -581.05 0.994 326.6 

265.7 0.986 43464. -696.77 1.006 311.2 
336.2 0.987 53475. -663.65 0.882 333.4 


(2O+5)M 
(10+10)M o 
(15+5)Af© 
(1O+5)M 
(5+5)M 


265.0 0.986 79956. -1114.80 0.831 259.7 
432.8 0.990 73167. -674.59 0.760 390.9 

359.2 0.990 96850. -958.04 0.662 277.7 

531.3 0.990 130780. -881.70 0.810 539.0 
865.6 0.987 224210. -828.64 0.538 896.0 


ET(2, 2.5) 


(2O+2O)M 
(2O+15)M 
(15+15)Af© 
(2O+1O)M 
(15+1O)M 


231.0 0.991 22372. -258.47 0.935 477.8 
263.5 0.992 28710. -302.99 0.770 425.5 

308.0 0.993 36351. -321.50 0.717 512.3 

305.1 0.993 41308. -423.25 0.756 473.1 
368.3 0.995 51338. -393.70 0.769 764.8 


(2O+5)M 
(10+10)M o 
(15+5)Af© 
(1O+5)M 
(5+5)M 


359.4 0.995 79070. -857.02 0.748 519.2 

462.0 0.995 71411. -420.76 0.668 722.3 
452.7 0.994 96788. -755.70 0.718 706.8 

610.1 0.993 129130. -607.98 0.665 910.1 
924.0 0.991 221910. -534.76 0.424 920.4 


ET(3,3.5,+2) 


(2O+2O)M 
(2O+15)M 
(15+15)Af© 
(2O+1O)M 
(15+1O)M 


212.1 0.990 22048. -356.02 0.997 367.7 
245.1 0.992 28516. -423.30 0.971 415.7 
285.8 0.992 36119. -450.40 0.775 408.3 
286.6 0.993 40717. -545.11 0.790 376.5 
344.5 0.993 51507. -563.26 0.785 515.2 


(2O+5)M 
(1O+1O)M 
(15+5)Af© 
(1O+5)M 
(5+5)M 


351.3 0.992 78355. -1057.40 0.763 402.3 

428.8 0.994 72187. -631.44 0.707 616.2 
433.7 0.994 96772. -982.67 0.757 572.4 
573.1 0.995 130830. -899.77 0.686 856.6 

847.9 0.986 225490. -892.59 0.552 914.8 


ET(3,3.5, -2) 


(2O+2O)M 
(20+15)Af© 
(15+15)Af© 
(2O+1O)M 
(15+1O)M 


207.1 0.990 21818. -386.23 0.848 300.4 

238.2 0.992 28247. -451.93 0.884 347.0 
274.0 0.992 36218. -502.72 0.903 452.4 
277.0 0.992 41148. -613.88 0.786 364.2 
330.6 0.992 51702. -623.17 0.822 501.2 


(2O+5)M 
(10+10)M o 
(15+5)M 
(1O+5)M 
(5+5)M 


345.9 0.991 78349. -1103.70 0.692 379.0 
411.0 0.994 72645. -709.64 0.685 499.9 

424.8 0.993 97086. -1052.50 0.846 600.5 

556.9 0.995 131730. -1003.30 0.699 821.5 
816.8 0.983 226430. -999.02 0.539 900.4 



■53 



PN model 



/end 



mn 



FF for 

•00 



projection 

•03/2 



onto the Fourier-domain 

/cut 
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"/cut 



detection template families 

/end mn ijlp "03/2 



„2/3 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



EP(2, 2.5) 



20+20 
20+15 
15+15 
20+10 
15+10 



M, 
M 
M 
M 



218.1 
249.1 
290.8 
289.8 
348.5 



0.991 
0.991 
0.991 
0.990 
0.991 



21315. 
28013. 
35947. 
40730. 
51920. 



-353.42 
-437.59 
-486.80 
-593.65 
-632.99 



0.773 
0.746 
0.672 
0.656 
0.637 



376.1 
380.6 
432.2 
391.0 
451.6 



(20+!: 
(10+10 
(15+5 
(10+5 
(5+5 



M 
M 
M 
M 
M 



345.8 
436.2 
433.1 
579.6 
872.5 



0.990 
0.994 
0.994 
0.993 
0.979 



79526. 

73183. 

98170. 
132250. 
226910. 



-1167.70 
-729.74 
-1099.60 
-1014.20 
-997.82 



0.709 
0.714 
0.630 
0.691 
0.571 



366.2 
645.1 
460.4 
868.7 
833.0 



EP(3,3.5, +2) 



20+20 
20+15 
15+15 
20+10 
15+10 



M 
M 
M 
M 
M 



219.7 
251.9 
293.1 
291.4 
353.7 



0.990 
0.991 
0.991 
0.993 
0.993 



22025. 
27970. 
35861. 
40598. 
51343. 



-329.13 
-368.53 
-409.25 
-512.88 
-527.79 



0.967 
0.888 
0.813 
0.820 
0.731 



398.3 
386.4 
452.7 
429.6 
495.5 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



M G 
M 
M 
M 



354.0 
439.6 
444.9 
582.5 
874.7 



0.993 
0.994 
0.995 
0.995 
0.989 



78344. 

71704. 

96416. 
130480. 
224370. 



-1027.30 
-579.45 
-934.82 
-855.36 
-820.10 



0.668 
0.719 
0.773 
0.685 
0.488 



376.6 
658.8 
608.4 
879.7 
916.2 



EP(3,3.5, -2) 



20+20 
20+15 
15+15 
20+10 
15+10 



A/, 

M 
M 
M 
Mr, 



214.4 
248.3 
287.0 
289.1 
347.3 



0.990 
0.992 
0.992 
0.993 
0.993 



22029. 
28185. 
35793. 
40653. 
51423. 



-349.92 
-400.30 
-429.31 
-537.88 
-558.41 



0.986 
0.849 
0.880 
0.869 
0.779 



384.7 
361.2 
510.6 
452.9 
494.7 



M 
M 
M 
M 



353.0 
430.5 
439.1 
575.7 
864.9 



0.992 
0.994 
0.994 
0.995 
0.988 



78099. 

71820. 

96411. 
130760. 
225110. 



-1035.10 
-613.97 
-960.71 
-899.02 
-886.01 



0.692 
0.718 
0.770 
0.696 
0.501 



400.3 
642.1 
591.0 
877.1 
909.7 



(20+5 
(10+10 
(15+5 
(10+5 



EP(3, 3.5, 0,4,0) 



20+20 
20+15 
15+15 
20+10 
15+10 



M 
M 
M 
M 
M 



318.1 
364.6 
432.0 
420.4 
510.8 



0.989 
0.988 
0.987 
0.987 
0.988 



20061. 
26379. 
34134. 
38610. 
49757. 



-192.06 
-249.89 
-293.98 
-385.18 
-426.26 



0.509 
0.437 
0.321 
0.455 
0.515 



379.7 
385.7 
422.8 
446.0 
493.4 



M e 
M @ 
M @ 

Me 



457.4 
647.2 
600.6 
831.6 
1292.2 



0.987 
0.990 
0.992 
0.995 
0.992 



76939. 

70495. 

95378. 
129410. 
223410. 



-936.06 
-502.74 
-866.93 
-792.01 



0.683 
0.585 
0.651 
0.680 
0.339 



450.0 
666.7 
601.3 
798.8 
1003.8 



EP(3, 3.5,0, -20,0) 



20+20 
20+15 
15+15 
20+10 
15+10 



M ( : } 
Me 
Me 
M @ 
Me 



118.9 
136.9 
158.5 
164.2 
192.4 



0.970 
0.983 
0.983 
0.985 
0.985 



26410. 
33451. 
41909. 
46550. 
56925. 



-787.54 
-868.80 
-921.14 
-1016.40 
-986.07 



0.964 
1.010 
1.045 
1.138 
1.096 



189.8 
238.9 
285.0 
321.2 
339.8 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



M Q 
M 
Me 
M e 
M 



215.1 
237.8 
258.0 
327.9 
476.1 



0.989 
0.983 
0.984 
0.977 
0.969 



83591. 

76684. 
101440. 
134130. 
226450. 



-1452.50 
-970.56 
-1323.20 
-1142.10 
-992.89 



1.087 
1.074 
1.158 
1.157 
1.167 



364.1 
373.8 
486.7 
589.3 
844.2 



EP(3, 3.5,0, -40,0) 



20+20 
20+15 
15+15 
20+10 
15+10 



Me 
Me 
Me 
M @ 
M Q 



94.0 
108.2 
125.3 
130.4 
152.5 



0.947 
0.962 
0.969 
0.970 
0.964 



29400. 
36837. 
45552. 
50375. 
61789. 



-1174.60 
-1268.40 
-1324.90 
-1423.70 
-1428.90 



1.097 
0.960 
1.010 
1.048 
1.077 



184.5 
169.2 
228.4 
252.8 
338.4 



M e 
Me 
Me 
Me 

Mr. 



174.0 
188.0 
206.7 
260.8 
376.1 



0.972 
0.959 
0.967 
0.957 
0.955 



88302. 

82469. 
105660. 
137720. 
228960. 



-1874.90 
-1437.30 
-1681.70 
-1431.90 
-1185.20 



1.073 
1.059 
1.357 
1.111 
1.122 



337.6 
411.8 
468.5 
537.6 
874.4 



EP(T 3 , 3.5,0, +40,0) 



20+20 
20+15 
15+15 
20+10 
15+10 



M: 

M Q 
Me 
M @ 
M @ 



349.5 
399.4 
465.3 
468.3 
558.6 



0.986 
0.989 
0.987 
0.989 
0.989 



19559. 
25098. 
32573. 
36812. 
47015. 



-43.77 
-58.70 
-86.76 
-153.63 
-159.41 



0.483 
0.387 
0.155 
0.243 
0.316 



374.1 
384.9 
341.5 
430.2 
652.3 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



Mr,. 

M @ 
M e 
M @ 
Me 



561.5 
699.0 
704.9 
935.2 
1398.0 



0.981 
0.988 
0.963 
0.989 
0.989 



72281. 

67699. 

92003. 
124940. 
219670. 



-542.92 
-246.28 
-570.09 
-469.29 
-517.04 



0.533 
0.166 
1.128 
0.458 
0.986 



549.7 
463.5 
522.3 
787.7 
928.5 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



EP(T 3 ,3.5, 0,-40,0) 



20+20 
20+15 
15+15 
20+10 
15+10 



Me 
M @ 
M @ 
M @ 
M @ 



95.0 
109.5 
126.9 
131.9 
154.1 



0.953 
0.968 
0.949 
0.978 
0.952 



28875. 
37319. 
44601. 
49188. 
60648. 



-1038.40 
-1203.50 
-1160.40 
-1252.90 
-1255.90 



0.998 
1.186 
1.069 
0.999 
1.017 



168.8 
244.7 
322.3 
207.7 
404.3 



M e 
M 
M @ 
M @ 
Me 



175.2 
190.3 
208.4 
263.3 
380.3 



0.973 
0.975 
0.975 
0.969 
0.965 



87007. 

77432. 
102210. 
135110. 
226990. 



-1721.30 
-1045.60 
-1406.10 
-1218.00 
-1027.60 



1.072 
0.648 
0.805 
1.231 
0.960 



348.7 
192.8 
214.6 
548.3 
883.4 



L(2, 0) 



20+20 
20+15 
15+15 
20+10 
15+10 



Me 
M @ 
M Q 
Me 
M Q 



87.0 
99.7 
116.0 
118.0 
140.0 



0.937 
0.953 
0.972 
0.974 
0.985 



18859. 
26088. 
34155. 
38075. 
48463. 



-726.78 
-939.25 
-1087.60 
-1201.00 
-1295.00 



0.997 
1.005 
0.999 
0.990 
0.996 



175.1 
175.0 
189.8 
191.1 
219.6 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



M @ 

M e 
M e 
M @ 
M @ 



148.9 
174.0 
181.6 
235.9 
348.0 



0.987 
0.990 
0.991 
0.991 
0.994 



72221. 

67126. 

89333. 
120130. 
207730. 



-1938.50 
-1420.30 
-1908.10 
-1869.90 
-2077.90 



0.970 
0.986 
0.996 
0.830 
0.709 



209.6 
252.5 
259.5 
274.6 
379.0 



L(2, 1) 



20+20 
20+15 
15+15 
20+10 
15+10 



Me 
M @ 
M @ 
M e 

Mr 



80.3 
92.1 
107.0 
109.4 
129.4 



0.935 
0.960 
0.969 
0.970 
0.969 



33179. 
41065. 
50159. 
55990. 
66431. 



-1379.20 
-1465.70 
-1486.50 
-1663.30 
-1519.40 



0.998 
0.997 
1.003 
0.994 
0.998 



136.5 
152.3 
164.3 
166.4 
187.9 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



M e 

M @ 
M e 
M @ 

Me 



140.1 
160.5 
169.6 
218.8 
321.0 



0.968 
0.969 
0.966 
0.964 
0.932 



99046. 

85317. 
114410. 
146040. 
244970. 



-2345.10 
-1293.90 
-1835.70 
-1373.00 
-1159.90 



0.996 
0.707 
0.673 
0.402 
0.743 



191.6 
167.5 
165.1 
194.6 
404.3 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



HT(2,2) 



20+20 
20+15 
15+15 
20+10 
15+10 



M Q 
Me 
Me 
M Q 
M @ 



389.2 
451.2 
507.1 
536.2 
632.6 



0.964 
0.937 
0.961 
0.960 
0.950 



6138. 
10015. 
12166. 
13624. 
16662. 



1091.40 
1120.00 
1236.10 
1378.00 
1468.50 



-0.539 
0.583 
-1.842 
-0.711 
-1.780 



242.5 
693.5 
322.1 
334.7 
378.9 



M e 
M @ 
M @ 
M @ 

Mr 



733.9 
758.8 
849.2 
1057.3 
1525.5 



0.928 
0.868 
0.905 
0.870 
0.937 



31397. 

34673. 

41087. 
109640. 
214890. 



1977.90 
1301.30 
1898.60 
351.04 
-317.77 



-0.634 
0.951 

-2.966 
0.939 
0.967 



981.8 
783.1 
1192.5 
899.7 
969.4 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



HP(2,2.5) 



20+20 
20+15 
15+15 
20+10 
15+10 



Me 
M @ 
Me 
Me 

M 



403.9 
459.0 
536.6 
530.3 
638.8 



0.923 
0.961 
0.921 
0.958 
0.914 



2544. 
1774. 
3321. 
6259. 
7474. 



1511.00 
1747.20 
1853.80 
1961.60 
2079.90 



0.547 
-1.790 

0.063 
-1.844 
-1.449 



459.7 
279.7 

1155.9 
331.6 

1193.8 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



Mr 

Me 
Me 
M @ 
M Q 



611.2 
816.6 
771.6 
1050.2 
1601.2 



0.918 
0.901 
0.892 
0.850 
0.921 



22867. 

10216. 

27498. 
107210. 
212810. 



2595.60 
2343.10 
2640.30 
707.41 
33.93 



-1.053 
-1.861 
-2.977 
0.893 
0.694 



1200.0 
509.6 

1200.0 
918.8 
916.8 



HT(3,3.5,+2) 



20+20 
20+15 
15+15 
20+10 
15+10 



Me 
Me 
M @ 
M @ 
Me 



358.4 
420.3 
484.0 
495.8 
586.9 



0.977 
0.975 
0.971 
0.983 
0.967 



16787. 
22751. 
29634. 
37522. 
46435. 



81.92 
13.30 
8.43 
-309.28 
-210.60 



0.187 
0.414 
0.016 
0.693 
0.586 



346.4 
728.7 
417.1 
731.8 
916.7 



(20+5 
(10+10 
(15+5 
(10+5 
(5+5 



Mr 
Me 
Me 
Me 
M e 



196.3 
726.0 
303.3 
970.3 
1433.6 



0.983 
0.964 
0.992 
0.992 
0.992 



83529. 

67085. 

98845. 
128810. 
221940. 



-1429.20 
-285.69 

-1096.60 
-755.33 
-679.22 



0.856 
0.594 
0.782 
0.526 
0.400 



232.3 
922.2 
395.9 
909.7 
907.7 



HT(3,3.5, -2) 



20+20 
20+15 
15+15 
20+10 
15+10 



Me 
Me 
Me 
Me 
Me 



316.4 
375.2 
426.3 
436.0 
514.5 



0.981 
0.980 
0.975 
0.986 
0.974 



17922. 
23737. 
31166. 
38125. 
47366. 



-37.78 
-95.75 
-123.77 
-390.93 
-316.76 



0.640 
0.603 
0.506 
0.538 
0.654 



498.6 
618.2 
587.6 
434.1 
806.4 



M e 
Me 
Me 
Me 
Me 



196.1 
639.4 
303.1 
868.1 
1273.2 



0.984 
0.972 
0.993 
0.992 
0.993 



83861. 

68270. 

98715. 
129460. 
223420. 



-1489.90 
-408.04 

-1142.80 
-848.59 
-812.58 



0.869 
0.640 
0.802 
0.675 
0.425 



232.7 
917.6 
389.2 
852.1 
883.0 



(20+!: 
(10+10 
(15+5 
(10+5 
(5+5 



HP(3,3.5,+2) 



20+20 
20+15 
15+15 
20+10 
15+10 



M @ 
Me 
Me 
M @ 
Me 



474.6 
539.6 
634.8 
598.9 
752.5 



0.968 
0.966 
0.955 
0.975 
0.948 



14652. 
20205. 
27087. 
36238. 
45078. 



236.51 
181.76 
170.17 
-213.15 
-109.24 



0.215 
0.071 
0.009 
0.438 
0.539 



863.4 
1076.9 
1200.0 
900.5 
911.2 



M e 
Me 
Me 
Me 
M Q 



196.4 
952.2 
304.1 
1212.8 
1921.0 



0.982 
0.948 
0.990 
0.991 
0.989 



83872. 

66050. 

98220. 
127870. 
220910. 



-1421.20 
-202.66 
-1035.20 
-682.01 
-608.88 



0.928 
0.548 
0.796 
0.555 
0.313 



261.4 
898.9 
405.4 
621.0 
925.7 
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PN model 




FF for projection onto the Fouric 

/end mil 1p -03/2 OtfcAt /cut 


-domain dete 


ction template families 

/o„d mil IpQ '03/2 Ct/ C „ t /cut 


HP(3,3.5, -2) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


363.2 0.973 16421. 113.38 0.384 525.9 
421.8 0.972 21952. 62.83 0.337 690.9 
489.8 0.968 28632. 62.71 0.000 422.7 
510.0 0.982 36893. -272.30 0.263 463.2 

591.3 0.959 45653. -168.47 0.469 924.6 


(2O+5)M 
(1O+1O)M 
(15+5)M 
(1O+5)M 
(5+5)M 


196.3 0.983 83747. -1435.60 0.996 289.7 
734.7 0.958 66819. -271.94 0.680 893.8 
303.7 0.992 98202. -1060.60 0.749 368.8 
998.3 0.991 128060. -722.63 0.491 887.7 
1445.9 0.991 221850. -685.53 0.390 930.8 



TABLE XV: Fitting factors for the projection of the target models (in the rows) onto the (ipo, ^3/2, ct, /cut) Fourier-domain 
detection template family. For ten choices of BBH masses, this table shows the minmax matches between the target (adiabatic) 
models and the Fourier-domain search model, maximized over the intrinsic parameters tpo, ^3/2, an d a,/ cu t, and over the 
extrinsic parameter a. For each intersection, the six numbers shown report the ending frequency / cn d (defined in Sec. IVI Bt of 
the PN model for the BBH masses quoted, the minmax FF mn, and the search parameters at which the maximum is attained. 







FF for projection onto 
mm M rj mn Af r\ 


P(2, 2.5), for 


< i) < 1 
mm M 77 mn M rj 


T(2, 2) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.997 35.53 0.35 0.994 35.55 0.35 
0.997 32.43 0.31 0.994 31.69 0.32 
0.997 28.45 0.29 0.993 27.54 0.32 
0.996 27.72 0.27 0.992 26.83 0.30 
0.995 23.68 0.28 0.988 22.95 0.30 


(2O+5)M 
(10+10)M o 
(15+5)M 
(1O+5)M 
(5+5)M 


0.988 22.04 0.21 0.979 22.51 0.20 
0.996 19.29 0.28 0.990 18.74 0.30 
0.993 18.31 0.23 0.985 17.94 0.24 
0.992 14.55 0.24 0.989 14.86 0.23 
0.994 10.60 0.23 0.992 10.73 0.22 


T(2, 2.5) 


(2O+2O)M 
(2O+15)M 
(15+15)M© 
(2O+1O)M 
(15+1O)M 


0.821 18.77 0.94 0.962 22.59 0.65 
0.862 16.60 0.94 0.966 19.40 0.68 
0.891 14.26 0.97 0.969 16.61 0.71 
0.905 13.67 0.94 0.974 15.95 0.69 
0.929 11.89 0.93 0.978 13.43 0.73 


(2O+5)M 
(10+10)M o 
(15+5)M 
f 10+5)Mn 
(5+5)M 


0.958 11.66 0.63 0.987 12.81 0.53 
0.948 9.96 0.89 0.984 10.71 0.77 
0.965 9.72 0.69 0.989 10.26 0.62 
0.971 8.58 0.60 0.987 8.67 0.59 
0.981 6.55 0.52 0.985 6.52 0.53 


T(3, 3.5, 0) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.997 37.04 0.31 0.994 36.96 0.31 
0.997 32.62 0.29 0.995 32.69 0.29 
0.998 28.89 0.28 0.994 28.07 0.30 
0.997 28.69 0.25 0.995 28.00 0.26 
0.998 24.35 0.26 0.994 23.72 0.27 


(2O+5)M 
(10+10)M o 
(15+5)M 
(1O+5)M 
(5+5)M 


0.996 23.66 0.18 0.990 23.28 0.19 
0.998 19.70 0.26 0.993 19.24 0.27 
0.997 19.18 0.21 0.993 18.82 0.21 
0.997 14.89 0.23 0.993 14.67 0.23 
0.999 10.16 0.24 0.997 10.27 0.24 


T(3, 3.5, 0) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.999 36.20 0.31 0.995 35.37 0.33 
0.999 31.40 0.31 0.997 31.22 0.31 
0.999 27.63 0.29 0.997 26.96 0.31 
0.999 27.42 0.27 0.996 26.84 0.28 
0.998 23.20 0.28 0.996 22.60 0.29 


(2O+5)M 
(10+10)M o 
(15+5)M 
(1O+5)M 
(5+5)M 


0.997 22.95 0.19 0.993 22.53 0.20 
0.999 18.67 0.29 0.995 18.20 0.30 
0.997 18.61 0.22 0.995 18.15 0.23 
0.998 14.13 0.25 0.994 13.87 0.26 
0.998 9.57 0.27 0.996 9.71 0.26 


P(3, 3.5, -2) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.998 35.30 0.33 0.996 34.73 0.34 
0.999 30.84 0.32 0.996 30.65 0.32 
0.999 27.01 0.31 0.996 26.38 0.33 
0.999 26.90 0.28 0.996 26.48 0.29 
0.998 22.76 0.29 0.995 22.18 0.31 


(2O+5)M 
(10+10)M o 
(15+5)M 
(1O+5)M 
(5+5)M 


0.998 23.03 0.19 0.995 22.62 0.20 
0.998 18.11 0.31 0.994 17.86 0.31 
0.997 18.43 0.22 0.994 17.96 0.23 
0.998 13.99 0.25 0.993 13.75 0.26 
0.997 9.50 0.27 0.996 9.63 0.27 


P(3, 3.5, +2) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.999 33.58 0.36 0.996 33.42 0.37 
0.999 30.03 0.33 0.997 29.70 0.34 
0.998 26.12 0.33 0.997 25.59 0.34 
0.999 26.38 0.29 0.997 25.84 0.30 
0.997 21.62 0.32 0.995 21.53 0.32 


(2O+5)M 
(10+10)M o 
(15+5)M 
(1O+5)M 
(5+5)M 


0.998 22.71 0.19 0.996 22.42 0.20 
0.999 17.87 0.31 0.995 17.36 0.33 
0.998 18.15 0.23 0.996 17.77 0.24 
0.998 13.59 0.27 0.994 13.31 0.28 
0.998 9.25 0.29 0.996 9.34 0.28 


ET(2, 2.5) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.994 26.75 0.56 0.989 25.10 0.65 
0.993 23.77 0.52 0.962 25.26 0.45 
0.991 20.87 0.50 0.970 21.86 0.45 
0.988 21.35 0.42 0.973 20.26 0.47 
0.987 17.99 0.44 0.969 17.28 0.48 


(20+5) M Q 
(10+10)M o 
(15+5)M 
(1O+5)M 
(5+5)M 


0.979 19.87 0.24 0.970 19.27 0.26 
0.989 14.75 0.43 0.983 14.93 0.43 
0.987 15.81 0.28 0.982 15.42 0.30 
0.994 11.98 0.33 0.990 11.70 0.34 
0.997 8.04 0.36 0.995 8.18 0.35 


ET(3, 3.5, -2) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.991 31.38 0.46 0.986 29.96 0.53 
0.989 28.48 0.40 0.978 26.83 0.47 
0.985 25.24 0.38 0.970 23.83 0.43 
0.977 25.09 0.34 0.955 23.62 0.39 
0.974 21.65 0.33 0.963 22.39 0.31 


(2O+5)M 
(10+10)M o 
(15+5)M 
(10+5)M o 
(5+5)M 


0.952 22.75 0.20 0.941 23.33 0.19 
0.977 18.69 0.29 0.971 18.03 0.32 
0.972 18.53 0.22 0.964 19.03 0.21 
0.983 15.04 0.22 0.979 14.79 0.23 
0.994 10.38 0.23 0.992 10.39 0.23 


ET(3, 3.5, +2) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.993 30.84 0.46 0.989 29.51 0.51 
0.991 27.38 0.43 0.981 25.94 0.48 
0.986 24.02 0.41 0.973 22.67 0.47 
0.981 24.19 0.36 0.963 22.66 0.42 
0.977 20.84 0.35 0.966 21.46 0.33 


(2O+5)M 
(10+10)M o 
(15+5)M 
(1O+5)M 
(5+5)M 


0.957 22.28 0.20 0.946 22.87 0.19 
0.983 16.95 0.35 0.976 17.42 0.33 
0.974 17.92 0.24 0.967 18.49 0.22 
0.984 14.43 0.24 0.982 14.28 0.24 
0.995 9.80 0.26 0.993 9.89 0.25 


EP(2, 2.5) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.988 30.91 0.48 0.977 28.86 0.58 
0.980 27.79 0.43 0.963 25.85 0.52 
0.972 24.47 0.40 0.947 22.77 0.49 
0.965 24.97 0.34 0.938 22.29 0.47 
0.963 23.00 0.29 0.951 21.93 0.32 


(2O+5)M 
(1O+1O)M 
(15+5)M 
(10+5)M o 
(5+5)M 


0.947 24.15 0.17 0.940 23.60 0.18 
0.975 18.50 0.30 0.964 17.90 0.32 
0.970 18.73 0.22 0.963 19.16 0.21 
0.984 15.15 0.22 0.980 14.80 0.23 
0.995 10.24 0.24 0.993 10.29 0.24 


EP(3, 3.5, -2) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.993 30.25 0.48 0.990 29.04 0.53 
0.990 26.86 0.45 0.981 25.54 0.50 
0.986 23.98 0.41 0.974 22.36 0.48 
0.982 23.79 0.37 0.964 22.56 0.42 
0.977 20.49 0.36 0.966 21.21 0.34 


(2O+5)M 
(1O+1O)M 
(15+5)M 
(10+5)M o 
(5+5)M 


0.958 21.90 0.21 0.947 22.61 0.20 
0.983 16.74 0.36 0.976 17.26 0.34 
0.975 17.83 0.24 0.967 18.24 0.23 
0.984 14.34 0.24 0.982 14.12 0.25 
0.994 9.74 0.26 0.993 9.84 0.26 


EP(3, 3.5, +2) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.994 29.47 0.50 0.991 28.39 0.55 
0.991 26.46 0.45 0.983 24.97 0.52 
0.986 22.97 0.44 0.975 21.73 0.50 
0.982 23.18 0.39 0.966 22.14 0.43 
0.978 19.94 0.38 0.968 20.58 0.36 


(2O+5)M 
(10+10)M o 
(15+5)M 
(1O+5)M 
(5+5)M 


0.960 21.84 0.21 0.948 22.30 0.20 
0.983 16.14 0.39 0.976 16.75 0.36 
0.977 17.52 0.24 0.968 18.08 0.23 
0.985 13.53 0.27 0.983 13.79 0.26 
0.994 9.54 0.27 0.993 9.55 0.27 
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FF for projection onto 
mm M rj mn M r/ 


P(2, 2.5), for 


< )) < 1 
mm M rj mn M rj 


HT(3, 3.5, -2) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.993 21.45 0.98 0.991 21.03 1.00 
0.986 19.86 0.84 0.982 18.48 1.00 
0.978 17.27 0.81 0.972 16.19 0.94 
0.965 20.87 0.49 0.949 18.74 0.66 
0.952 17.74 0.49 0.944 16.36 0.59 


(2O+5)M 
(10+10)M o 
(15+5)M Q 
(10+5)M o 
(5+5)M 


0.995 26.36 0.15 0.986 25.79 0.15 
0.964 15.24 0.43 0.958 14.57 0.48 
0.988 19.17 0.21 0.980 19.60 0.20 
0.978 14.07 0.25 0.975 13.93 0.26 
0.987 9.61 0.27 0.986 9.55 0.27 


HT(3, 3.5, +2) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.982 20.21 1.00 0.960 20.04 1.00 
0.984 17.81 0.98 0.967 17.53 1.00 
0.977 15.20 1.00 0.962 16.06 0.89 
0.964 19.18 0.57 0.950 17.04 0.77 
0.954 15.66 0.61 0.943 16.70 0.54 


(20+5) M 
(1O+1O)M 
(15+5)M Q 
(10+5)M o 
(5+5)M 


0.997 25.94 0.15 0.990 25.48 0.16 
0.965 13.39 0.55 0.959 13.95 0.51 
0.991 18.63 0.22 0.984 19.08 0.21 
0.980 13.23 0.28 0.975 13.51 0.27 
0.986 9.03 0.30 0.985 8.93 0.31 


HP(3, 3.5, -2) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.962 19.87 1.00 0.946 20.16 1.00 
0.971 17.46 1.00 0.960 17.69 1.00 
0.963 15.66 0.93 0.960 15.14 1.00 
0.961 17.81 0.69 0.950 19.45 0.55 
0.947 16.31 0.56 0.941 15.77 0.60 


(2O+5)M 
(1O+1O)M 
(15+5)M 
(10+5)M o 
(5+5)M 


0.997 25.87 0.15 0.990 25.26 0.16 
0.962 12.92 0.59 0.957 13.34 0.55 
0.992 18.51 0.22 0.982 17.98 0.23 
0.978 13.04 0.29 0.975 13.27 0.28 
0.984 8.97 0.30 0.982 9.02 0.30 


HP(3, 3.5, +2) 


(2O+2O)M 
(2O+15)M 
(15+15) Af 
(2O+1O)M 
(15+1O)M 


0.915 19.33 1.00 0.887 20.18 0.94 
0.942 17.26 1.00 0.921 17.71 0.96 
0.938 15.03 0.99 0.933 14.89 1.00 
0.959 16.40 0.81 0.947 17.94 0.65 
0.949 12.43 0.99 0.937 13.30 0.86 


(2O+5)M 
(10+10)M o 
(15+5)M 
(10+5)M o 
(5+5)M 


0.998 25.69 0.15 0.992 25.21 0.16 
0.957 12.04 0.67 0.953 11.32 0.76 
0.993 18.25 0.23 0.985 18.61 0.22 
0.978 12.90 0.29 0.976 12.65 0.31 
0.982 8.62 0.33 0.982 8.70 0.32 



TABLE XVI: Fitting factors for the projection of the target models (in the rows) onto the P(2, 2.5) detection template family. 
For ten choices of BBH masses, this table shows the maxmax (mm) and minmax (mn) matches between the target models 
and the P(2,2.5) search model, maximized over the intrinsic parameters of the search model. For each intersection, the triples 
(mm,M,?7) and (mn,M,r7) denote the maximized matches and the mass parameters M = mi + m-z and n — mimz/M 2 at 
which the maxima are attained (maxmax and minmax matches give rise to slightly different optimal values of M and rj). In 
computing these matches, the search parameter r\ was not restricted to its physical range < r\ < 1/4, but it was allowed to 
move in the range < rj < 1, for which the energy-balance equation 1311 is still formally integrable. With few exceptions, this 
table shows that maxmax and minmax matches are very similar, so we generally use the more conservative minmax matches. 







FF for projection onto E 
mm M rj mmc M rj f cu t 


P(3, 3.5,0), fo 


rO<i|<l 
mm M rj mmc M rj / cu t 


T(2, 2) 


(2O+2O)Af 
(2O+15)M 
(15+15) AT© 
(2O+1O)M 
(15+10)M© 


0.984 51.05 0.14 0.984 51.38 0.14 171.7 
0.981 44.12 0.14 0.981 44.11 0.15 199.7 
0.977 37.15 0.16 0.978 37.16 0.16 236.7 
0.974 35.85 0.15 0.974 35.62 0.16 246.7 
0.976 27.95 0.20 0.976 27.82 0.20 316.6 


(2O+5)M 
(10+10)M© 
(15+5)M© 
(10+5)M© 
(5+5)M 


0.981 25.34 0.16 0.981 25.32 0.16 347.1 
0.984 22.16 0.21 0.985 22.15 0.21 395.5 
0.986 20.90 0.18 0.987 20.89 0.18 424.3 
0.992 16.17 0.20 0.999 16.20 0.20 368.4 
0.996 11.05 0.21 0.999 11.12 0.21 553.1 


T(2, 2.5) 


(2O+2O)M 
(20+15) M Q 
(15+15) Af s 
(2O+1O)M 
(15+10)M© 


0.948 24.94 0.57 0.985 24.09 0.60 202.5 
0.956 21.50 0.58 0.990 20.77 0.62 241.9 
0.965 18.08 0.62 0.986 18.76 0.57 279.0 
0.965 17.57 0.58 0.992 17.00 0.62 285.7 
0.974 14.85 0.61 0.994 14.28 0.66 329.9 


(2O+5)M 
(10+10)M© 
(15+5)M© 
(10+5)M© 
(5+5)M 


0.975 14.60 0.41 0.975 14.52 0.42 567.2 
0.983 11.72 0.65 0.995 11.48 0.69 415.0 
0.983 11.38 0.51 0.994 11.50 0.50 448.1 
0.986 9.44 0.50 0.993 9.23 0.53 629.2 
0.989 6.89 0.48 0.990 6.93 0.47 787.5 


T(3, 3.5, 0) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.979 53.09 0.12 0.979 52.83 0.12 166.4 
0.971 45.36 0.13 0.972 45.28 0.13 194.2 
0.969 37.79 0.15 0.969 37.79 0.15 232.8 
0.961 34.57 0.17 0.963 34.57 0.16 254.3 
0.971 28.74 0.18 0.971 28.63 0.18 306.5 


(2O+5)M 
(10+10)M© 
(15+5)M© 
(10+5)M© 
(5+5)M 


0.965 27.22 0.13 0.966 27.27 0.13 322.7 
0.979 22.77 0.19 0.979 22.79 0.19 384.8 
0.976 21.94 0.16 0.978 22.09 0.15 398.7 
0.985 16.47 0.19 0.985 16.47 0.19 533.1 
0.994 10.69 0.22 0.999 10.66 0.22 474.6 


T(3, 3.5, 0) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.948 52.47 0.11 0.979 50.43 0.13 174.3 
0.967 43.35 0.14 0.968 43.41 0.14 202.7 
0.963 33.41 0.20 0.966 36.09 0.16 243.5 
0.963 33.33 0.17 0.964 33.17 0.18 267.0 
0.971 27.25 0.20 0.972 27.25 0.20 321.4 
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0.977 21.88 0.21 0.980 21.59 0.21 408.5 
0.973 19.98 0.19 0.976 21.32 0.16 411.8 
0.985 15.21 0.22 0.998 15.38 0.21 346.9 
0.994 10.14 0.24 0.999 10.13 0.24 522.3 


P(2, 2.5) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.941 55.20 0.11 0.956 56.36 0.10 152.1 
0.940 44.34 0.14 0.940 44.34 0.14 198.2 
0.946 37.08 0.16 0.948 37.27 0.15 236.2 
0.943 37.06 0.14 0.943 37.13 0.14 236.9 
0.945 30.16 0.16 0.948 30.40 0.16 289.2 


(20+5) Mq 
(10+10)M© 
(15+5)M© 
(10+5)M© 
(5+5)M 


0.937 27.97 0.13 0.938 27.93 0.13 315.3 
0.958 22.26 0.20 0.958 22.20 0.20 395.3 
0.959 22.06 0.15 0.961 22.26 0.15 395.1 
0.977 15.53 0.21 0.998 15.83 0.20 284.1 
0.992 10.42 0.23 0.999 10.37 0.23 408.0 


P(3, 3.5, -2) 


(2O+2O)M 
(20+15) M e 
(15+15)M 
(2O+1O)M 
(15+1O)M 


0.979 49.53 0.14 0.979 49.58 0.14 179.2 
0.972 42.49 0.15 0.972 42.49 0.15 206.8 
0.962 33.09 0.20 0.970 35.17 0.17 250.1 
0.965 32.49 0.18 0.964 32.67 0.18 273.4 
0.974 26.67 0.21 0.974 26.65 0.21 328.4 


(20+5)M© 
(10+10)M© 
(15+5)M© 
(10+5)M© 
(5+5)M 


0.955 25.30 0.16 0.959 26.63 0.14 330.9 
0.982 21.04 0.22 0.982 21.04 0.22 416.3 
0.977 21.04 0.17 0.978 21.07 0.17 416.6 
0.987 15.39 0.21 0.998 15.27 0.21 356.4 
0.996 10.10 0.25 0.996 10.11 0.25 795.8 


P(3, 3.5, +2) 


(20+20)M o 
(2O+15)M 
(15+15)M© 
(2O+1O)M 
(15+1O)M 


0.976 49.06 0.14 0.980 48.36 0.15 181.7 
0.972 41.31 0.16 0.974 41.27 0.16 213.0 
0.971 33.77 0.19 0.971 33.77 0.19 259.5 
0.966 31.61 0.19 0.966 31.66 0.19 277.0 
0.975 25.80 0.22 0.975 25.81 0.22 338.4 


(20+5)M© 
(10+10)M© 
(15+5)M© 
(10+5)M© 
(5+5)M 


0.956 26.37 0.14 0.957 26.37 0.14 333.4 
0.983 20.41 0.24 0.983 20.41 0.24 429.2 
0.976 20.81 0.17 0.976 20.72 0.17 423.6 
0.988 15.07 0.22 0.988 15.07 0.22 580.0 
0.996 9.84 0.26 1.000 9.81 0.26 566.3 


ET(2, 2.5) 


(2O+2O)M 
(2O+15)M 
(15+15)M 
(2O+1O)M 
(15+10)M© 


0.998 35.41 0.31 0.999 35.10 0.31 244.8 
0.999 30.78 0.30 0.999 30.78 0.30 280.3 
0.998 26.43 0.31 0.998 26.53 0.30 324.7 
0.998 27.01 0.26 0.998 27.04 0.26 324.1 
0.998 22.16 0.29 0.998 22.17 0.29 393.8 
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(15+5)M© 
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0.995 22.57 0.19 0.995 22.62 0.18 392.5 
0.999 17.42 0.32 0.999 17.42 0.32 492.4 
0.996 17.89 0.22 0.996 17.93 0.22 488.0 
0.997 13.19 0.27 0.997 13.09 0.28 658.5 
0.999 8.61 0.32 0.999 8.60 0.32 996.6 
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0.941 
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31.21 
17.14 
21.56 
15.05 
9.81 



0.10 
0.35 
0.16 
0.22 
0.26 



0.962 
0.982 
0.988 
0.993 
0.992 



30.26 
17.25 
21.56 
15.02 
9.83 



0.11 
0.34 
0.16 
0.22 
0.26 



287.2 
493.8 
407.3 
582.1 
798.6 



0.932 
0.973 
0.987 
0.992 
0.989 



30.88 
14.84 
21.20 
14.15 

9.23 



0.10 
0.45 
0.17 
0.25 
0.29 
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0.974 
0.987 
0.992 
0.989 



29.95 
14.74 
21.15 
14.15 
9.28 



0.11 
0.46 
0.17 
0.25 
0.29 



292.5 
553.2 
416.8 
615.7 
754.7 



0.930 
0.970 
0.986 
0.991 
0.987 



30.83 
15.51 
21.04 
13.83 
9.20 



0.10 
0.41 
0.17 
0.26 
0.29 



0.953 
0.970 
0.987 
0.991 
0.987 



29.87 
15.54 
20.98 
13.85 
9.23 



0.11 
0.41 
0.17 
0.26 
0.29 



294.6 
531.0 
420.3 
633.8 
640.9 



0.951 
0.964 
0.986 
0.988 
0.985 



29.72 
14.15 
20.84 
13.44 



0.11 
0.49 
0.17 
0.27 
0.31 



0.951 
0.964 
0.987 
0.988 
0.986 



29.70 
14.28 
20.76 
13.65 
8.99 



0.11 
0.48 
0.17 
0.26 
0.30 



296.4 
559.8 
423.1 
634.4 
724.4 



TABLE XVII: Fitting factors for the projection of the target models (in the rows) onto the EP(3, 3.5, 0) detection template 
family. For ten choices of BBH masses, this table shows the maxmax matches between the target models and the EP(3, 3.5, 0) 
search model, with (mmc) and without (mm) the time-domain cut discussed in Sec. IVIll The matches are maximized over the 
intrinsic parameters of the search model (over M and rj for the mm values; over M, rj and / cu t for the mmc values). For each 
intersection, the triple (mm,M,?)) and the quadruple (mm,M,?7,/cut) denote the maximized matches and the mass (and cut) 
parameters at which the maxima are attained. In computing these matches, the search parameter rj was not restricted to its 
physical range < rj < 1/4, but it was allowed to move in the range < rj < 1 for which the energy-balance equation 1)3 l|l is 
still formally integrable. This table shows that the addition of the time-domain cut can improve the fitting factors considerably, 
especially for the higher M's in the in the left half of the table, and for the models whose orbital evolution is ended within the 
range of good interferometer sensitivity. 
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